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Abstract 

An intriguing new duality between planar MHV gluon amplitudes and light-like Wilson 
loops in AT = 4 super Yang-Mills is investigated. We extend previous checks of the duality by 
performing a two-loop calculation of the rectangular and pentagonal Wilson loop. Further- 
more, we derive an all-order broken conformal Ward identity for the Wilson loops and analyse 
its consequences. Starting from six points, the Ward identity allows for an arbitrary function 
of conformal invariants to appear in the expression for the Wilson loop. We compute this 
function at six points and two loops and discuss its implications for the corresponding gluon 
amplitude. It is found that the duality disagrees with a conjecture for the gluon amplitudes 
by Bern et al. A recent calculation by Bern et al indeed shows that the latter conjecture 
breaks down at six gluons and at two loops. By doing a numerical comparison with their 
results we find that the duality between gluon amplitudes and Wilson loops is preserved. 
This review is based on the author's PhD thesis and includes developments until May 2008. 
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1 Introduction 



In the standard model of elementary particles, the strong interactions are described by quantum 
chromodynamics (QCD). It is a non-Abelian Yang-Mills gauge theory, with the quarks being 
in the fundamental representation of the gauge group SU(3). In contrast to quantum electro- 
dynamics, the non-Abelian nature of the gauge group allows the gauge bosons to interact with 
each other. It is this property which leads to asymptotic freedom and the confinement of quarks. 
Despite the simplicity and elegance of the QCD Lagrangian, many open problems remain, such 
as understanding the transition from the high-energy region (short distances), where perturba- 
tion theory is valid, to the low-energy region (long distances) of confined quarks. 

't Hooft proposed to consider Yang- Mills theories with gauge group SU(N), for N very large, 
while keeping the 't Hooft coupling g 2 N (with g being the Yang-Mills coupling constant) fixed 
PQ. In this way, one may hope to see simplifications for N large, and eventually get insight 
into QCD by computing terms in an 1/N expansion. The latter resembles that of the genus 
expansion of string theory. For N — > oo, non-planar diagrams are suppressed, which is why 
the 't Hooft limit is also referred to as the planar limit. Hints that simplifications may occur 
in the 't Hooft limit appeared in [21 [3], where integrable structures were found when studying 
high-energy scattering in gauge theory. Nevertheless, QCD in the 't Hooft limit is still far from 
being solved. It seems natural to study supersymmetric Yang-Mills theories, which are much 
simpler. In particular, the maximally supersymmetric Yang-Mills theory, N = 4 SYM, has 
many special properties. It was shown long ago that its (3 function vanishes, and hence its cou- 
pling constant does not run. It is an interacting superconformal field theory, where the coupling 
constant is a free parameter. Moreover, through the AdS/CFT correspondence it is expected to 
be dual to type IIB superstring theory on AdS$ x 55 [3J. This is a realisation of 't Hooft's idea 
of a gauge/string duality in the N — ► 00 limit. The nature of this duality, which relates field 
theory at strong coupling with string theory at weak coupling, implies that the perturbative 
series of 'observables' (such as e.g. correlation functions of gauge-invariant operators) in N = 4 
SYM has to reproduce perturbative string theory results. In order for this to happen, one may 
suspect that the perturbative expansion of these quantities must have special properties. In- 
deed, remarkably simple structures have been observed in N = 4 SYM, mainly in two domains: 
firstly, anomalous dimensions of composite operators and secondly, on-shell n-particle scattering 
amplitudes. 

The AdS/CFT correspondence identifies string states with composite, gauge-invariant oper- 
ators on the gauge theory side. In testing the correspondence, computing scaling dimensions of 
such operators in N = 4 SYM plays an important role. The first tests were done for 'protected' 
operators, whose scaling dimension can be shown to be coupling independent. It therefore equals 
their classical (tree-level) dimension, which has to agree with the string theory prediction. For 
generic operators, the scaling dimension depends on the coupling constant, and a comparison 
with string theory is difficult because of the weak/strong nature of the duality, i.e. it would 
require summing the complete perturbation series. 

Important progress in computing anomalous dimensions in M = 4 SYM was made with the 
discovery of integrability in the planar (large N) limit. In QCD, integrable structures at one 
loop were observed some time ago [HI [6], El El 13 EH]- In certain cases, the dilatation operator, 
which measures the conformal dimension of a given operator, turns out to be described by an 
integrable spin chain. In J\f = 4 SYM, the integrability of the complete dilatation operator at 
one loop was found in [T2| [T3"l [T3] . Most importantly, in some sectors it was shown to extend 
to higher loop levels PSI USl HZl H51 US1 12D1 12X1 1221 1231 125]- If this integrability persists to 
arbitrary loop levels, then one can hope to make contact with strong coupling results. Indeed, 
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on the string theory side of the AdS/CFT correspondence, the classical Green-Schwarz super- 
string action for AdS§ x S$ , constructed in , is integrable [27] . If the integrability survives 
quantisation [281 [29], it might be possible to identify the integrable structures on both sides of 
the correspondence. 

Integrability, or in certain cases the assumption of integrability, allowed the computation of 
the anomalous dimensions of various operators in M = 4 SYM. Spectacular progress was made 
by the all-order conjecture of Beisert, Eden and Staudacher [30] (see also [3T] [32] [331 [MIES]) for 
the asymptotic Bethe ansatz H describing the anomalous dimensions of composite operators in 
the SL(2) sector of M = 4 SYM. Their proposal leads to an integral equation (BES equation) 
for the cusp anomalous dimension [36] E7J [38] , valid to all orders in the coupling constant. It cor- 
rectly reproduces the known perturbative values of the cusp anomalous dimension [38 [ [TBI 139] 130] 
and also agrees with results at strong coupling [UJ [321 [33] (for the strong-coupling expansion of 
the BES equation see [3311351136]). 

The second domain where unexpected simplicity was found is that of on-shell gluon scatter- 
ing. Indeed, the scattering amplitudes turn out to be much simpler than one would expect on 
general grounds. Even at tree-level, the number of Feynman diagrams contributing to a given 
amplitude increases factorially with the number of external gluons. Nevertheless, some classes 
of tree-level amplitudes are known for an arbitrary number of legs. For example, the maximally 
helicity- violating (MHV) amplitudes are given by very simple one-line expressions [471 135] , and 
all tree-level amplitudes satisfy recursion relations derived by Britto, Cachazo and Feng (and 
Witten) [491 150] . Moreover, tree-level amplitudes have simple properties in twistor space [51] . 

Furthermore, enormous progress has been achieved to compute loop level amplitudes, mainly 
using unitarity-based techniques [52J [53] , which are inspired by the Cutkosky rules [53] . These 
have allowed the computation of a large class of one-loop amplitudes, and in the case of N = 4 
SYM even some two-, three-, and four-loop amplitudes [55, 56| 1391 [30] . Anastasiou, Bern, Dixon 
and Kosower (ABDK) noticed that the four-gluon amplitude at two loops has a remarkable 
iterative structure [56]. For the infrared divergent terms of the amplitude, such an iteration is 
expected on general grounds, but a similar iteration was found to hold also for the finite part. 
Following years of intensive studies of gluon scattering amplitudes [521 1551 157] and based on the 
observation of ABDK, the conjecture was put forward by Bern, Dixon and Smirnov [39] that the 
maximally helicity-violating (MHV) planar gluon amplitudes in M = 4 SYM have a remarkably 
simple all-loop iterative structure. In general, these amplitudes have the following form: 

In A< MHV) = PR divergences] + 4 MHV) (pi, ...,p n ;a) + 0(e) . (1) 

Here A^1 MHV ^ is the colour-ordered planar gluon amplitude, divided by the tree amplitude. 
The first term on the right-hand side describes the infrared (IR) divergences and the second 
term is the finite contribution dependent on the gluon momenta pi and on the 't Hooft cou- 
pling a = g 2 N/(8ir 2 ). The structure of IR divergences is well understood in any gauge the- 
ory ^W\Wi^MiMi^Mi\MM\\M\M\M^\Mn^. In particular, in theories with 
a vanishing (3 function like M = 4 SYM, the leading IR singularity in dimensional regularisa- 
tion is a double pole, whose coefficient is the universal cusp anomalous dimension appearing 
in many physical processes [381 EH [751 02] ■ Interestingly, the latter is also predicted by inte- 
grable models, as mentioned above. The BDS conjecture provides an explicit expression for 
the finite part, i^ MHV ^ = F^ DS \ for an arbitrary number n of external gluons, to all orders 
in the coupling a. Remarkably, the dependence of F^ DS ^ on the kinematical invariants is de- 
scribed by a function which is coupling independent and, therefore, can be determined at one 

2 Asymptotic means that it is valid only for 'long' operators, i.e. valid to 0(g 2L ~ 2 ), where L is the number of 
elementary fields constituting the composite operator. 
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loop. At present, the BDS conjecture has been tested up to three loops for n = 4 [39] and up 
to two loops for n = 5 [76] . The results of this thesis are relevant for the case n = 6 at two loops. 

The aforementioned scattering amplitudes in M = 4 have even more interesting features. 
Drummond, Sokatchev, Smirnov, and the present author found that all integrals appearing in 
the four-gluon amplitude up to three loops are of a special type [77J. In dual momentum variables 
defined by 

Pi = x i+1 - Xi , (2) 

all integrals can be seen to have broken conformal properties. This is unexpected since this 
broken dual conformal symmetry is not, at least not in an obvious way, a consequence of the 
usual (super-)conformal symmetry of N = 4 SYM. This surprising observation was confirmed 
at four loops [40] and was used as an assumption for constructing the four-gluon amplitude at 
five loops [75]. 

In an important recent development in the study of the AdS/CFT correspondence, Alday and 
Maldacena proposed [79J the strong coupling description of planar gluon scattering amplitudes 
in M = 4 SYM and were able to make a direct comparison with the BDS prediction based 
on weak coupling results for the same amplitudes. According to their proposal, certain planar 
gluon amplitudes at strong coupling are related to the area of a minimal surface in AdSs space 
attached to a specific closed contour C n , 



lnA4 n = A min (C n ) . (3) 

The contour C n is a polygon with light-like edges defined by the gluon momenta 

through relation ([2]), and with the cyclicly condition x n+ \ = x\. Notice that the Xj, which 
correspond to the n cusps of the polygon C n , are the same dual momentum variables that were 
introduced previously for discussing the broken conformal properties of integrals appearing in 
the gluon amplitudes! The prescription of [79] is insensitive to the helicity configuration of the 
gluon amplitude under consideration. 

For n = 4 the minimal surface j4 m in(C4) was found explicitly in [79], by making use of the 
conformal symmetry of the problem. With the appropriate AdS equivalent of dimensional reg- 
ularisation, the divergent part of In A^4 has the expected pole structure, with the coefficient in 
front of the double pole given by the known strong coupling value of the cusp anomalous dimen- 
sion. Most importantly, the finite part of ln.M4 is in perfect agreement with i^ BDS ) from the 
BDS ansatz. For n > 5 the practical evaluation of the solution of the classical string equations 
turns out to be difficult, but it becomes possible for n large [50]. In the limit n — > oo the strong 
coupling prediction for lnA4 n disagrees with the BDS ansatz. This indicates [80] that the BDS 
conjecture should fail for a sufficiently large number of gluons and/or at sufficiently high loop 
level. 



Alday and Maldacena also pointed out [79] that their prescription ([3]) is mathematically 
equivalent to the strong coupling calculation of the expectation value of a Wilson loop W(C n ), 
defined on the light-like contour C n [81, 82], 

W{C n ) = -i(0| TrPexp (ig £ dx^A^j |0> . (4) 

This should not come as a total surprise, since the intimate relationship between the infrared 
divergences of the scattering of massless particles and the ultraviolet divergences of Wilson loops 
with cusps is well known in QCD |38[ [741 175) . Inspired by this, Drummond, Korchemsky and 
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Sokatchev conjectured that a similar duality relation between planar MHV gluon amplitudes 
and light-like Wilson loops also exists at weak coupling [83J. They illustrated this duality by 
an explicit one-loop calculation in the simplest case n = 4. This was later extended to the 
case of arbitrary n at one loop by Brandhuber, Heslop and Travaglini |84J. The duality relation 
identifies the two objects up to an additive constant and non-planar corrections, 

In _M^ MHV) = In W{C n ) + const + 0(e, 1/N) . (5) 

This means that upon a specific identification of the regularisation parameters and the kine- 
matical invariants, the infrared divergences of the logarithm of the scattering amplitude In Ai n , 
match the ultraviolet divergences of the light-like Wilson loop In W(C n ), and, most importantly, 
the finite parts of the two objects also coincide (up to an additive constant and non-planar cor- 
rGctions) 

f (muv) = F (WL) + congt + Q(l/N) . (6) 

While the former property follows from the known structure of divergences of scattering ampli- 
tudes and of Wilson loops in generic gauge theories [381 ESI El] , the property © is extremely 
non-trivial. 

In this report we present work in collaboration with Drummond, Korchemsky and Sokatchev 
that provided further evidence in favour of the duality relation ©. To begin with, we carried 
out explicit two-loop calculations of \nW{C n ) for n = 4 [86J and n = 5 |87j. Our results are 
in perfect agreement with the two-loop MHV gluon amplitude calculations [561 ES] , and hence 
with the BDS ansatz for n = 4, 5. 

Furthermore, in [SB] we argued that we can profit from the (broken) conformal symmetry of 
the light-like Wilson loops. Due to the presence of a cusp anomaly in the Wilson loops, conformal 
invariance manifests itself in the form of anomalous Ward identities. In [86J we proposed and 
in [57] we proved an anomalous conformal Ward identity for the finite part of the Wilson loops, 
valid to all orders in the coupling constant H. It reads 



X i,i+2 n 



< M , (7) 



J2 ""i.i+l 

where and r cusp (a) is the (coupling-dependent) cusp anomalous dimension. This 

identity uniquely fixes the functional form of the finite part of lnW(C n ) for n = 4 and n = 5, 
up to an additive constant, to agree with the conjectured BDS form for the corresponding MHV 
gluon amplitudes. For n > 6, ([7|) gives partial restrictions on the functional dependence on the 
kinematical variables. Quite remarkably, the BDS ansatz ([1]) for the n-gluon MHV amplitudes 
satisfies the conformal Ward identity for arbitrary n [86j. 

However, for n > 6 the conformal Ward identity allows F^ L ^ to differ from the BDS ansatz 
^t(bds) ^ an arbitrary function of conformal invariants (for n = 6 there are three such in- 
variants) [86J. This result provided a possible explanation of the BDS conjecture for n = 4, 5 
(assuming that the MHV amplitudes have the same conformal properties as the Wilson loop), 
but left the door open for potential deviations from it for n > 6. To verify whether the BDS 
conjecture and/or the proposed duality relation ([6]) still hold for n = 6 to two loops, it was 
necessary to perform explicit two-loop calculations of the finite parts of the six-gluon MHV 
amplitude i*g MH , and of the hexagon Wilson loop F^ L \ 



3 Later, similar Ward identities were also obtained at strong coupling using the AdS/CFT correspondence in 

Refs. IHH1IH3- 
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By doing an explicit two-loop calculation we were able to derive a (multiple) parameter in- 
tegral representation for i ? g WL \ which we evaluated numerically. We found that -Fg WL ^ differs 
from the BDS ansatz at two loops by a non-trivial function of conformal cross-ratios. We stud- 
ied this function numerically and found that it was consistent with the collinear limit (of gluon 
amplitudes) [90]. Later, when the results of the two-loop calculation of the six- gluon MHV am- 
plitude became available, we compared the two results numerically and found agreement with 
the duality relation ([6]) [HJ E2]. The BDS ansatz fails at two loops and at n = 6, but it fails 
just in such a way as to preserve the duality with Wilson loops! We take this as strong evidence 
that the duality relation ([6]) holds for arbitrary n and to an arbitrary number of loops. 

We would like to point out that a weaker form of the duality © has already been observed 
in QCD in the special, high-energy (Regge) limit s — t > for the four-gluon amplitude up 
to two loops [75]. The same relationship holds in any gauge theory ranging from QCD to N = 4 
SYM. The essential difference between these theories is that in the former case the duality is only 
valid in the Regge limit, whereas in the latter case it is exact in general kinematics. We would 
like to mention that a thorough analysis of the Regge limit of planar multi-gluon amplitudes in 
N = 4 SYM was recently performed in [93l US] and it provided further evidence that the BDS 
ansatz needs to be corrected [94J. 
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This report is organised as follows: 

Sections [2H6] are mainly introductory. In section El we recall some notions of conformal sym- 
metry and its implications for correlation functions of gauge invariant operators in M = 4 SYM. 
section [3] briefly introduces the reader to gluon scattering amplitudes, focusing on multi-loop 
results for maximally helicity-violating (MHV) amplitudes and the BDS conjecture. Following 
that, section [J] summarises the prescription of Alday and Maldacena to compute gluon scattering 
amplitudes at strong coupling. This will lead us to consider Wilson loops, and after introducing 
them in section [5] and summarising their renormalisation properties, we discuss in section [6] the 
main point of study of this thesis: the duality between gluon amplitudes and Wilson loops. The 
remaining sections present original work by the author and collaborators. In particular sections 
ESlandOare based on [77], and sections U\B E and appendix[B]are based on [861 IH71I9U1 M\ [95] . 

The only exceptions to the introductory nature of sections [2H6] are sections 12.61 and 1 3 . 21 where 
we present original work on off-shell conformal four-point integrals, and their relevance for on- 
shell gluon scattering amplitudes, respectively. We study the known (to four loops) four-gluon 
scattering amplitude in N = 4 SYM and find that it is given in terms of 'pseudo-conformal' 
integrals. 

In section [7J we verify the Wilson loop /gluon amplitude duality in the non-trivial cases of 
four and five points/gluons at two loops by an explicit Feynman graph calculation of the rect- 
angular and pentagonal Wilson loops. 

In section [HJ we derive all-order Ward identities for the light-like Wilson loops and discuss 
their consequences. 

From this discussion of the Ward identities in section [8] it will become clear that a crucial test 
of the duality is at two loops and for six points/gluons. We therefore compute the hexagonal 
Wilson loop at two loops in section [9] and study its properties in the collinear limit (of gluon 
amplitudes) . We present a numerical comparison with recently available results for the six-gluon 
MHV amplitude at two loops. 

In section [TOl we present our conclusions. 

There are two appendices. Appendix [X] contains an alternative proof of an identity between 
two conformal integrals derived in section 12.61 using the Mellin-Barnes technique. In appendix 
|B] we present the previously unpublished details of the calculation of the four-point Wilson loop 
at two loops. 

2 (Super-) conformal symmetry 

Conformal symmetry in quantum field theory was extensively studied beginning from the late 
1960's, see e.g. |96l 197] and references therein. The conformal properties of correlation functions 
were studied. Conformal symmetry has very strong consequences in two dimensions, where the 
conformal group has an infinite number of generators. In four-dimensional field theories, which 
are more relevant for particle physics, the conformal group has only 15 parameters, and as a 
consequence, it is less restrictive than in two dimensions. Moreover, in generic field theories, 
conformal symmetry is broken by quantum corrections. An example is (massless) QCD. Never- 
theless, studying the deviation from conformal invariance can be useful in practice, as discussed 
in the review [98J. The situation is much better in a field theory which is conformally invari- 
ant also at the quantum level. As we will see in section 12.41 the maximally supersymmetric 
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Yang- Mills theory in four dimensions, N = 4 SYM, which was discovered in |99L 100], has this 
remarkable property. 



2.1 Definition of conformal transformations 

The presentation of conformal symmetry follows roughly chapter four of the book |101j . where 
the reader can find more details. Let us consider a d-dimensional space with flat metric r] pu (the 
treatment for Euclidean and Minkowski space is identical). By definition, conformal transfor- 
mations leave the metric invariant up to a local rescaling 

x^ -» x" 1 , dx^dx^ -» A(x)dx' pdx"* . (8) 

Geometrically, ([8]) means that conformal transformations preserve angles. In order to find the 
most general solution to ©, consider an arbitrary infinitesimal coordinate transformation 

x m _> X 'M = x n + e M( x ) (9) 

It is then easy to show that in d > 2 dimensions, the most general form of compatible with 
(JBJ is given by 

e M = + m^x 1 ' + Ax M + 2(6 • x)x M - 6^x 2 , m Mi/ = -m UfJi , (10) 

where we used the notations b-x = b v x v and X — X X i/ • The two-dimensional case is special and 
we will not treat it here, since we are interested in d = 4. For more information on conformal 
symmetry in two dimensions, see for example [101 j . The infinitesimal transformations in (|10p 
corresponding to the parameters and m^ v are translations and rotations, respectively. Thus 
the Poincare group is a subgroup of the conformal group. The transformations corresponding 
to A and b v are dilatations, and special conformal transformations, respectively. Let us define 
generators for the infinitesimal transformations according to 

x' p = (1 + ioTP^ + im^M^ + iXD + iWKp)x p . (11) 

By comparing (jlip and (|10p it is easy to see that the generators of the conformal group are 
given by 

Pp = -id,j, , 



„2 



D = -ix^d^, 



K p = -i{2x^x u d u -x z d^). (12) 

Here we used the shorthand notation dp = gj^j-. From (|12p we can see that the conformal group 
in four dimensions has 15 generators, i.e. four translation generators P„, six rotations Mav, 
a dilatation D and four special conformal boosts Kp. In addition to the usual commutation 
relations of the Poincare algebra, 

[P M , P v ] = , [A V ,P p ]=i (vupP„ ~ V^pPu ) , 

[Mp,,, M pa \ = i(-r]^M up + r) va M w + r} w M va - r] up M pa ) , (13) 

the generators f 1 1 2 [> have the following commutation relations: 

[D, Pp] = iPp , [K p , Mp U ] = i(rj plt K v - 7] P vKp) , 

[D, Kp] = -iKp , [Kp, P v ] = 2i( Vtlu D - A-V) . (14) 

Relations (|13p and (|14[) define the conformal algebra. 
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Let us now consider finite conformal transformations. For translations P p , rotations M^ v 
and dilatations D, the form of the finite transformations are easy to find. Furthermore, it is 
possible to show that the finite transformation corresponding to the infinitesimal conformal 
boosts is given by 

r m = — — — — — 

l-26-x + 6 2 x 2 ' 1 ; 

For practical purposes it is convenient to introduce another finite transformation, the inversion, 



/ : x» s . 



(16) 



The special conformal transformation (|15p can be obtained by composing an inversion, a trans- 
lation by — and another inversion. Because of this, it will often be convenient to check 
the conformal invariance (or covariance) of a given translation invariant expression by doing 
inversions (|16|) rather than the more complicated special conformal transformations (|15|) . In 
particular, inversions will be useful to us in order to determine the consequences of conformal 
symmetry on correlation functions. Note that the inversion (|16p is an element of the conformal 
group not connected to the identity. In other words, there is no infinitesimal generator corre- 
sponding to (fT6|) . 



So far we discussed the action of the conformal group on coordinates. In a quantum field 
theory, one has to define the action of the conformal generators on fields as well. The conformal 
generators acting on fundamental fields (fii(x) with conformal weight A and Lorentz indices / 
are @ 



M^fa 
Bfa 



{x»3 v - x v 3 p )fa + (m^/fa , 
x • d fa + A fa , 
d» fa, 

(2x p x • d - x 2 d p ) fa + 2x^ Afa + 2x 1/ (m> lv ) I J (j)j 



(17) 



Here m^ v is the generator of spin rotations, e.g., = for a scalar field and {m lll ')\ p = 
g up 5^ — g pp ^\ for a gauge field. From (|17j) one can in principle determine the variation of 
fields under finite transformations. For example, a scalar conformal field has the transformation 
property 



(f>(x) — > 4>{x') 



dx 1 



dx 



-A/d 



where A is the conformal dimension of fax). Here 

dx' 



dx 



A(x) 



-d/2 



(18) 



(19) 



is the Jacobian of the conformal transformation of coordinates, c.f. 
called primary. 



Fields obeying (|T8l) are 



As was already said, we will consider conformal symmetry in four dimensions only. Examples 
of classically conformal field theories are theories which have only dimensionless parameters in 
the action (i.e. dimensionless coupling constant, no masses). An explicit example is the scalar 
fa model in four dimensions, whose action reads 



S 



d*x 



-d p fax)d,fax) + ^fa(x) 



(20) 



4 The generators G determine the infinitesimal transformations with parameters e: 4>'(x) = 4>{x) + £ 
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Under translations and rotations, the scalar field transforms as 4>'(x') = 4>(x), so that (|20p is 
invariant under Poincare transformations. Further, one can see that (|20|) is invariant under 
dilatations and special conformal transformations, if 4>(x) transforms according to (|18p with 
weight A = 1. Note that the coupling constant g in (|2(jp is dimensionless. Dimensional param- 
eters in the Lagrangian, as for example m corresponding to a mass term m 2 (j) 2 (x), would spoil 
conformal symmetry. 

Let us stress that the conformal invariance of (|20p is only valid classically, and it is broken 
by quantum corrections. The reason for this is that due to ultraviolet divergences, the coupling 
constant g and the elementary field <fi have to be renormalised. This is done in the usual way 
by introducing appropriate renormalisation factors Z^^Zg [102J. The net result is that the 
renormalised coupling constant g depends on the renormalisation scale fi, i.e. the /3 function is 
nonvanishing, 

"| = *>*0. (2D 

This is the generic situation in quantum field theory: the conformal symmetry is lost in the 
quantum theory, and may reappear only at certain fixed points g c where (3{g c ) = 0. 

2.2 Conformal correlation functions 

In this section, we briefly review the implications of conformal symmetry for correlation func- 
tions. For a comprehensive review of conformal symmetry in quantum field theory, see e.g. 

MM- 

In a conformal field theory one usually considers primary operators because they have sim- 
ple transformation properties under conformal transformations, see (|18p . Requiring conformal 
covariance of these operators under conformal transformations according to (I18p restrains the 
functional form that their correlation functions can have. As we will see, conformal symmetry 
is most restrictive for two- and three-point functions of scalar primary operators. 



2.2.1 Two-point functions 



Take for example two scalar operators 0\ and 2 possessing conformal dimension Ai and 
A2, respectively. Under conformal transformations their two-point function has to transform 
according to (fl~8j) . i.e. 

' {0 1 (x' 1 )0 2 (x' 2 )) , (22) 



(O l (x 1 )0 2 (x 2 )) 



dx' 




dx' 


dx 


X=X\ 


dx 



X=X2 



where d is the space-time dimension. From invariance under the Poincare group (i.e. under 
translations and rotations) we immediately deduce that their two-point function can only depend 
on the translation and rotation invariant variable x\ 2 = [x^ — x 2 ) 2 -, 

(0 1 (x 1 )0 2 (x 2 )) = f{x\ 2 ). (23) 

Requiring covariance under dilatations with conformal weights Ai and A2 (c.f. equation (|18p ). 
respectively, we find 

(0 1 (x 1 )0 2 (x 2 )) = C[xl 2 ]- {Al+A2)/ \ (24) 

where C is an arbitrary normalisation constant. Finally, covariance under conformal boosts 
requires Ai to be equal to A2 (this is easiest seen by doing an inversion), and therefore we have 

(Oi(xi) 2 (x 2 )) =C [x\ 2 ] ~ A , for Ai = A 2 = A , otherwise. (25) 

Thus, the functional form of (0\(x{) 2 {x 2 )) is completely fixed, and the only dynamically 
determined quantities are the overall normalisation C and scaling dimensions Ai and A2. 
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2.2.2 Three-point functions 

The same reasoning as in the previous section leads to the most general form of the three-point 
function of scalar operators, 

in I \n I \/n I \\ n\ 2 1 -(Ai+A 2 -A 3 )/2 r 2 1 -(A 2 +A 3 -Ai)/2 r 2 i -(A 1 +A 3 -A 2 )/2 
(Oi(xi)0 2 (x 2 )0 3 (x 3 )) = C [x 12 J [x 23 J [x 13 J (26) 

We can immediately check that under an inversion (|16p . (|26|) has the correct conformal weight 
Aj at points i = 1,2, 3. An example of relation (|26p is the well-known star-triangle identity [103J. 

Three-point functions of operators with spin are also constrained. For example, consider a 
spin one operator of conformal dimension Ai = 3 (for example a conserved current). Its 
three-point function with two scalar fields of conformal dimension A 2 = A 3 = 2 is 

(v»( Xl ) o(x 2 ) o(x 3 )) = 2 2 2 y/; 23 , (27) 



x 13 x 12 x 23 
A 1 A 4 

x 12 x 13 



where the vector 



is conformally covariant at point 1 and invariant at points 2 and 3. The generalisation of (J2j 
to tensor fields of arbitrary spin and to arbitrary conformal dimensions is straightforward (see 
for example P-04J). An application to twist-two operators in N = 4 SYM can be found in [105J. 
Three-point functions of spin one operators |106| were studied in the context of anomalies in 

2.2.3 Four-point functions 

An interesting new phenomenon occurs starting from four points. It is then possible to write 
down conformal invariants in the form of cross-ratios, which in general read 

4^ • ( 29 ) 

X ik X jl 



At four points, there are two independent cross-ratios, 

,71 2 aj-.2 2 /~vi 2 

_ X 12 X 34 _ X 14 X 23 



2 2' 22 
X 13 X 24 X 13 X 24 



(30) 



Thus a four-point correlation function contains in general an arbitrary function of u and v, 

(Oi(xi) 2 (x 2 ) 3 (x 3 ) 4 (x 4 )) = — 2^ — 2-/(u, v) , (31) 

x 13 x 24 x 12 x 34 

and the prefactor of / on the r.h.s. of (|3ip carries the overall conformal weight of the external 
points (we have chosen Ai = A 2 = A 3 = A4 = 1 for simplicity). 

Let us remark that correlation functions of operators with spin are also constrained, see e.g. 
[108, 109J. As an example, the most general form of the four-point function of two scalars and 
two spin-one operators is found to be 



(0 1 (x 1 )0 2 (x 2 )Kf(x 3 )^(x 4 )) = 4 



X 12 



+*3V& 2 /3(«, V) + y 3;l 2 n;l 3 /4(^, v) + Y^Yfafsfr, v) 

Here, the conformal tensor 



(32) 



1% = V - 2^ (33) 
x 12 x 12 



carries spin 1 and dimension 1 at points one and two. 
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2.3 Supersymmetric extension of the algebra 

One can extend the Poincare algebra (I13p by supplementing it with fermionic generators satis- 
fying the anti-commutation relations (for an introduction to supersymmetry, see e.g. |11Q|. 1111] ) 

{Qa,Qa} = -2a^P fl . (34) 

The supersymmetry generators Q a and Qd commute with P^, and they transform under the 
representation (1/2,0) and (0,1/2) of the Lorentz group, respectively. It is possible to give a 
representation of the Super-Poincare algebra on a superspace. This space consists, in addition 
to the usual commuting variables x^, of two anticommuting spinors 9 a and 9 a . Then, we can 
define a supersymmetry transformation with anticommuting parameters £ a , £ a by 

5 t x ll = i {t*„6 - 9a , 8<:6 = Z, 5^9 = £ . (35) 

From (I35p we can see that the infinitesimal supersymmetry generators are defined by 

Qa = ^ ~ * ^ ■ Q« = " ^5 + i (^)a ^ ■ ( 36 ) 

They satisfy the anticommutation relation (|34H . 

If one combines supersymmetry generators and conformal generators (note that from com- 
muting them one gets additional generators), one obtains the superconformal group. In section 
12.41 we will study a field theory which has an J\f = 4 extended superconformal symmetry, with 
symmetry algebra PSU(2, 2|4), see e.g. |112| for more details on the algebra. 

We remark that superconformal symmetry constrains correlation functions of (superconfor- 
mal) primary operators in a similar way as discussed in section [2~2l for conformal symmetry. For 
more details, see e.g. [113J . 

2.4 The superconformal field theory J\f = 4 SYM 

Let us introduce the field theory studied in this thesis. N = 4 SYM is the maximally supersym- 
metric Yang-Mills theory in four dimensions. 

2.4.1 Action in components and supersymmetry transformations 

Its action was first found by dimensional reduction of J\f = 1 SYM in ten dimensions [100J. The 
original action in ten dimensions is 

C = Tr (-±F MN F MN + ighr N V N ^ . (37) 

The (trivial) dimensional reduction consists in requiring that the fields in (|37H do not depend 
on six of the ten spacetime dimensions, i.e. 

d*+ m A N = 0, d 4+m $> = 0, d 4+m ^ = 0, m = l,...6. (38) 

Then, one splits the ten-dimensional indices M, N up into four- and six-dimensional ones and 
makes the following definitions. One has to make a choice for the representation of the ten 
dimensional T matrices in terms of four and six dimensional 7 matrices, e.g. (see [1 12J ) 

r M = 7a ,®1, foryu = l,...,4, (39) 
r 4 +m = 75®fm, form=l,...6, (40) 
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where 7^ and 75 are the standard Dirac matrices in four-dimensional Minkowski space, and T r , 
are Dirac matrices in a six-dimensional Euclidean space, which can be written as 



r, 



a 



a m 

- 1 



The scalars are defined as the last six components of the ten dimensional gauge field A^: 



A 



4+m 



for m = 1, ... 6 . 



Finally, with the help of the six dimensional a m matrices we can define 



1 



(41) 



(42) 



(43) 



This leads to the following form of the four-dimensional action 

+ig\i [\j,4> ij ] + ig\ { [X j , fa] + 9 2 \ [<t>ij^kl] 



(44) 



The theory contains a gauge field A^, four complex fermions X at (i = 1,2,3,4), and six real 
scalars (j) 13 = —(ft 1 . All fields are in the adjoint representation of the gauge group SU(N). 

The action following from (I44D can be seen to be invariant under J\f = 4 on-shell supersym- 
metry transformations (which follow from the N = 1 supersymmetry of (|37[) . 



SA 



><t>ij 
SX, 



i&V/iX 1 - iXia^C 
£,iXj — ijXi + €ijki£ X 

--ia^&F^ + ■2ini-P j ,o i ^ + 2ig 



(45) 



where the transformation parameters £ and £ are chiral and antichiral spinors, respectively. The 
algebra of the supersymmetry transformations (|45p closes up to a gauge transformation and 
on-shell (i.e. using the equations of motion). 



2.4.2 Finiteness 

The form and relative factors in (|44p are completely fixed by N = 4 supersymmetry, and for the 
same reason there is just one coupling constant g. A special property of N = 4 SYM is that its 
superconformal symmetry is not broken by quantum corrections. The (3 function of A/ - = 4 SYM 
was shown to vanish up to three loops by direct calculations |114^115|I116| . Furthermore, there 
exist several arguments for the vanishing of the f3 function to all loops |117[ 11181 [119, 120, 1121] . 
For more details and further references, see the review [112] . 

Let us remark on different formalisms in which M = 4 SYM can be studied. Unfortunately, 
an off-shell M = 4 formalism is not available to date, and there are reasons to believe that it 
does not exist. For this reason one has to resort to less supersymmetric formalisms. For exam- 
ple, one can very well use the action (|44p . after the usual gauge fixing procedure, for practical 
calculations. However, it may be advantageous to utilise a manifestly supersymmetric setup, 
either for calculational convenience, or because of conceptual advantages. For example, one 
can take a formulation of N = 4 SYM in terms of M = 1 superfields. We will give the nec- 
essary definitions in the next section, and present a sample calculation in section 12.5.41 It is 
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also possible to use a manifestly j\f = 2 supersymmetric formalism, which employs harmonic 
superspace |122j . For examples of recent calculations in this setup, see [l'2'6\ 1124] , Further, 
there is even a j\f = 3 harmonic superspace formalism [122] . and its quantisation was consid- 
ered in |125j . However up to now, no supergraph calculations were carried out in this formalism. 

Let us clarify a point that might otherwise lead to confusion. j\f = 4 SYM is often referred 
to as a 'finite' field theory. In a supersymmetric formalism, it is indeed true that propagators 
and the coupling constant acquire no or only finite corrections in perturbative calculations. On 
the other hand, in a non-super symmetric gauge (e.g. the Wess-Zumino gauge), the gauge de- 
pendent propagators do get divergent corrections and need to be renormalised by appropriate 
wavefunction renormalisations, and the same is true for the coupling constant. It is only the (5 
function (which is gauge independent) that vanishes. 

Furthermore, even in a finite superspace setup, divergences can and usually do arise when 
one studies composite operators. The latter are traces of products of several operators at the 
same space-time point. Such operators generically have short distance singularities which have 
to be regularised. We will see an example in section 12.51 

2.4.3 Action of J\f = 4 SYM in M = 1 superspace 

The field content of j\f = 4 SYM can be realised in j\f = 1 superspace by introducing three chiral 
superfields $ and one real gauge superfield V . The gauge fixed action in the Feynman gauge, 



Here the dots stand for other vertices involving three and more gluons and also ghosts. We do 
not display them since they will not appear in the one-loop calculation we present later. The 
f a bc are the structure constants of the gauge group SU(N). 

For the study of conformal field theories, it is suitable to write correlators in coordinate 
(super-)space, rather than in momentum space. The coordinate space propagators following 
from (|46p are 




S 




id 



(46) 




(47) 



and 




(48) 



where 




e 



3 ' 




(49) 



The Feynman rules can be read off from (|46p . 
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Except for the definition of the coupling constant, which is 2g here compared to g in 126 
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2.5 Consequences of conformal symmetry for perturbative calculations 

In this section, we give a number of examples to show how the abstract conformal correlation 
functions discussed in section 12,21 are implemented in J\f = 4 SYM. A subtlety arises because 
the gauge fixing procedure usually (e.g. when using a covariant gauge) breaks the conformal 
invariance of the action following from (|44|) . Therefore, the conformal predictions for correlation 
functions of section [2~2l only apply to gauge invariant quantities, in which the gauge dependent 
variation of the gauge fixing term drops out. 

In the next section we give examples of gauge-invariant operators. Then, we introduce the 
concept of anomalous dimensions and briefly discuss operator mixing in a conformal field theory. 
Finally, in section 12.5.41 we show how a particular four-point function can be used to compute 
anomalous dimensions using the operator product expansion. 

2.5.1 Gauge invariant operators 

From what was just said it is clear that in order for conformal properties to be manifest, one 
should consider gauge independent quantities. There are different possibilities. For example, 
one can take a trace over several elementary fields multiplied together at the same space-time 
point. For the gauge group SU(N), the minimal number of fields is two, since Tr(t a ) = 0. 

Let us give two examples which often appear in the literature on J\f = 4 SYM. The first, 

= Tr (W - \5 ij 4>hA , (50) 

is in the representation 20' of SU(4). It is the lowest component of the energy momentum 
supermultiplet. By applying supersymmetry transformations, one can generate the conserved 
currents of N = 4 SYM. For example, in the supersymmetry variation of (|50p the SU (4) current 
of J\f = 4 SYM appears. It reads 

jf = Tr{V»cf> AC <j> CB - <PacV^ cb - \ (\ A ^X B - -f A \ C 0»\ C ) } , (51) 

and satisfies the conservation equation d^J 1 ^ = 0. From the point of view of representation the- 
ory of PSU(2, 2|4), Q20' belongs to a short representation of fixed conformal dimension A = 2 
|127j . Thus, its conformal dimension stays at its classical value to all orders in perturbation 
theory. Furthermore, it is known (see e.g. the review |128j ) that two- and three-point function 
of (foT)|) do not receive quantum corrections. 

The second example, 

K. = ^Tr , (52) 

is a singlet of SU(4:). It is the lowest component of the Konishi supermultiplet |129|, 130j . 

Finally, one can also define non-local gauge invariant operators. An example is the Wilson 
loop 

W[C] = ^TrPexp (ig f dx^A , (53) 

which is a functional of the (closed) contour C. Wilson loops defined on particular polygonal 
contours are in fact one of the main objects of study of this thesis. We will discuss them in 
much more detail in section [SJ 
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2.5.2 Anomalous dimensions 

In a quantum field theory, the conformal dimension A of an operator depends in general on the 
coupling constant, and one writes 

A(g) = A + 7 ( 5 ) , (54) 

where Ao is the classical, i.e. tree-level dimension and 7(5) is the anomalous dimension of the 
operator. 

Let us start with the operators Q 2 "o'' ^ or which = Ao- For practical calculations, it 

is convenient to employ the M = 1 superfield formalism presented in section 12.4.31 Since the 
SU(4) symmetry of M = 4 SYM is not manifest in this formalism, we have to decompose 
SU(4) — > SU(3) x U(l). Under this decomposition, Q 1 ^, gives rise to the following operators 
(among others): 

C U = 'p r , c ] u = Tr . (55) 

Here I, J = 1, 2, 3. It is convenient to use them for quantum calculations since they are gauge 
invariant without the need to include gauge links like exp(2^y). Their lowest component is 



C IJ = Tr (0V) , Cjj = Tr {fify . (56) 

Let us compute the two-point function of C 11 and C\± in order to check the absence of pertur- 
bative corrections at one loop. We find that the sum of the graphs contributing to the one-loop 






Figure 1: Feynman graphs contributing to {C 11 {x\)C\ 1 {x2)) and {K{x\)K{x2}) at one loop. Solid lines 
denote the chiral superfield propagator (|T7|) . wiggly lines the gluon superfield propagator 



propagator corrections, Fig. H](b,c), vanishes identicalljo The graph in Fig. Q] (a) leads to an 
integral that can be shown to vanish for x\2 7^ 0. Thus, up to contact terms, we have 

(C u ( Xl ) Cl(x 2 )) = + O(g') , c = . (57) 

As we already said, ([5?]) is in fact valid to all orders in the coupling constant. Comparing ((57)) 
to (|25p . one can see that C 11 has the conformal dimension A = 2, which is just the sum of the 
classical dimensions of its constituents (f) 1 . 

The second example is the Konishi operator (|52)) . In the M = 1 formalism it is given by 

K = |tr (e-W^eW® 1 ) . (58) 

The exponentials in (|58|) are needed to make tC gauge invariant. It has the same classical 
dimension Aq = 2 as Q20" unn ke Q20" ^ does get renormalised. Starting from one loop0, 



6 This is true in the Feynman gauge. 

7 By a calculation to n loops, we mean up to order g 2n in the coupling constant. For composite operators, the 
pictures would suggest a different loop order. 
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Figure 2: Additional Feynman graphs contributing to (K.(xi)K.(x2)) at one loop. Solid lines denote the 
chiral superfield propagator (|T7|) . wiggly lines the gluon superfield propagator 



/C has UV divergences typical for composite operators. They come from diagrams (a) and (b) 
in Figure EJ Evaluating them in dimensional regularisation, one finds (for simplicity, we will 
compute the lowest, i.e. 9 = 9 = 0, component of (K,{x\) K.{x2)) only) 

FigUa) + FigH&) = , 2 \ 2 _ 2e g 2 {x\ 2 n 2 ) 

The first factor on the r.h.s. of (|59[) accounts for the engineering dimension of (ZC(xi) K,{x 2 )) 
in dimensional regularisation, and in addition one gets a factor of g 2 (x 2 2 fi) e for each loop, with 
fx being the dimensional regularisation scale. Furthermore, diagrams of the same topology as 
those in Fig. Q] contribute, but the only difference is the orientation of some matter lines, and 
they still vanish for x\% ^ 0. The graph in Fig. 0c) leads to a finite contribution. 



1{N 

€ 



1)N 



(2tt) £ 



+ 0(6°) 



(59) 



In dimensional regularisation, one defines renormalised operators by multiplying by a suitable 
renormalisation factor Z [102J. This means that the renormalised operators [K] are defined by 

[K]=ZK, Z = l + ^V + 0( 5 4 ). (60) 

The Z-factor subtracts the poles in e, so that one is left with a finite result. From ([59]) we 
can see that one should take Z\-\ = (3N)/(8tt 2 ). Expanding in e, we find for the renormalised 
operator [K], to one loop 



<[/q(x 1 )[/q(x 2 )> = ^-^4- 

l Z7 V X 12 



N 



+ 0( 5 4 )+0(e), (61) 



with c(g 2 ) = l/3 + g 2 N/(8Tr 2 ). One might think that the appearance of logarithms is surprising 
in a conformal field theory. The reason is that to a given order in the coupling constant, the 
expansion of (|25p leads to these logarithms. Indeed, one can rewrite (|6ip in the manifestly 
conformal form (1251). 



N 2 - 1 

([/C](x0[/C](x 2 )) = -j-^- [c(g 2 ) (x 2 2 y^ (^)--j +0 (g*)+0(e), (62) 
with the conformal dimension of K, and anomalous dimension 7^ at one loop being 

A K (g 2 ) = 2 + lK {g 2 ) : 7c = ^ 6 + 0( 5 4 ). (63) 
2.5.3 Operator mixing 

Let us briefly remark on the topic of operator mixing (cf. e.g. [102] ) in the context of a conformal 
field theory [1044 11311 1132} 11331 1134] . Operator mixing occurs, generally speaking, when there 
are several operators that have the same quantum numbers. A particular example that often 
appears in the literature are operators in the s/(2)-sector of N = 4 SYM, see e.g. [17} I124j . Of 
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these, the twist-two operators^] are built from two elementary fields, e.g. (j) 1 and an arbitrary 
number j of (covariant) derivatives T>^. For j = 0, we already encountered one of these operators 
in the form of C 11 = (c/> 1 (/> 1 ). For general (even) j, they read 

j 

3 = CjkO jk , O jk = tr (V^ ...TP* 1 . . . P«> 1 ) . (64) 

fc=0 

Here, the curly brackets stand for traceless symmetrisation of the vector indices, so that Oj is in 
a representation of spin j. In a generic situation, there is more than one possible distribution of 
the derivatives on the two fields 4> l , which is reflected by the sum in (|64|) with a priori arbitrary 
coefficients Cjk- As a consequence, the operators Ojk for a given j mix under renormalisation. 

In a generic quantum field theory, the mixing matrix Cjk at order g n can be determined by 
a calculation at order g( n+2 h In a conformal field theory, this mixing problem can be resolved 
using conformal methods that significantly reduce the complexity of the calculations |104} 1131] . 
We developed these methods further in [135]. As a particular example, we were able to deter- 
mine the mixing matrix Cjk at order g 2 by evaluating only order g Feynman graphs, compared to 
a conventional calculation at order g . Similarly, the computation of the anomalous dimensions 
of the twist two operators in (|64|) can also be reduced in loop order (by one unit of g 2 ). 

One might object that the resolution of the mixing problem is not relevant, since it is known 
that mixing matrices are scheme dependent. This scheme dependence, however, is under control 
and is governed by a renormalisation group equation. Moreover, one can define a conformal 
scheme, in which the properties of conformal correlation functions are realised (see e.g. the 
review [98]). In contrast, in a generic scheme, such as for example minimal subtraction, the 
conformal properties discussed in section [2T2l are in general not present. 

The anomalous dimensions of the operators in (|64l) have received considerable attention 
recently. They are known up to two loops and conjectured at three loops [136] . They are 
also predicted by a conjectured integrable model [30J. It would be desirable to confirm these 
conjectures to three loops and beyond. A promising technique to achieve this is the determina- 
tion of the anomalous dimensions through the operator product expansion (OPE) of conformal 
four-point functions. We conclude this introductory section with a one-loop example of this 
technique. 

2.5.4 Conformal four-point functions 

Let us consider the four-point function of Q^q, . We already saw in section 12.5.21 that Q l ^ , does 
not need to be renormalised. Its conformal dimension is equal to its classical dimension, A = 2. 
In this case, from formula (I3ip we find for the four-point correlator 

(c 11 {x l )C 22 (x 2 )C\ 1 {x z )C^ 22 (x i )) = aa^a f(u,v;g 2 ) (65) 

The correlator (|65p is known up to two loops in perturbation theory \1'67\ [J26J. One reason 
for being interested in it is that the function f(u,v;g 2 ) allows to determine the anomalous 
dimensions of the twist-two operators (|64p via the operator product expansion (OPE). This 
relation has been worked out in [138] 1139] . 

As was already discussed, the two-point functions of chiral primary operators that could in 
principle contribute to the disconnected piece of (I65p vanish (up to possible contact terms, see 
e.g. [140J for a discussion). The only nonvanishing Feynman graph contributing to (|65p at one 

8 The twist of an operator is denned as the difference between its classical dimension and its spin. 
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Figure 3: The only connected Feynman graph contributing to (C 11 (xi) C 22 (x2) C T n (#3) 22(^4)) at 
one loop. Solid lines denote the chiral superfield propagator (jTf)) . The dots represent the (chiral) 
superspace integrations at points x§ and x§. 

loop is shown in Fig. O For (|65|) we need its lowest component only, i.e. we can set all external 
#'s and 9's to zero. In that case, using the superfield propagator (14"T|) . we obtain 



Fig.ES = 



g 2 N(N 2 - 1) 1 



d x§d xq I d 9§d 9§ exp {iO^d^Oo) 



,.— ^yi 2 ,-y>2 ^yi 2 ™2 

x 15 x 25 x 56 x 63 x 64 



• (66) 



The notation d§§ means that the derivative acts on the term l/xg 6 only. The integrals over 9 
and 9q 'see' only the quadratic term in the exponential, i.e. 

d 2 9 5 d 2 9 6 exp (i9 5 d 56 9 6 



where □ = d^d^. Then, we use that 



□ - 



x- 



-4ttV 4 )(x) 



to undo one of the space-time integrals in (j66H . We arrive at 



g 2 N(N 2 - 1) 1 
2(2*)" *?2 X §4 



, X 2 ,X 3 ,X4) , 



with 



X2,X 3 ,X 4 ) 



d A x 5 



2 2 2 2 

zy»ij /-y*4J ry**- 1 ry**' 

x 15 x 25 x 35 x 45 



2 ^2 



-d>«< 



X-\nX 



(67) 



(68) 



(69) 



(70) 
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Here x™ = Xj — Xj and the conformal cross-ratios u and v were defined in (|30|) . The fact that 



Figure 4: The one-loop ladder integral. Each line represents a propagator with the integration point given 
by a solid vertex. The reason for the names ladder and box is clearer in the momentum representation 
of the same integral. 

the integral is characterised by a single function of two variables follows from its conformal 
covariance [141] . Indeed, performing a conformal inversion on all points, 



x' J 



x- 



x~ 



10 



1 2 ' 

xfxj 



d x 5 



d x 5 



.8 ' 



(71) 
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we find that the integral transforms covariantly with weight one at each point, 



h^'(xi,X2, X 3 ,X4) > x\x\x\x^ X ' {xi,X2,X3,X4). (72) 

Since rotation and translation invariance are manifest, we conclude that the integral is given by 
a conformally covariant combination of propagators multiplied by a function of the conformally 
invariant cross-ratios ([30]) . in agreement with (j65l) . 



The function §>W(u, v) has been calculated in [1421 1143] , where it was also shown that the 
same function appears in a three-point integral. The latter can be obtained from the four-point 
one by sending one of the points to infinity [141] . We can multiply equation (fTOl) by xf 3 , say, 
and then take the limit X3 — > 00. This gives, 

/d^x 1 
-9 2~^~2~ = -2~ $(1) (^^)' ( 73 ) 
X lA X lb X 24 

where the cross-ratios u and v have become u and v in the limit, 



u ^u=^, v ^y = rtk. (74) 

L 24 x 24 



.1 



Thus the three-point integral contains the same information as the four-point integral, i.e. the 
same function of two variables. The reason is that one can use translations and conformal 
inversion to take the point X3 to infinity and the function of the cross-ratios is invariant under 
these transformations. The explicit formula for $W is [144] ^ 



A 



2(Li 2 (—pu) + Li 2 (—pv) ) + In - In h ]n(pu) ln(pv) + — 

u 1 + pu d 



(75) 



where 



A(u, v) = -u-v) 2 - Auv , p(u,v) =2(1 -u-v + X)' 1 . (76) 

In order to extract the anomalous dimensions of twist two operators, we need ^^(u, v) in the 
OPE limit X12 — ► 0, i.e. u — > 0, only. In this limit, one finds (taking < u < 1 for simplicity) 

$(!) (« «) = [In (u) ln(ti) + 2Li 2 (1 - v)} +0(u), (77) 

1 — v 

Using the results of [139] , one can deduce from ([7?]) the well-known expression for the one-loop 
anomalous dimensions of the twist- two operators (I64p . 

7(j)=4a£^ + 0(« 2 ), a=0. (78) 
fc=l 

One can in principle determine the anomalous dimensions at higher orders in perturbation theory 
from the four-point function (|65|) . It has been computed to two loops in perturbation theory 
[1371 1126] . It would be very interesting to determine its three- loop, i.e. order g 6 , correction 
in order to test the conjectured formula [136] for the twist-two anomalous dimensions at three 
loops. Moreover, one may wonder whether the conjectured integrability of J\f = 4 SYM contrains 
the four-point function itself, beyond the contraints coming from conformal symmetry? 



'Valid for it, v > 0, see [145] for a discussion of the analytic continuation. 
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2.6 Conformal four-point integrals 



In section 12.5.41 we saw that from the four-point function of chiral primary operators (|65p one 
can deduce the anomalous dimensions of the twist- two operators (|64p . Due to the absence of 
divergences of the chiral primaries, such four-point functions are finite and do not require any 
regulator. Because of its conformality, a calculation at a given loop order (e.g. with Feynman 
diagrams) must eventually yield an expression in terms of (finite) conformal integrals. In |77j . 
we took this as a motivation to study an infinite class of conformal four-point integrals. We 
found that there are simple identities between integrals corresponding to different diagrams. We 
describe these 'magic identities' in this section. 

At the same time, let us stress that in this report, the interest in the conformal integrals 
mainly comes from an entirely different motivation. We found that, very surprisingly, conformal 
integrals also appear in the completely unrelated context of on-shell gluon scattering amplitudes 
in momentum space. As will be explained in section 13.21 their appearance suggests that the 
latter have a (broken) conformal symmetry in a dual space defined through the gluon momenta. 
We hope that the example of (|65p . which illustrates where conformal integrals normally appear 
(off-shell, finite, in x-space) helps clarifying the difference to the integrals that appear when 
analysing on-shell gluon scattering amplitudes (which are on-shell, IR-divergent, and defined in 
momentum space) in section [3~2l 

2.6.1 Proof of conformality 

We will discuss an infinite class of conformal four-point integrals in four dimensions^], each of 
which is essentially described by a function of two variables. We already discussed the simplest 
example, the one- loop ladder integral (I70p . in section T2.5.41 It is the first in an infinite series of 
conformal integrals, the n-loop ladder (or scalar box) integrals, which have all been evaluated 
[144] . In particular the 2-loop ladder integral is given by 

h {2) (x 1 ,x 2 ,x 3 ,x 4 ) = x\ A I l X T f ? o 2 = zJ-JT® {2) ^ v ) ■ ( 79 ) 

J x 15 x 25 x 45 x 56 x 26 x 46 x 36 x 13 x 24 

The prefactor x| 4 is present to give conformal weight one at each external point. 



2 




4 



Figure 5: The two-loop ladder integral. The dashed line represents the numerator x\±. 

Again conformal transformations can be used to justify the appearance of the 2-variable 
function <I>( 2 ). The r.h.s. of (|79p is invariant under the pairwise swap x\ < — > x 3 < — > X4, 
hence 

h ( - 2 \x 2 ,x 1 ,x 4: ,x 3 ) = /i (2) (xi,x 2 ,x 3 ,x 4 ) . (80) 

This symmetry is not immediately evident from the integral. It is its conformal nature which 
allows this identification. At three loops we consider two conformal integrals, the three-loop 

10 In this section we consider and prove identities for Euclidean integrals. The corresponding Minkowskian 
version of the identities can be obtained through Wick rotation of the integrals. In the Euclidean context we 
consider integrals with separated external points, Xij 7^ 0. This is the Euclidean analogue of the off-shell regime, 
Xij ^ 0, for a Minkowskian integral. 
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4 4 



Figure 6: The two-loop turning identity obtained from the pairwise point swap, x% < — ► X2, X3 < — ► 2:4. 
ladder, 

/on 4 /" dx* > dx§dx r r 1 ,fsl/ 

/l v - , (Xi,X 2 ,X 3 ,X 4 ) = X 24 / -5 o 2 2 2 2 2 2 2 T = ~2 2~ $ ( M > u ) > ( 81 ) 

J x 15 x 25 x 45 x 56 x 26 x 46 x 67 x 27 x 47 x 37 x 13 x 24 

and the so-called 'tennis court' [39] . 

(»\, n f d 4 X5(i 4 X6(i 4 X7 1 T (TW 

Q K >(xi,X2,Xz,Xi) = X 24 / -g 5 2 2 2 2 2 2 2 2~ = ~2 2~* K V ) > ( 82 ) 

J x 15 x 25 x 45 x 56 x 57 x 67 x 26 x 47 x 36 x 37 x 13 x 24 

which are shown in Fig. [7J Notice the presence of the numerator x§ 5 in the integrand of the 



2 2 




4 4 



Figure 7: Two examples of three-loop conformal four-point integrals, the three-loop ladder and the 
'tennis-court'. 

tennis court. It is needed to balance the conformal weight of the five propagators coming out of 
point 5. 

2.6.2 'Magic identities' between conformal integrals 

We will show that the three-loop ladder and the tennis court are in fact the same, i.e. we will 
prove <3?( 3 ) = \E'( 3 ). First we shall present a diagrammatic argument. We consider the n-loop 
ladder as being iteratively constructed from the (n — l)-loop ladder by integrating against a 
'slingshot' (the '0-loop' ladder is a product of free propagators). For example we write the 
three-loop ladder as 

/l (3) (xi,X 2 ,X 3 ,X 4 ) = X 24 / 2 2 I2 (^24 / 2 x 2 „2 I2 J „2 „2 ) ' ( 83 ) 

J x 15 x 25 x 45 J x 56 x 26 x 46 x 67 x 27 x 47 x 37 

where inside the parentheses we recognise the two- loop ladder integral ()79|) . We can then show 
the equality of the three-loop ladder and the tennis court by using the turning symmetry (|80|) 
on the two-loop ladder sub-integral. Then the tennis court integral (|82p can be recognised as 
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4 



4 



Figure 8: The three-loop ladder expressed as the integral of the two-loop ladder against the 'slingshot'. 
The empty vertex is the point x§ which must be identified with the point x$ from the two-loop ladder 
sub-integral before being integrated over. 



the turned two-loop ladder integrated against the slingshot, 



/l (3) (xi,X 2 ,X 3 ,X 4 ) = X24 



d 4 x 5 



r 2 

x 24 



2 2 2 
x 15 x 25 x 45 

d 4 x 5 

111 

zy»^ fy**-^ rv*^- 1 

x 15 x 25 x 45 



,X 2 ,X 3 ,X4) , 
,X 5 ,X4,X 3 ) , 



= 5- (3) (xi,X2,X3,x 4 ) . (84) 
This proof can be more easily seen in the diagram (Fig. [9]). In using the turning identity (|80f) 








Figure 9: Diagrammatic representation of the proof of equality of the tennis court and the three-loop 
ladder. The identity follows from the turning identity (|80j) for the two-loop subintegral. 



we have ignored the possibility of contact terms. These could, in principle, spoil the derivation 
of identities like $( 3 ' = as the proof (I84p involves turning a subintegral. Contact terms 
could then generate regular terms upon doing one further integration. We now give an argument 
why this cannot happen for any conformal four-point integral. We again use the example of the 
3-loop ladder and tennis court identity. 



Consider inserting the n-loop subintegral (the 2-loop ladder in this case) into an H-shaped 
frame with a dashed line across the top, as illustrated below. This generates an (n + 2)-loop 
integral which is conformal with weight 1 at each external point (provided the subintegral is 
conformal with weight 1 at each external point). When inserting the 2-loop ladder in this way 
the 4-loop integral one obtains is 



f ( - 4 \x 1 ,x 2 ,x 3 ,x i ) 



T 2 
x 34 



d 4 x^d 4 XQ 



ry 2 ™2 rf, 2 ™ 2 ™2 

x 15 x 45 x 56 x 26 x 36 



r 2 
x 35 



d 4 Xjd 4 X8 



2 2 3 2 2 2 2 
x 67 x 57 x 37 x 78 x 58 x 38 x 48 



x 13 x 24 



f(u,v) 



(85) 
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Figure 10: The 2-loop ladder inserted into an H-shaped frame, generating a 4-loop integral. 

As usual, the second equality follows from conformality. Now we consider the action of Di on 
the above integral using 

□J- = -47rV 4) (:c). (86) 
x z 



On the integral one obtains 



lT 2 ^§4^13 f d*xed*xjd*X8 _ 4tt 2 x| 4 (3) 

2 / 2222222222 422 v") • V u ' / 

ry" I ry " ry" /T"* ry" ry" ry" ry-" ry " ry" ry " ry^ ry " ryi ^ 

x 14 J x 26 x 16 x 36 x 67 x 17 x 37 x 78 x 18 x 38 x 48 x 13 x 14 x 24 



On the functional form of (|85p one uses the chain rule to derive the action of a differential 
operator on the function /. In this way we find the differential equation, 



ry 6 ry ■ ^ ~* ~y — ■ 't ' 

X 13 X 24 x 13 x 14 x 24 

(2) 

The operator AuJ is given explicitly by 

A$ = ^ + ^ + (« + u - + 2<9 U + 28 v . 
Similarly we can act with O2 on the 4-loop integral to obtain the following integral, 

_ 47r 2^i f d A x 5 d 4 x 7 d A x 8 x 2 35 4tt 2 x 2 u (3) 

2/ 2222222222 224 V ' / ' 

ry" I ry-" ry" ry " ry" ry" ry" ry" ry" ry" ry" ry" ry" ry^ 

x 23 J x 15 x 25 x 45 x 57 x 58 x 78 x 27 x 48 x 37 x 38 x 23 x 13 x 24 
and the corresponding differential equation, 



2^2 



X ^A^f(u,v) = -^4^^\u,v). (91) 



fi 4 "itu J 2 2 4 

,'v» \j 'y * ~ ry" ry" • y * ■ 

x 24 x 13 x 23 x 13 x 24 



From (|88I91|) it follows that $( 3 ) = ^( 3 ), the point being that one obtains the same differential 

(2) 

operator A st under the two □ operations. The argument has the obvious generalisation of plac- 
ing any conformal integral (in any orientation) inside the frame. This argument indirectly shows 
that the previous argument (I84p based on turning the subintegral cannot suffer from contact 
term contributions. 
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The identity we have obtained at three loops is just the first example of an infinite set of 
identities which all come from the turning symmetry of subintegrals. We generate (n + l)-loop 
integrals by integrating n-loop integrals against the slingshot in all possible orientations. The 
resulting integrals are equal by turning identities of the form (f5Uj) . At two loops we get just 
one integral (the two-loop ladder). At three loops we have already seen two equivalent integrals 
(ladder and tennis court). At four loops we generate two equivalent integrals from the three-loop 
ladder and three equivalent integrals from the tennis court. Finally, all five four-loop integrals 
obtained in this way are equivalent by the three-loop identity for the ladder and tennis court 
(see Fig. ITT]) . 

In general it is more common to give the diagrams in the 'momentum' representation (which 
has nothing to do with the Fourier transform) where we regard the integrations as integrals 
over loop momenta rather than coordinate space vertices. This representation is neater but 
the numerators need to be described separately as they do not appear in the diagrams. The 
momentum-space version of the four generations of integrals from Fig. [TTJ is given in Fig. [T7[ 
The transition between the two notations will be discussed in more detail in section [3] As was 
already mentioned, integrals related to the ones discussed here will reappear in the analysis of 
loop corrections to gluon scattering amplitudes in section [321 




Figure 11: The integrals in a given row are all equivalent. They generate the integrals in the next row 
by being integrated in all possible orientations against the slingshot attached from above. The ladder 
series is in the left-most column. 



3 Gluon amplitudes in M = 4 SYM 

3.1 Introduction 

What are the motivations for computing scattering amplitudes in A^ = 4 SYM? There are several 
answers to this question. Firstly, being four-dimensional gauge theories, N = 4 SYM and QCD 
share some properties, but calculations are generally much harder in the latter. Having efficient 
calculational methods for computing QCD amplitudes is important for collider physics. The 
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Figure 12: The momentum notation for our integrals up to four loops. The slingshot translates into 
the top box in each diagram, beneath which are the integrals at one loop lower, arranged in all possible 
orientations. The ladder series is again in the left-most column. 

reason is that in searches for new physics, the theoretical uncertainties for the standard model 
background processes are often high, in particular in proton-proton collisions [146J . N = 4 
SYM can be used as a testing ground for new calculational methods. Moreover, a part of QCD 
amplitudes can be predicted from their N = 4 counterparts. For example, tree-level gluon 
amplitudes in pure Yang-Mills are identical to those in J\f = SYM. This follows simply from the 
Feynman rules. Further, at one loop, gluon amplitudes in pure Yang-Mills can be decomposed 
in the following way: 

A g = ( A g + 4Af + 3A s ) -A(A f + A.) + A s . (92) 

JV=4 JV=1 

Here the abbreviations g, f, s stand for gluon, fermion, and scalar, respectively. They denote 
the particles circulating in the one-loop diagrams. In this way, a generic one-loop Yang-Mills 
amplitude can be written in terms of two supersymmetric amplitudes and an amplitude with a 
scalar circulating in the loop, which is much simpler than the original amplitude. 

Another motivation comes from the fact that, as we will discuss in this section, scattering 
amplitudes contain the cusp anomalous dimension r cusp , which has received considerable atten- 
tion over the last years in the study of the AdS/CFT correspondence. Its value is predicted (in 
principle at any given order) from conjectured integrable models that describe the spectrum of 
anomalous dimensions in N = 4 SYM. Therefore, knowing r cusp to high orders in perturbation 
theory is important to test and fine-tune these models. The three- and four-loop values of r cusp 
were indeed computed via four-gluon scattering amplitudes. 

Gluon scattering amplitudes in M = 4 are also studied within the relation between M = 4 
SYM and N = 8 supergravity |147j . see |148j for a review. This may give new insights into 
the ultraviolet properties of J\f = 8 supergravity, see [149] and references therein for a recent 
discussion. 
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Pi Pr, 

Figure 13: n-gluon scattering 



Finally, the scattering amplitudes themselves reveal interesting properties, on which we will 
focus in this report. As will be described in this section, an iterative structure for the scattering 
amplitudes in N = 4 SYM was uncovered by Anastasiou, Bern, Dixon and Kosower (ABDK) 
[56j and generalised to higher loops by Bern, Dixon and Smirnov (BDS) [39]. In particular, it 
turns out that their finite part seems to be much simpler than could be expected on general 
grounds. These discoveries give hope that scattering amplitudes in the large N limit of N = 4 
SYM may be solvable. 

Let us consider an n-gluon scattering amplitude in N = 4 SYM, as shown in Fig. [TBI All 
gluons are treated as ingoing, furthermore they are massless and on-shell. They are charac- 
terised by their momenta pi (with pf = 0), their helicity h{ = ±1 and their colour index (in 
the adjoint representation of the gauge group SU(N)). There is momentum conservation, i.e. 
Yl?=iPi = 0- Taking the gluon momenta on-shell introduces infrared (IR) divergences that 
have to be regularised. For this one usually uses a supersymmetry preserving regulator, such 
as for example dimensional reduction (or a variant of it, the four-dimensional helicity scheme 
|15m ri51] ) . We will discuss the structure of IR divergences in section [3.1.21 

A generic contribution to an n-gluon scattering amplitude can consist of a vast number of 
terms. For this reason it is very convenient to classify the amplitudes according to their quantum 
numbers. Let us start by considering the colour structure of the amplitudes. It is clear that 
the gauge theory factor of a given contribution to the amplitude can be written as a trace or 
a product of several traces over colour matrices. We are interested in the large N limit, so we 
neglect multiple traces. It is convenient to define partial amplitudes which correspond to one 
particular colour structure, 

An = Tr[T a ^T a * . . .T a "]A^-^( Pl , P2 , . . . ,p n ) + . . . . (93) 

The (planar) partial amplitudes An 1 ' 2 '"'' n (pi, P2, ■ ■ ■ ,Pn) depend on the helicity configuration 
and on the momenta of the gluons only. 

3.1.1 Helicity structure 

The next simplification comes from classifying the different helicity configurations. In order to 
do this, let us first introduce the necessary spinor helicity notation. The massless gluon momenta 
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p^ = <y^ ld p aa can be written as a product of commuting spinors, 

p™ = XfX$. (94) 

We also introduce a shorthand notation for spinor products, 

(pq) = \p\%e a p = \p\ qa , [pq] = X$X%e. $ = X p ^ . (95) 

This allows to write for example a scalar product between two light-like vectors p^ and q^ as 

(p" + q»f = 2p-q = 2{pq)[pq]. (96) 

The gluons in the partial amplitudes An 1,h2, '" ,h ™ (pi,P2, ■ ■ ■ ,Pn) are projected according to their 
helicities hi = ±1 with the following polarisation vectors: 

= iki P \ ' e M (?) = t^)" J? ^ i • ( 97 ) 

V2(gp) V2[gp] 

Here g M is an arbitrary (up to g' 1 ^ p^) reference momentum. A change in q^ corresponds 
to a gauge transformation. More information on the spinor helicity notation can be found for 
example in [1521 1153] . 

In = 4 SYM, one can define analogous scattering amplitudes where some of the gluons are 
replaced by fermions or scalars. By requiring invariance of the S'-matrix under supersymmetry 
transformations and the invariance of the vacuum state 1 0} , |154} 1155] , one can derive relations 
between scattering amplitudes involving different particles and/or helicities. The action of the 
M = 4 supersymmetry generators on gluons g^, fermions f i , and scalars Sij is [152J 

[Qi(q,e),g + (k)) = 9[kq]ff, 

[Qi(q, 0), ff{k)] = 65^ (kq) g + + 9[kq] Sij , 

[Qi(q, 0), s jk (k)} = OSij (kq) f+ - 95 ik (kq) f+ + 9[qk]e ijkl f~ , 

[Qi(q,9)Ji(k)) = eS i:j [qk)g- + 1/2 e(qk)e ijkl s kl , 

[Qi(q,6),g-(k)) = 8(qk)fr. 

Supersymmetry can be used to derive relations between some amplitudes. For example, acting 
with Qi(q,0) on /0|/r(l) 5 +(2)...< ? +(n)|0V one finds 



= 5i d (ql) ((%+(l) g+ (2) . . . g + (n)\0) + [ql] <0|a o -(l) g + (2) . . . g + (n)\0) . (98) 

Here we have used that amplitudes with two f + vanish because of fermion helicity conservation. 
From ([HED it follows that both terms in it vanish. We are mostly interested in gluon scattering, 
so we note 

((%+(!) g + {2) . . . g+(n)\0) = A++-+ = . (99) 



In a similar manner, it can be shown by acting with Qi(q, 6) on \Q\g (1) 7^(2) 5 + (3) • • • g + (n)\0 
that 

A -+... + = Q ^ 100 ^ 

Furthermore, relations between amplitudes with different particles can be obtained from su- 
persymmetry. Given (|99p and (jlOOp . the first non-zero pure gluon amplitude has two flipped 
helicities, and is called maximally helicity- violating (MHV), 

A (MHV) = A — +...+ ; A -+-+...+ ^ _ _ (1Q1) 



For n > 3, because for n = 3 the prefactor (ql) vanishes for any choice of q. 
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We remark that for four and five gluons, all gluon amplitudes are MHV (or MHV, which have 
all helicities plus except for two, e.g. A§ h++ ), e.g. 

{A++-, Ai~ + -, ...}, {A+++-, A+-+-, ...}. (102) 

Amplitudes with different helicity ordering are related by cyclicity. It can be shown that in 
J\f = 4 SYM, at four and five gluons there is only one independent MHV amplitude. 

Amplitudes with more flipped helicities are called NMHV (next to maximally helicity- 
violating), NNMHV, and so on. NMHV amplitudes (where three helicities are flipped) start 
appearing at six points, e.g. 

At ++ — , ^ 6 + - ++ ,.... (103) 

MHV amplitudes have the special property that for them, there is only one possible helicity 
structure, and hence it can be factorised. Using the shorthand notation {piPj) = (ij), we have 
for n > up to an inessential normalisation 

4 tree) (l + , • • • , I" , (l + 1) + , • • • , m" (m + 1) + , . . . , n+) = i ^^p , (104) 

and all loop corrections are proportional to A^ tTee \ It is then natural to define the loop correc- 
tions to MHV amplitudes by factorising out the tree amplitude, i.e. 



A n = A^ M n = At ee) \l + M« +Jl4 2) + ... 



(105) 



Here M n is a function of the momentum invariants pi ■ pj and of the dimensional regulator only. 
We already said that four- and five-gluon amplitudes are special because all of them are MHV. In 
contrast, amplitudes with more than two flipped helicities have a more complicated form. Their 
loop corrections involve in general new helicity structures, and hence they are not proportional 
to the tree level term. It remains an interesting open question how to extend the duality between 
Wilson loops and MHV gluon amplitudes discussed in this report to those amplitudes. 

3.1.2 Infrared divergences 

As we already mentioned, the on-shell amplitudes suffer from infrared divergences. In dimen- 
sional regularisation, the divergences appear as poles in the regulator e. The structure of these 
divergences is well understood in any gauge theory [58l [59], EOj EH E21 [631 Ell E51 ESI EH EE1 ESI 
[701 HH EH]- In order to explain the general structure, let us proceed by giving an example. 
Consider the one-loop corrections to the four-gluon amplitude. They can be written as [155] 
(with a = g 2 N/(8TT 2 )) 

\-^stlW( s ,t;e) + 0(a 2 )\ , (106) 



A A = A\ iee 



where the one-loop scalar box integral 1^ is defined by 

2 3 

/ 



j(D 



° I k 2 {k-p l ) 2 (k-p l -p 2 ) 2 (k+p A ) 2, (10?) 



where C = fi 2e e e7E (47r) 2+6 and s = (pi + pi) 2 and t = (p 2 + P3) 2 are the usual Mandelstam 
variables. 1^ has infrared divergences coming from two regimes: the soft regime where k 2 ~ 0, 



2 One can also define on-shell three-particle scattering for complex momenta. 
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and the collinear regime where ~ p^, with p^ being one of the on-shell momenta entering the 
scattering process. 



Evaluating ([1U7P one finds, up to terms vanishing as 0(e), 

1 



.4, 



A 



tree 



e 2 V -a 



a /VV 
e 2 V -* 



1 + a [ - In* - + 4C 2 



+ 0(a 2 ,e) 



(108) 



We can see from (|108|) that the loop corrections to A^ factorise into divergent pieces, each 
depending on one kinematical variable only, and a finite part. This factorised structure is 
known to hold to higher loop orders as well (in the planar case). For non-MHV amplitudes, 
the only difference compared to (|108|) is that the finite part in general depends on the helicity 
structure of the amplitude as well. As we already pointed out, at four points there are only 
MHV amplitudes, and this problem does not appear. Their general form is known to factorise 
, ED dSU [T37] according to 



M 4 = 4 4 M t 4 ree = div f^- \ div (^\ fin (i) + 0(e) 



(109) 



where div(x) contains all poles in e and fin(rr) is finite as e — > 0. Notice that the finite part 
was chosen such that it is independent of the dimensional regularisation parameter ^u. The 
dependence on it is contained in the divergent part and in O(e) terms only. The divergent 
factors in (|109|) are governed by an evolution equation [158] |6"8"| 11591 116U|, [70] . For theories with 
vanishing (3 function, such as Af = 4 SYM El the factors containing the divergences take the 
particularly simple form 



div(x) = exp 



{-k£ j > 

i i=i 



rggp go 

(Ze) 2 le 



(110) 



where 

r cusp (a) = Y, a'r« p , G(a) = £ a'0« p . (Ill) 
I i 

Here r cusp is the cusp anomalous dimension (of Wilson loops), and G is the scheme-dependent 
collinear anomalous dimension. Their perturbative expansion in the coupling constant up to 
two loops is (in the DRED scheme) 



r cusp (a) = 2a - 2C 2 a 2 + 0(a 3 ) , G{a) = -( 3 a 2 + 0(a 3 ) . 



(112) 



We can see that ()108|) is indeed of the factorised form expected on general grounds, see formula 
(|109p , and the coefficient in front of the double pole agrees with the known expression (|112p at 
one-loop for the cusp anomalous dimension. 



By definition [361 E7J [38] , t ne cusp anomalous dimension r cusp (a) describes specific ultravi- 
olet divergences of a Wilson loop evaluated over a contour with a cusp. Its appearance in the 
infrared divergent part of the scattering amplitude (jllOp is not accidental. It has its roots in 
the deep relation between scattering amplitudes in gauge theory and Wilson loops evaluated 
over specific contours in Minkowski space-time, defined by the particle momenta |37 [ 1159} 175]. 
A distinguishing feature of this contour is that it has cusps. As was found in Refs. |37 ] 1159] 175], 
the infrared divergences of planar scattering amplitudes are in one-to-one correspondence with 
the cusp ultraviolet divergences of light-like Wilson loops. It should be mentioned that this 
relation is not specific to N = 4 SYM and it holds in any gauge theory, including QCD. We 

"This necessarily implies that a regulator is used where the (3 function is indeed zero 
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explain this relationship in a one-loop example in section 16.11 

The two-loop expression for r cusp (a) in a generic (supersymmetric) Yang-Mills theory was 
found in Refs. [3H1 IE]- We remark that in M = 4 SYM theory, r cusp (a) is known at weak 
coupling to four loops [ID], and it is conjectured to be governed to all loops by the BES equation 
[30] . At strong coupling, the solution of the latter produces a strong coupling expansion of 
Tcusp (°0 [441 [4"5l 06] . The first few terms of this expansion are in agreement with the existing 
quantum superstring calculation of Refs. [411 14"2| 143] . The non-universal collinear anomalous 
dimension G(a) is known to four loops at weak coupling [39, 161J and there exists a prediction 
at strong coupling [79] . 

For planar MHV amplitudes (where the common helicity structure can be factored out), the 
generalisation of (|1(J9|) to an arbitrary number n of gluons is 

M n = In M^ MHV ) = Z n + F n + O(e) , (113) 

with 

i « / r (0 r (l)\ n ( ,[2]\- 1 " 

where tf^ = s^j+i = (pi + Pi+i) 2 is the invariant mass of two adjacent gluons with indices % 
and i + 1, and the periodicity condition i + n = i is tacitly implied. The finite part F n = 
F n (t^) is a dimensionless function of the momentum invariants = (p» + • • .pj+j-i) 2 of the 
scattering process 0. From (|113f) we can see that MHV amplitudes contain two interesting 
pieces of information: the cusp anomalous dimension r cusp in the double pole term, and an a 
priori arbitrary function of the momentum invariants in the finite part. What makes the latter 
particularly interesting is that it turns out not to be arbitrary, but seems to have a very simple 
form, as we will discuss in section 13.1.41 In some sense, the finite part F n is more interesting 
than r cusp because at a given loop order, it is a function (in general of many variables), whereas 
Fcusp is merely one number. Finally, note that G is scheme dependent because one can always 
arbitrarily redefine the dimensional regularisation scale \x. For the same reason F n contains in 
general a scheme-dependent constant, but its functional form is scheme independent. 



3.1.3 Calculational technique 

The loop corrections to M4 can be represented in terms of scalar loop integrals, as for example 
the scalar box integral 1^ encountered in the previous section. How do these integrals arise? 
In a standard approach, one would write down all Feynman graphs contributing to a given 
amplitude, involving gluons, fermions, scalars, and also ghosts. Then, one has to work out the 
numerator algebra in the integrals corresponding to these graphs. Usually one also has to use 
integral reduction identities in order to express the amplitude in terms of master integrals [162J. 
One can think of the scalar integrals as the result of going through this (lengthy) procedure. 
With increasing loop order (or number of gluons) this approach involves writing down a large 
number of Feynman graphs and would be very cumbersome if not impossible. For example, even 
the tree level formula (|1U4|) for large n would require writing down an astronomical number of 
Feynman graphs, since their number grows factorially with n. 

This is why one often employs a different technique based on the unitarity of the 5*-matrix, 
see [163] and references therein. The latter, which describes the scattering of particles in a 

14 Note that F n does not depend on all independently, because (more than four) four-dimensional vectors 
have to satisfy Gram determinant constraints. 
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Figure 14: Graphical representation of the unitarity relation 2 Im (T) = TT^ for the four-gluon 
amplitude. The sum is over a complete set of intermediate states. 



quantum theory, can be written as 

S = l + iT. 

Requiring unitarity of the S'-matrix, S< = S^ 1 , one obtains 

2Im(T) = i(T f -T) = TT 1 " 



(115) 



(116) 



If one inserts a complete set of states between T and on the right-hand side of (|116p . one gets 
the diagrammatical relation shown in Fig. [T4l There, the sum is over all allowed intermediate 
states and the integral is over their on-shell momenta. 



Equation (|116|) relates objects at different orders in perturbation theory. For example, the 
imaginary part of a one-loop amplitude on its left-hand side is expressed through tree-level am- 
plitudes on its right-hand side. One can think of the r.h.s. as being a one-loop amplitude that 
was cut into two tree-level pieces. 

In the unitarity based method of Bern et al |52l [53] , one does not compute the integral on 
the right-hand side of (|116p directly. Rather, the idea is to reconstruct the amplitudes from their 
cuts. More precisely, one uses (1 11 6 [) in order to construct the loop integrand of integrals that 
have the correct behaviour under unitarity cuts. A generic one-loop amplitude can be expressed 
in terms of master integrals by integral reduction techniques [162J. Therefore it is sufficient to 
know the coefficients with which the master integrals contribute to the amplitude. Moreover, 
in supersymmetric field theories, it is known from integral reduction methods, combined with 
power counting arguments, that only a certain subclass of scalar integrals can appear in the 
amplitude [S3]. All of them have either s- or t-channel cuts, so they can be detected by the 
corresponding cuts. Therefore, one-loop amplitudes in M = 4 SYM are 'cut-constructible'. 

Let us give a one- loop example of the cutting technique, following [153] . Consider the four- 
gluon amplitude A^(l~ , 2 _ , 3 + , 4 + ) at one loop. From the s-channel cut, we have 

A4(l-,2- 3+,4+)| s _ cut = ^|^0 J 2vr5W(/?)2vr5W(/ 2 2 ) (117) 

xAT e (-h, 1~,2~, l 2 )AT e (-l 2 , 3+ 4+ h) 

Here the cut momenta are l\ = k and l 2 = k — p\ — p 2 , with k being the loop momentum. The 
sum in <\117\) is over all allowed particles and helicity states for the cut lines. In this example, due 
to supersymmetric Ward identities, the contribution of scalars and fermions to (|117p vanishes. 
Hence, the only non-vanishing contribution comes from gluon tree amplitudes with the following 
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helicities 



15 



j4.^ ree (— ^2 , 3 + , 4 + , ^ 



(-ill) (12) (2l 2 )(l 2 - h) 

(-hh) 4 
'H 2 3)(34)(4/ 1 )(Z 1 -Z 2 ) 



(118) 



We now replace the delta functions 27rj( + )(Zj) by full propagators i/lf, while still using the 
conditions if = in the numerator. In this way, we construct the integrand of an integral that 
has the correct behaviour in the s-cut, as was explained above. 



A 4 (l-,2-,3 + ,4+)| s _ cut 



d D k 1 1 



(12) / 



(2vr)^/f/ 2 2 ( _ /ll)(12) {2h) {h _ h 



(-hhY 



(-Z 2 3) (34) (4Zi) (h - l 2 ) 



(119) 



By construction, (|119l) has the correct properties in the s-channel. Since we want to arrive at 
usual loop integrals, we make propagators appear in the denominator by writing e.g. 



1 



\2h 



(2z 2 ) (k- Pl y 

Using identities like (|120l) . we see that we can factor out the tree amplitude, 

d D k [hi] (14) [4Zi] (iii 2 ) [i 2 3] (32) [2Z 2 ] (hh) 



(120) 



^4;i(l",2-,3 + ,4+)| 



s— cut 



-iA\ iee 



(121) 

After some algebra, the numerator in the integrand can be simplified to — st(k— Px) 2 (k>+Pi) 2 ■, up 
to a term involving the totally antisymmetric Levi-Civita tensor, which drops out of a four-point 
amplitude. Thus, the final answer is 



A 4; i(l",2-,3 + ,4^ 



-st A 



tree 



dPk 



(2vr) D k 2 (k - Pl ) 2 (k -px- P2) 2 {k+ Pi ) 2 



■ (122) 



The t-channel cut can be evaluated in a similar manner. There, supersymmetry allows scalars 
and fermions to contribute as well. Taking all contributions into account, one again finds the 
scalar box integral, just as in the s-channel. By consistency, the coefficients obtained in both 
channels have to match. Thus, we conclude that in Af = 4 SYM 0, 



A 4;1 (l- 2- 3+ 4+) = -stAT e f 
in agreement with (|106l) . 



d D k 



1 



(2vr) D k 2 (k - Pl ) 2 (k -px- p 2 ) 2 (k + p A ) 2 ' 



(123) 



One-loop amplitudes in supersymmetric theories were shown to be 'cut-constructible' in the 
above sense. We remark that in non-supersymmetric theories, there are rational terms in the 
amplitude that cannot be detected by cuts, and have to be determined by other means. This 
problem can be solved, however, by using L>-dimensional cuts. The reason is that in D di- 
mensions, the rational terms produce cuts that can be detected. This approach further uses a 
supersymmetric decomposition of one- loop amplitudes (|92p . This means that only cuts of scalar 

15 Note the change of helicity in Zi,2 depending on whether the lines are ingoing or outgoing 

16 Strictly speaking, the preceding calculation proves (|123p up to terms vanishing as O(e) only, since we used 

4-dimensional cuts. To extend the derivation to 0(e) terms one could use D-dimensional cuts. However, in [53] 

it was argued that (|123[) is valid to all orders in e (see also [55]). 
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loops have to be evaluated in D dimensions. These techniques have allowed the computation of 
several QCD amplitudes at one loop. 

In J\f = 4 SYM the cutting techniques were successfully applied to higher loop amplitudes as 
well. They allowed the computation of amplitudes that would have been practically impossible to 
determine on the basis of Feynman diagrams. We will summarise the results for MHV amplitudes 
in J\f = 4 SYM in the next section. Finally, we mention another important development, that 
of 'generalised cuts', where one requires more than two propagators to be on shell, see [164] . 



3.1.4 Results of perturbative calculations 

Over the last thirty years, considerable efforts have been undertaken to compute loop corrections 
to scattering amplitudes in QCD and its supersymmetric extensions, see e.g. the review [57] , 
Most progress has been made for MHV amplitudes in J\f = 4 SYM. Up to now, the following 
MHV amplitudes have been computed in N = 4 SYM: 

• four gluons at two loops [55j [56] and three loops [39] ; the four-gluon amplitude is known 
in terms of loop integrals at four loops [ID] , and the integrals were evaluated (numerically) 
up to the second pole term, which allowed to determine r cusp at four loops; a conjecture 
for the integrals appearing at five loops was put forward in [75] 

• five gluons to two loops |165|, [75] 

• six gluons to two loops [92] 

• n gluons to one loop [52] 

Similar computations in QCD are much harder and up to now are restricted to lower loop orders. 



In [56], Anastasiou, Bern, Dixon and Kosower (ABDK) noticed an interesting iterative struc- 
ture in the two-loop result for the four-gluon amplitude in N = 4 SYM. Based on [56] and the 
result of their three-loop calculation, Bern, Dixon and Smirnov (BDS) formulated a conjecture 
for the n-gluon amplitudes in M = 4 SYM to all orders in the coupling constant [39]. As we 
already discussed, the divergent part of the scattering amplitudes is well understood, so their 
conjecture is essentially a statement about the finite part. Their conjecture implies that the 
finite part F n , c.f. (|113l) has the same functional form as at one loop, Fn \ multiplied by the 
(coupling-dependent) cusp anomalous dimension, 



1 



Fn = ^r cuap (a)i^ 1 ) + const 



(124) 



For example, in application to the four- and five-gluon amplitudes this means that 



Fa 
F 5 



1, 



cusp 



(«)[ 



In 2 



+ const 



cusp(^) 



5> 

.1=1 



— ■ In h const 



s i+l,i+2 j+4 



(125) 
(126) 



where s,t and Sjj+i = (pi + Pj) 2 are the momentum invariants of the scattering processes. 



Why should the finite part of the amplitudes be so simple? One might speculate that some 
symmetry could be at the origin of this simplicity. In the next section, we will study the 
four-gluon amplitude in more detail. We will see that indeed one can observe 'dual' conformal 
properties of the integrals contributing to it, thus providing hints for a broken conformal sym- 
metry governing the gluon amplitudes. We will see in section [HJ that the conjectured duality with 
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Figure 15: One-loop box integral 1^ and two loop double box integral 1^ and their dual 
pictures. 

Wilson loops, if true, can explain these observations, and even make them more quantitative, 
by providing broken conformal Ward identities. 

We should stress that very recently, the BDS conjecture was found to fail at six gluons and 
two loops [92]. We will see in section [9] that it is modified in a way predicted by the duality 
§5§ between Wilson loops and gluon amplitudes. The duality, if true, also implies that the BDS 
conjecture is correct for four and five gluons to all loops (which is consistent with all available 
data) . 

3.2 Dual conformal properties of the four-gluon amplitude 

Let us focus for now on four-gluon amplitudes, where results from high-order calculations in 
perturbation theory are available. We will see in this section that the integrals contributing to 
the four-gluon amplitude are all from a set of 'pseudo-conformaF integrals, in the sense that will 
be defined presently. 

Let us come back to the scalar box integral (|107|) . An important observation made in [77J is 
that it has particular properties in a dual 'coordinate' space, defined by the change of variables 

p 1 = Xi - X 2 = X 12 , P2 = X 2 3 , P3 = X34 , Pi = X41 . (127) 

We stress that (]127p is simply a change of variables, and not a Fourier transform. Under the 
identification (|127p and that of k = X15, one obtains the expression 

' {1) M;e) = C / ^ ^ , 9/> , „ = C / 2 \ x 2 \ . (128) 



k 2 {k - pi) 2 (k - pi - P2) 2 (k +Pi) 2 J xf 5 xl 5 xl 5 xl 5 

The integral together with its dual picture is shown in Fig. [JSl The point is that in the dual 
'coordinate' space defined by (112T|> . the integral has (broken) conformal properties. Indeed, 
if it was not for the on-shell conditions pf = which lead us to use dimensional regularisation, 
would be given by the finite off-shell conformal integral <&W encountered in section 12.5.41 
Recall that $W is conformal because under a conformal inversion, the variation of the integral 
measure is precisely canceled by the variation of the four propagators going to the vertex. In 
D = 4 — 2e, we see that the cancelation is no longer exact, and hence the (dual) conformal 
symmetry is broken. So we conclude that 1^ has broken conformal properties in the dual co- 
ordinate space. 

In order to avoid confusion with (broken) conformal Ward identities for Wilson loops dis- 
cussed in later parts of this thesis, we will call integrals of this type 'pseudo-conformaF |4U| . 
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At this stage one might think that the appearance of the pseudo-conformal integral 1^ is 
just a one-loop coincidence. However, we will see presently that this pattern continues to hold 
to higher loops. At two loops, 



Mf > - T 



' s 2 tl^(s,t;e) + st 2 l^(t,s;e) , (129) 



one encounters one new integral, the double box integral 

J( 2) (at . e)=( «f d^pd^q 

Going to dual coordinates (c.f. Fig. [T5l) . we see that it is the dimensionally continued version 
of <I>( 2 ) considered in section [276l Notice the dotted line between external points one and three, 
which is there to make sure that all external points have the same conformal weight. There is 
an easy way to implement the change of variables (|127p by using the following pictorial rule: 

• for each loop integration in p-space, draw a dot inside the loop, which corresponds to an 
integration vertex in the dual coordinate space. 

• the external points in the dual coordinate space are situated between two consecutive 
external momenta, denote them by a point. 

• draw a straight line between points that are separated by a momentum line. 

The numerator lines have to be added by going through the change of variables (|127p . For 
pseudo-conformal integrals they are usually easy to find since they have to be chosen such that 
each integration vertex has conformal weight four, and from the convention that external points 
have conformal weight one. 

We remark that the notion of dual graphs exists for planar graphs only [166J. For this reason 
it is not clear how to apply this notion to non-planar corrections to gluon scattering amplitudes. 



At three loops, the four-gluon amplitude is given by two integrals (see Fig. [16 



M 4 



(3) 



1 



sHli 3)a (s, t; e) + st 3 li 3)a (t, s; e) + 2st 2 lf ] \s, t; e) + 2s 2 tlf )b (t, s; e) , (131) 



which are again both 'pseudo-conformal'. The two three- loop integrals appearing in the four- 
point amplitude (|13ip . and depicted in Fig. [TBI are 



r( 3 W„ /.^ _ n'A I d P d rd 1 



li 3)a (s,t;e) = C 3 J 



p 2 (p — hi) 2 (p — k\ — k^) 2 

1 



(p + r) 2 r 2 (q — r) 2 (r — k^ — k^) 2 q 2 (q — ki) 2 (q — ks — k^) 2 ' 

and 

j4-2e^ j4-2e^, j4-2e„ / , \2 



(132) 



(3) bf + . , n3 f d^pd^rd^q(p + r) 



lP b (s,t;e) = C*J 



p2 g2 r 2 ^ _ ^ 1 ^2 (jy _|_ r _ ^ 1 ~)2 

1 



{p + r - ki - k 2 ) 2 (p + r + fe 4 ) 2 (q - /c 4 ) 2 (r +p + q) 2 (p + qf 



. (133) 



The first integral is the triple box integral, which in the dual space corresponds to the di- 
mensionally continued version of $( 3 ). The other three- loop integral, the so-called 'tennis court', 
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or rather 'squash court' integral, exhibits a new interesting feature. Its expression contains a 
numerator (p + r) 2 (which is not depicted in the momentum space picture), given by x\ 5 in 
the dual space and depicted by a dotted line. We see that five propagators are connected to 
integration point five, which is one to many to cancel the conformal weight coming from the 
transformation of the measure under a conformal inversion. This extra weight is precisely can- 
celed by the presence of the numerator. 




with their dual pictures. 

At four loops, the number of pseudo-conformal integrals contributing to the amplitude in- 
creases compared to three loops, and one finds eight of them |40| . They are shown in Figs. 
1171181 along with their dual conformal pictures. Let us comment that the integrals shown in 
Fig. [TTJ were predicted by the rung rule [39] , whereas the integrals shown in Fig. [18] cannot be 
obtained from it. The rung rule came from the observation that when constructing (re + l)-loop 
integrals from re- loop integrals by the unitarity method, the two-particle cuts have an iterative 
structure, and (re + l)-loop integrals can be obtained by adding a 'rung' to re-loop integrals in all 
possible ways, (see also [167] ). It is also interesting to note that in the calculation of [40J, two 
integrals that one might (naively) call pseudo-conformal came with coefficient zero, i.e. they did 
not appear. A possible explanation for this was found in [83] , where it was shown that those 
latter integrals are not finite in four dimensions for off-shell gluons. It is striking that of all 
10 pseudo-conformal integrals, the two integrals that did not appear are both ill-defined in this 
sense. 

Based on the 'conformal hypothesis', i.e. that the expectation that the amplitude should be 
given by a linear combination of pseudo-conformal integrals, an amplitude in terms of integrals 
was conjectured at five loops |78j . Strikingly, it was found from unitarity cuts that integrals 
that are ill-defined in the sense discussed in the previous paragraph are absent in the five-gluon 
amplitude. 

We remark that, compared to a generic gauge theory, such as for example QCD, the number 
of integrals appearing in the final expression for the amplitudes is very small. It seems to be 
a special property of J\f = 4 SYM that the amplitudes can be reduced to such a small set of 
integrals. In particular, in calculations done so far, it was observed that no diagrams (represent- 
ing integrals) with triangle or bubble subgraphs appear If true in general, this seems to be a 
special feature of M = 4 SYM. We remark that if one demands or conjectures that the amplitude 
should be built from pseudo-conformal integrals only, then triangle and bubble subgraphs are 
automatically excluded. 

Based on [77J SO], in [83J it was conjectured that the four-gluon amplitude is a sum over 
17 In the four-loop computation [3U] this was used as an assumption. 
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(d) (e) (f) 

Figure 17: The rung-rule dual diagrams. A factor of st has been removed. Figure taken from 

m- 

pseudo-conformal integrals 

M 4 = 1 + al{1) oily (134) 
j 

where the sum runs over all pseudo-conformal integrals (which would be finite off-shell in four 
dimensions), I, a is the coupling constant and 1(1) is the loop order of a given integral. According 
to (I134p . the four-gluon amplitude would be fixed by the relative coefficients ax coming with 
the pseudo-conformal integrals. Can these coefficients be predicted as well? One constraint on 
them comes from the known exponential form of the divergences of the amplitude, c.f. (I109p . 

oo „ 

lnM 4 = -^a'-p^ + 0(l/e). (135) 
i=i ^ e ' 

The e expansion of individual integrals at a given loop order I typically contains divergences up 
to e~ 21 . This means that the coefficients in (|134p have to be such that when one takes the log- 
arithm, the leading divergences cancel, and at most double poles remain. This relates integrals 
appearing at different loop orders. A practical problem is that in order to use this relation, 
one has to be able to extract divergences of integrals. It would be highly desirable to have an 
efficient way of determining at least the leading singularities of a given integral in order to fix 
some of the relative coefficients (c.f. [1681 I169| 1167] for related work). 

Recall that r cusp , the coefficient of the double pole, is predicted from a (conjectured) inte- 
grable model [30]. This implies that so should be the coefficients ax in ()134p . since they implicitly 
determine r cusp . At the moment, integrability is mainly applied to determine the spectrum of 
anomalous dimensions. Solving the problem of determining the coefficients ax would be an 
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(d 2 ) (f 2 ) 

Figure 18: The non-rung rule dual diagrams. A factor of st has been removed. Figure taken 
from jU3]. 



important step towards extending the use of integrability to the computation of scattering am- 
plitudes. 

We studied four-gluon scattering amplitudes. For more than four gluons, one can also de- 
fine the notion of pseudo-conformal integrals, and the integrals appearing in five- and six-gluon 
MHV scattering amplitudes up to two loops seem to have this property (76j[92]. This is certainly 
worth further study. 

To summarise, the discovery of this hidden structure in (four-)gluon amplitudes in M = 4 
SYM is one of the results of this thesis. One can say that the 'conformal rule' (meaning that the 
amplitude is constructed from a linear combination of 'pseudo-conformal' integrals only) is an 
improvement over the 'rung-rule' [39], which fails starting from four loops @. The conformal rule 
works at least to four loops and possibly at five loops [78]. Will this new rule also fail at some 
loop order? We like to think that the 'conformal rule' is correct to all loops. Support for this 
expectation comes from the duality between Wilson loops and gluon amplitudes which will be 
discussed in sectional If this duality holds, then the gluon amplitudes inherit the dual conformal 
properties of the Wilson loops. This might explain why only 'pseudo-conformal' integrals appear 
in the amplitude. Moreover, it would fix (at least some of) the relative coefficients with which 
the integrals appear in the amplitude, since then, in addition to (|135l) . the (logarithm of the) 
sum of the integrals would have to satisfy the same conformal Ward identity as the Wilson loop 
(see section [8} . 



4 Wilson loops and gluon amplitudes at strong coupling 

In this section we summarise work by Alday and Maldacena |79l [HO] that motivated the idea of 
a duality between Wilson loops and gluon amplitudes in M = 4 SYM. 



4.1 Prescription for computing scattering amplitudes at strong coupling 

In [TjJ], Alday and Maldacena proposed a way to calculate scattering amplitudes in J\f = 4 SYM 
at strong coupling via the AdS/CFT correspondence [HI] . According to [79], gluon scattering 
at strong coupling is described by the scattering of open strings ending on a D-brane in AdS^ 
space with the metric 



ds 2 = R 2 



dx\ +1 + dz 2 



(136) 



It still predicts correctly the appearance of the rung-rule diagrams, together with theirs prefactors, but it 
misses other diagrams, as was explained above. 



38 



in 

Figure 19: Figure illustrating the minimal surface, whose area A enters Eq. (|139H . 



The .D-brane is located at z = zm, where ztr can be thought of as an infrared regulator. The 
scattering takes place at fixed angle and at large momentum, and one finds that the solution is 
dominated by a saddle point of the classical action (see also [170 J for a similar computation in 
flat space). 

It turns out that the boundary conditions are easier to describe after changing variables to 
T-dual coordinates y^, for which the boundary conditions are simply 

Ay^ = 2vr^ . (137) 

Notice that this is precisely the relation between momenta of gluons and the dual coordinate 
space variables 11271 from section [3T2l Importantly, after defining r = R 2 /z, one finds again an 
AdS$ metric in the T-dual space, 



ds 2 = R 2 



dyndy^ + dr z 



(138) 



So, in summary, the computation of the gluon amplitudes in the original space was related to 
a computation in a dual coordinate space, which also has an AdS§ metric. The new boundary 
conditions imply that one should compute a minimal E surface that ends on a polygon defined 
by the coordinates y^ (which are in turn defined by the momenta of the gluons through (|137p ) , 
as shown in Fig. [TTJ1 

R 2 . \ ( V\ 



A oc exp ^"^r^minj = exp ^"^r^minj • (139) 

The proportionality symbol in (|139|) stands for the dependence of the amplitude A on the he- 
licity information of the scattering process. 



It was further pointed out [7S] that the computation of the minimal surface in (|139|) is 
formally similar to the computation of a light-like Wilson loop (defined in the dual AdS$ space 
with coordinates y^) at strong coupling. However, from the AdS (and strong coupling) point of 
view, one might expect the relation between gluon scattering amplitudes and Wilson loops to 
break down at subleading orders in l/\/A- 



4.2 Computation of 4-cusp Wilson loop at strong coupling 

The computation of the minimal surface in (|139p is a non-trivial problem and so far has been 
achieved explicitly only in the case of n = 4 gluons. Following [79J, we first review the one- 
cusp solution obtained in [82], and then present the solution for four gluons. The quantities 

19 Strictly speaking, one should say extremal, see a discussion in [171] . 
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Figure 20: Graphical representation of the one-cusp solution (|144|) . Figure from [79J. 



entering in (|139|) are IR divergent, and hence have to be regularised. It turns out that for our 
purposes it is sufficient to compute the minimal surface for the unregularised string action, and 
then compute the regularised minimal area by using this unregularised solution. Therefore we 
present here only the solution of the unregularised minimal surface. 

4.2.1 Single cusp 

Let us consider Wilson lines forming a single cusp on the boundary of AdS$. It is sufficient to 
consider a subspace AdS^ of AdS§, with the Wilson lines extending along the light-cone, 

r = 0, yi = ±jto, y >0- (140) 

The metric is 

dh = -d 2 yo + dy 1 + d *r (141) 

The Nambu-Goto string action is 

S = ^j d Sl ds 2 ^det (d a X»d b X»G^). (142) 
In a nonparametric form, where r is expressed through yo and yi, r = r(yo,yx), it reads 

S = 27 / d yodyi^\/l + (dor) 2 + (d 1 r) 2 . (143) 
It is straightforward to see that 



r(y , yi ) = y/2y/ y*-yl (144) 

solves the Euler-Lagrange equations following from (|143p with the boundary conditions (|140p . 
The minimal surface defined by (j!44H is shown in Fig. [20j 



4.2.2 Rectangular Wilson loop 

Let us now turn to the rectangular Wilson loop. In the case s = t, the minimal surface has to 
end on a square on the boundary. By scale invariance, the size of the square can be chosen to 
be 1, and hence the boundary conditions are 

r(±l,y 2 ) = r(j/ 1 ,±l)=0, y (±l,y 2 ) = ±y 2 , y ( yi , ±1) = ±y x . (145) 
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Inspired by the one-cusp solution, Alday and Maldacena found that 



3/0(2/1,2/2) = 2/12/2, r(y 1 ,y 2 ) = yj (1 - y\){l - y\) , (146) 

satisfies the equations of motion following from (|142p and solves the boundary condition (|145p . 
By applying conformal transformations on the solution (I146p . one can find the solution for 
arbitrary s and t \79 \ 1172] . For its explicit expression and more details, we refer the reader to 
[79]. 



4.2.3 Analog of dimensional regularisation at strong coupling 



A subtlety arises when one tries to compute the value of the action (|142p by plugging e.g. 
the solution (|146f) into it. One finds that the action diverges. Therefore the action has to 
be regularised. In order to be able to compare to gauge theory results, which are usually 
derived in dimensional regularisation, Alday and Maldacena introduced an analog of dimensional 
regularisation at strong coupling. Their regularisation consists in changing the metric (|138p to 



ds 1 



DCD 



dy 2 + dr" 2 

r 2+e 



(147) 



with Ad = A 



(4-7re" 



and c D = 2 4e vr 3 T(2 + e) [73]. This leads to the modified action 
V A_dcd f £<e=o 



S 



2tt 



(148) 



where £ t= o is the Lagrangian density for AdS§, and e < regularised the infrared divergences 
O The Euler-Lagrange equations following from (|148p now depend explicitly on e. So, strictly 
speaking, the solutions of the Euler-Lagrange equations following from ()142p obtained in the 
previous section are no longer valid. In particular, since conformal invariance is broken by the 
regulator, conformal arguments cannot be used to relate the solution at s = t to that of arbitrary 
s and t. 



On the other hand, one can imagine the solutions obtained in this way are still useful if one 
is only interested in the first few terms of the perturbative expansion in e of the minimal area. 
Indeed, Alday and Maldacena argue that one get the accurate value of the minimal surface by 
plugging the unregularised solution following from (|142|) into the regularised action, neglecting 
terms of order 0(e). More precisely, since the divergences come from the cusps, Alday and 
Maldacena modify the solution near the cusps by an 0{e) term, and in this way they obtain the 
correct contribution to the finite part of (|139jl . Since this constant is scheme-dependent, we will 
not be interested in it here. For a more detailed discussion, see [79] EE 

One finds \79\ 



2ir 



' 2vr 



A div + - In 2 



+ C 



(149) 



where C is a scheme-dependent constant. The divergent part is given by A^ v = 2Adi V ,s+2A c n Vit) 
where 



A 



div,s 



-e- 2 - e- 1 



1 - In 2 



-6/2 



(150) 



20 Note that in [173] it was found that the dimensional regularisation procedure used in [7S] has to be modified 
at subleading orders in 1/vA. 

21 See also [174] for a related discussion of regularisation prescriptions for Wilson loops at strong coupling, when 
using a cut-off. 
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Together with the prescription (11390 . we see that (j!49H is in agreement with the general form of 
the divergences of gluon scattering amplitudes, c.f. equation (|114|) , It also reproduces correctly 
the known strong-coupling value of r cusp , 

r cusp (a) = V2v^ + 0(l,Vv^). with a = ^ = A. (151) 

As was already explained in section [3j the subleading e^ 1 term and the constant in the finite 
part of (I149p are scheme-dependent, and therefore we are only interested in the functional form 
of finite part. Most importantly, it is in agreement with the BDS ansatz for the four-gluon 
amplitude at strong coupling! This constitutes strong evidence that the BDS ansatz holds for 
the four-gluon amplitude. 

The following comments are in order: Extending the result (|149p to n > 5 points turns 
out to be difficult, and so far no explicit solution for the minimal surface is known. In a later 
paper, Alday and Maldacena found that the problem simplifies when the number of gluons 
n is very large. By considering a special contour they were able to relate this problem to a 
known spacelike contour, and found a disagreement with the BDS ansatz for n — > 00 [80J. This 
suggests that the BDS ansatz should break down at a certain loop level and for a certain number 
of gluons. Finally, note that the strong-coupling calculation of [79] is insensitive to the helicity 
configuration of the scattering amplitude under consideration. 



5 Wilson loops 

Wilson loops were introduced in [175] in the context of quark confinement in strongly coupled 
QCD. For a Yang- Mills theory with gauge group SU(N), they are defined as 

W{C) = 1 (0| Tr P exp (ig £ da? A A |0> . (152) 

where A^(x) = A^(x)t a is the gauge field, t a are the SU(N) generators in the fundamental 
representation normalised as ti(t a t b ) = ^S ab , and P indicates the ordering of the SU(N) indices 
along the integration contour C n . W(C) is a gauge invariant functional of the contour C, along 
which the gauge field A% is integrated. 



5.1 Renormalisation properties 

5.1.1 Wilson loops denned on smooth contours 

Let us summarise the renormalisation properties of the Wilson loops (|152|) . In a generic gauge 
theory with ultraviolet divergences, the gauge field A% is renormalised by multiplicative coun- 
terterms, and so is the coupling constant g. These counterterms do not depend on the contour 
C. Let us assume that this standard renormalisation has already been carried out and focus on 
the divergences inherent to the Wilson loop (|152p . 

For smooth contours, the only divergence of (|152p is linear (in the cut-off) [361 1176] . and 
proportional to the length of the contour. It can be absorbed in an overall factor, 

W{C) = e - KL ^ x finite , (153) 

which can be interpreted as the mass renormalisation of a test particle moving along the contour 
C [361 EE]. In (fTSHD . L{C) is the length of the contour C. 
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For example, the one-loop correction to (|152l) is given by rn 




oc 



C JC 



dx^dx'v- 
(x — x') 2 ' 



(154) 



where the double line in the picture denotes the integration contour C, and the wiggly line is a 
free gluon propagator. The integral in (|154f) has a divergence at x = x', which can be regularised 
by introducing a cut-off a -1 , 



oc 



dx^dx'^ 



cJc - x') 2 + a 2 



i(s) • x(s + i) dsdi 
[x(s + i) -x(s)] 2 + a 2 



(155) 



Following [36], we choose a parametrisation of the contour with x(s) ■ x(s) = const, such that 
x(s) ■ x(s) = 0. Since the divergence in (|155f) comes from the integration near t = 0, we can 
write 



Wg(C) 



ds x 2 (s) 



1 



+ finite 



'-A x 2 (s)t 2 + a 2 

= - / dsy/x 2 ^) + finite = -L(C) + finite. (156) 
a J a 

For smooth contours, there are no divergences apart from the linear divergence |176l 1177] . 




Figure 21: An example of a Wilson loop defined on a contour C having a cusp at the point x, the cusp 
angle being 7. 



5.1.2 Contours with cusps 

The situation gets more interesting when the contour C is not smooth, but contains one or 
several cusps, as shown in Fig. [2TJ In that case, UV (short-distance) divergences associated to 
each cusp appear [36 ] 1178] . For example, a one-loop calculation (for gauge group SU(N)) gives 
the result 

W(C) = 1 - 2g 2 C F [ 7 cot( 7 ) - 1] In (^\ , (157) 

where L is the length of the contour, a is a short-distance cut-off, and 7 is the cusp angle. Cp = 
t a t a = (N 2 — 1)/(2N) is the quadratic Casimir of SU(N) in the fundamental representation. 
The result of |178j is that to any order in the coupling constant, the new divergences can be 
removed by a multiplicative renormalisation, 

W R {C) = Z(^)W{C) . (158) 

The divergences of Wr{C) depend only locally on the contour C, through the cusp angle 7. The 
locality of the counterterms persists when there are several cusps with cusp angles 7^. In that 
case, the multiplicative renormalisation factor 

Z( 7 i,..., 7n ) = Z( 7 i) ••• Z( 7n ), (159) 
22 In (|154l) and the following we neglect a trivial normalisation factor . 
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factorises into a product of renormalisation factors, one for each cusp. Thus, there are no 
'anomalous' divergences that depend non-locally on the contour Q 

In |18U[ [3"B] a renormalisation group equation for Wilson loops with cusps was derived. It 
reads 

^ajl + ^ 9R) dg^ + rcusp(7 ' 9r) ) Wr ( l ^ 7 ' W ' 5r) = • (160) 

Here <?r is the renormalised coupling and 



d 



r (7>5R) = \rm e ^ Q Zn—Z 1 • 
From (|16ip one can derive a renormalisation group equation for Z , 



d 
dgR 



r(7,5R,e) 



Z-\ 7 ,g R ,e) = 0. 



Its solution is 



Z 1 (T: 5R) e ) = exp 



9R 



(161) 



(162) 



(163) 



In = 4 SYM in dimension D = 4 — 2e, and when using a scheme where 0(ffR, e) = — <?rc, the 
integral in (1163|) can be done explicitly, 



Z(7,SU,e) = exp 



^ 2n e 



.71=1 



(164) 



where we expanded r(7, <?r, e) = X^^o fl ,2n r*- n ^ (7, e) • One can define a Z-factor in analogue with 
the MS-scheme by requiring that r(7,gR,, e) = T(j,ffn). So we have seen that the renormalisa- 
tion group equation (|162p implies the exponentiation of the divergent part of the Wilson loop, 
and in N = 4 SYM it takes the particular simple form (|164p . 



The Wilson loops studied in this report are defined in Minkowski space, with cusp having 
an infinite cusp angle. These will be studied in more detail in section 15.41 



We mention that additional divergences can appear if one considers instead of a closed con- 
tour C an open line. In that case one looses gauge independence, and there are logarithmic 
divergences associated to the endpoints |176| . 

We also mention the interesting case of crossed contours (i.e., with self- intersections) |178} 
I181j . For these, there are additional divergences associated to the self-intersection. It is found 
that the loop with the self-intersection can mix with the product of two similar loops [178J, 
which differ from the initial one by a different ordering of the colour indices. 

5.2 Wilson loops in the AdS/CFT correspondence 

The study of Wilson loops in the AdS/CFT correspondence started with the paper [H] (see also 
[182J). There, Maldacena showed how to compute a particular kind of Wilson loop in N = 4 
SYM at strong coupling, using the AdS/CFT correspondence. The Maldacena- Wilson loop is 
similar to (|%|), but comprises an additional term involving scalars, 

W{C) = ^(0|TrP exp (ig j> ds [x^A^ + |0> . (165) 

23 This can be also understood in the formalism of |179l 1177] . where the nonlocal Wilson loop is reformulated 
in terms local fermionic degrees of freedom. 
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The 6 1 in (fT65l) parametrise 5 5 of AdS 5 x 5 5 , i.e. 6> 7 6>j = 1. 



Depending on the shape of the contour, the Maldacena- Wilson loop can be invariant under 
a certain number of supersymmetries. For example, half of the super symmetries are preserved 
if i M = const, i.e. if the contour C is a straight linecj. Such Wilson loops were studied in [183J. 
Perturbative calculations up to order g 4 and a strong coupling calculation using AdS/CFT 
suggest that 

W(C) straight line = 1 , (166) 

i.e. that the loop corrections vanish identically. This is in close analogy with two- and three- 
point functions of 1/2-BPS operators, see section [2~5"1 which are also invariant under half of the 
supersymmetry. For Wilson loops that are invariant under less super symmetry, see e.g. [184] . 

In this report we will study Wilson loops with light-like contours, i.e. \x\ = 0. In this case 
(|165jl reduces to the usual definition of the Wilson loop (J2|). Therefore we will not expand on 
the discussion of the Maldacena- Wilson loop. For reviews, see [185] 1186] and also section 10 of 

[TBU. 

5.3 Loop equations 

Another development concerning Wilson loops is the study of loop equations. These describe 
the behaviour of the loops under small deformations of the contour C. We refer the interested 
reader to the literature [188] 189J. 

5.4 Light-like Wilson loops 
5.4.1 Definitions 

In this report, we consider M = 4 SYM theory with SU(N) gauge group in Minkowski space. 
The central object of our consideration is the light-like Wilson loop defined as 

W{C n ) = ^(0\TrPexp(i£ da^A M (x)J |0) , (167) 

where C n is a polygon with n cusps. The n segments of this polygon C n = UILi ^ joining the 
cusp points xf (with i = 1, 2, . . . , n) 

U = {x»(Ti) = nxf + (1 - Ti)x1 +l \ n e [0, 1]} , (168) 

are all light-like, i.e. xf i+l = 0. Such Wilson loops were considered for the first time in [85] for 
n = 4 points in the special kinematics xf 3 = —x\±. 

5A.2 Cusp singularities 

We already saw in section [5.1. 2l that cusps cause specific ultraviolet divergences (UV) to appear 
in Wilson loops. The fact that the cusp edges are light-like makes these divergences even more 
severe [85]. In order to get some insight into the structure of these divergences, we start by giv- 
ing a one-loop example. Later on in this section we give arguments which lead to the all-order 
form of the divergences. 

One-loop example 



One can imagine that this Wilson 'loop' is closed at infinity. 
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Xi-1 



Figure 22: The Feynman diagram contributing to the one-loop divergence at the cusp point Xi. The 
double line depicts the integration contour C n , the wiggly line the gluon propagator. 



Let us first discuss the origin of the cusp singularities in the n-gonal Wilson loop W(C n ) to 
the lowest order in the coupling constant. According to definition (I167p . it is given by a double 
contour integral, 

W(C n ) = 1 + \{ig) 2 C F I dx» I dy» D^{x-y)+0(g A ). (169) 

Here Cp = (N 2 — 1)/(2N) is the quadratic Casimir of SU(N) in the fundamental representation 
and D^ v {x — y) is the gluon propagator in the coordinate representation 

{A«(x)A b ,(0))=g 2 5 ab D^(x). (170) 

To regularise the ultraviolet divergences of the integrals entering (|169|) . we use dimensional 
regularisation with D = 4 — 2e and e > 0. Also, making use of the gauge invariance of the 
Wilson loop (|167|) and for the sake of simplicity, we perform the calculation in the Feynman 
gauge, where the gluon propagator is given by | 
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D^(x)= V ^D(x), D(x) = -Iii-^(-x 2 + iO)- 1+£ (^V 7B ) e - (171) 

The divergences in ()169p originate from the integration of the position of the gluon in the vicinity 
of a light-like cusp in the Feynman diagram shown in Fig. [22] (the other Feynman diagrams at 
one loop vanish or are finite, see section [6T3]) . The calculation of this diagram is straightforward 
|85j . Let us compute it for a gluon attached to the edges pi := Xi + \^ and Pi-\ := Xii— i, as 
shown in Fig. 



\ 9 2 C F , _ 7 2 ^rvi \ f 1 dxdy{pi-p., 

V[pi-i,Pi) = — ^2~(~ e P ) r ( 1_e ) J 



i-i 



o [{piX+pi-xy) 



211— c 



9 2 CF ( , ^,^(l-6 



H 2 xUl,i-i) -^TT 2 - ■ ( 172 ) 



4vr 2 V /- -t+X,t-lJ 2e 2 

In (|172[) we have dropped the — iO prescription in the gluon propagator because we are only 
interested in the kinematical region where all i _ 2 < 0. Adding together the contributions 
of all cusps we obtain the following one-loop expression for the divergent part of W(C n ): 

W{C n ) = 1 + ^_C F I -— 2 M-M+i V 2 Y + 0(e°) | + 0( 5 4 ) , (173) 



25 As in |86| . we redefine the conventional dimensional regularization scale as /j, 2 ire lB i— > /x 2 to avoid dealing 
with factors involving ir and the Euler constant 7b- 
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where the periodicity condition 



x ij — x i+n,j ~ x i,j+n — \ x i x j) K^'^J 

is tacitly implied. We see from (|173p that the leading divergence in W(C n ) is a double pole in 
e, compared to the single pole in (|164p . This is due to the light-likeness of the edges defining 
C n , which creates additional divergences. We will discuss the general structure of divergences 
of light-like Wilson loops in the next section. 



All-loop structure 

Generalising (|173|) to higher loops we find that the cusp singularities appear in W(C n ) at 
the l-th. loop-order as poles (afi 2e ) 1 /e m with m < 21. Furthermore, W(C n ) can be split into a 
divergent ('renormalisation') factor Z n and a finite ('renormalised') factor 

\nW{C n ) = Z n + F^ . (175) 

From the studies of renormalisation properties of light-like Wilson loops it is known [85, 190J 
that cusp singularities exponentiate to all loops and, as a consequence, Z n has the special form 

E 

^ = - iX>T (-4W) fc + — ) . (176) 

with a = g 2 N/(87r 2 ). Here Tcusp and are the expansion coefficients of the cusp anomalous 
dimension and the so-called collinear anomalous dimension, respectively, 

2 

r c us P (a) = Y, a 1 r« p = 2a - \a 2 + 0(a 3 ) , (177) 
l>i 6 

T(a) = al r( ° = - 7 <3a 2 + 0(a 3 ) . 
l>\ 

We have already seen at one- loop order that the divergences in W(C n ) are due to the presence of 
the cusps on the integration contour C n . They occur when the gluon propagator in Fig. [22] slides 
along the contour towards a given cusp point X{. The divergences arise when the propagator 
becomes singular. This happens either when a gluon propagates in the vicinity of the cusp point 
(short distance divergences), or when a gluon propagates along the light-like segment adjacent 
to the cusp point (collinear divergences). In these two regimes we encounter, respectively, a 
double and single pole. 



Going to higher orders in the coupling expansion of W(C n ), we immediately realise that the 
structure of divergences coming from each individual diagram becomes much more complicated. 
Nevertheless, the divergences originate from the same part of the 'phase space' as at one-loop: 
from short distances in the vicinity of cusps and from propagation along light-like edges of the 
contour. The main difficulty in analysing these divergences is due to the fact that in diagrams 
with several gluons attached to different segments of C n various regimes could be realized simul- 
taneously, thus enhancing the strength of poles in e. This implies, in particular, that individual 
diagrams could generate higher order poles to W(C n ). 



26 For convenience of the presentation, the following formulae are given in the planar limit, with a = g^s". For 
general N, the only difference is that T cusp and F in ()177|) depend on iV as well as on a. 
27 Formula (|176[) follows from the evolution equation (8) of [85] . 
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The simplest way to understand the form of In Z n in (I176|) is to analyse the divergences of 
the Feynman diagrams not in the Feynman gauge, but in the axial gauge |191j defined as 

n-A(x) = 0, (178) 

with being an arbitrary vector, n 2 ^ 0. The reason for this is that the contribution from 
individual Feynman diagrams become less singular in the axial gauge and, most importantly, 
the potentially divergent graphs have a much simpler topology [1921 1193] r 8 l The axial gauge 



gluon propagator in the momentum representation is given by the following expression, 

5 -« = -ra ^ , < t > = ^- fe T^ £ + ^(4- (179) 



The polarisation tensor satisfies the relation 



t u d^{k) = {D-2) + ^, (180) 



(to be compared with the corresponding relation in the Feynman gauge, which is g^ v d (k) = D) 
from which we deduce that for k 2 = 0, it describes only the physical polarisations of the on-shell 
gluon. El 

Let us consider a graph in which a gluon is attached to the i-ih segment. In configuration 
space, the corresponding effective vertex is described by the contour integral J dripf A^(xi — 
PiTi). In the momentum representation, the same vertex reads J drip^ A^(k) e lk ^ Xi ~ pm ^ where 
the field A^(k) describes all possible polarisations (2 longitudinal and D — 2 transverse) of the 
gluon with momentum k^. Let us examine the collinear regime, when the gluon propagates along 
the light-like direction p^ . The fact that the gluon momentum is collinear, k^ ~ pf , implies 
that it propagates close to the light-cone and, therefore, has a small virtuality k 2 . Furthermore, 
since pfA^(k) ~ k^A^{k) we conclude that the contribution of the transverse polarization of 
the gluon is suppressed as compared to the longitudinal ones. In other words, the most singular 
contribution in the collinear limit comes from the longitudinal components of the gauge field 
A^(k). The properties of the latter depend on the gauge, however. To see this, let us examine 
the form of the emission vertex in the Feynman and in the axial gauge. In the underlying 
Feynman integral, the gauge field A^(k) will be replaced by the propagator D^ik), with v 
being the polarisation index at the vertex to which the gluon is attached. In this way, we find 
that 

hP-A (k) i ( k u , Feynman gauge 

pf ZV(fc) ~ k^ w (k) = = -p x { k2 fe _ ^ J q1 _ (181) 



(fcn) 2 (fen) 



axial gauge 



Since k 2 — > in the collinear limit, we conclude that the vertex is suppressed in the axial gauge, 
as compared to the Feynman gauge. This is in perfect agreement with our physical intuition - the 
propagation of longitudinal polarisations of a gluon with momentum k^ is suppressed for k 2 — > 0. 
This property is rather general and it holds not only for gluons attached to the integration con- 
tour, but also for the 'genuine' interaction vertices of the N = 4 SYM Lagrangian [192| 1193"] . It 
should not be surprising now that the collinear divergences in the light-like Wilson loop come 
from graphs of very special topology that we shall explain in a moment. 



28 Note that the contribution of individual Feynman diagrams is gauge dependent and it is only their total sum 
that is gauge invariant. 



_><( 



That is the reason why the axial gauge is called physical. 
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(a) (b) (c) 

Figure 23: Maximally non-Abelian Feynman diagrams of different topology ('webs') contributing to 
In W(C n ). In the axial gauge, the vertex-type diagram (a) generates simple pole; the self-energy type 
diagram (b) generates double pole in e; the diagram (c) with gluons attached to three and more segments 
is finite. 



Another piece of information that will be extensively used in our analysis comes from the 
non-Abelian exponentiation property of Wilson loops [176^ 11941 1195] . It follows from the combi- 
natorial properties of the path-ordered exponential and it is not sensitive to the particular form 
of the Lagrangian of the underlying gauge theory. For an arbitrary integration contour C it can 
be formulated as follows: 



(w(c))=i+y: (h) wik) 
k=i v / 



exp 



.k=l 



4tt 2 J 



^ I — - - I <""'ir 



(182) 



Here denote the perturbative corrections to the Wilson loop, while c^w^ are given by 

the contribution to from 'webs' w^ k ' with the 'maximally non-Abelian' color factor c^ k \ 

To the first few orders, k = 1,2,3, the maximally non-Abelian color factor takes the form 



-(*) 



CpN k 1 , but starting from k = 4 loops it is not expressible in terms of simple Casimir 



operators. Naively, one can think of the 'webs' w^ k ' as of Feynman diagrams with maximally 
interconnected gluon lines. For the precise definition of 'webs' we refer the interested reader to 



Let us now return to the analysis of the cusp divergences of the light-like Wilson loops and 
take advantage of both the axial gauge and the exponentiation (|182p . A characteristic feature of 
the 'webs' following from their maximal non-Abelian nature is that the corresponding Feynman 
integrals have 'maximally complicated' momentum loop flow, e.g. they cannot be factorised 
into a product of integrals. When applied to the light-like Wilson loop, this has the following 
remarkable consequences in the axial gauge: 

• the 'webs' do not contain nested divergent subgraphs; 

• each web produces a double pole in e at most; 

• the divergent contribution only comes from 'webs' localised at the cusp points as shown 
in Fig.E|a) and (b). 

These properties imply that the divergent part of lniy(C n ) is given by a sum over the cusp 
points with each cusp producing a double and single pole contribution. Moreover, the cor- 
responding residue depend on i+1 - the only kinematical invariant that one can built out 
of three vectors xf_ 1 , and xf +1 satisfying xf_ li = xf i+1 = 0. In this way, we arrive at the 
known relation (|176p for the divergent part of the light-like Wilson loop. 
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Our consideration relied on the analysis of Feynman diagrams at weak coupling. It was 
recently shown in Refs. [S2\ 11961 H97| that the structure of divergences of In W(C n ) remains the 
same even at strong coupling. 



6 Duality between Wilson loops and gluon amplitudes 

In this section we formulate and discuss the main point of this thesis - the proposed duality 
relation between planar MHV amplitudes, .M^ HV \ and light-like Wilson loops, W(C n ). As 
was already mentioned in section [3] it is known that the leading infrared (IR) divergences of 
gluon amplitudes are in direct equivalence with the leading ultraviolet divergences of Wilson 
loops. We will illustrate this relationship by a one-loop example in section 16.11 In section 16.21 
we state the duality relation between gluon amplitudes and Wilson loops. 



6.1 IR divergences and their relation to Wilson loops 

Infrared divergences of scattering amplitudes can be understood in terms of Wilson loops. For 
non-Abelian gauge theories, this property was found in [198, 199, 37J. It is known [68], [70] that 
the IR divergent part of planar scattering amplitudes factorises into a product of form factors, 
and therefore it is sufficient to understand the relationship between IR divergences of form fac- 
tors and a particular kind of Wilson lines. 

Here, we will give a one-loop example, following to some extent [73] Let us take for 
simplicity a form factor of massive quarks, with on-shell momenta p± and — p 2 , p\= p\= m 2 . 
At tree level, this form factor is simply given by 



n = fabcV2(p2)l a u(p X ). (183) 

At one loop, IR divergences come from the vertex diagram shown in Fig. Oonly. Its contribution 




Figure 24: Contribution to the quark form factor at one loop. 

is given by 

V a _ f i_ r [ dDk vjPih^ifo- h + m)T{h- h + m)Yu{pi) ~ ( , - 

1 - Jabc 2 J (2?r)D ^ _ k)2 _ m2+ . Q] [(pa _ k)2 _ m2+ . Q] JV W > U» 4 i 

where D^ik) is the momentum space gluon propagator, e.g. D^ u {k) = 7]^ /k 2 in the Feynman 
gauge. The divergence in (|184p originates from the soft region of loop momenta, k <C m. In 
dimensional regularisation, it manifests itself as a pole in e, f^d D k/k 4: ~ 1/(4 — D), where S 
denotes the integration region for small k^\ In that region, we can neglect jk in the numerator, 
and make the approximation 

(Pi ~ k) 2 - m 2 ps -2(px • k) , (185) 



30 I am indebted to Gregory Korchemsky for enlightening explanations, on which this section is based. 
31 Note that for massless quarks, m 2 = 0, there are additional collinear divergences, which is why the leading 
pole is double in that case. 
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and similarly for the other quark propagator. Using in addition the Dirac equations of motion 
for the on-shell quarks, e.g. — m)u(pi) = 0, we obtain 



rr = f^C F ^ 2h °u( Pll /^p^j^W). (186) 



We can interpret the integral in (I186p as coming from the vertex diagram 1241 with the quark- 
antiquark-gluon vertex replaced by the eikonal vertex 



P " (187) 




—pi ■ k + iO 



Notice that the quark helicity is not modified by this eikonal vertex, and r" IR is proportional 
to the tree level form factor in (|183|) , Prom now on, we will drop the tree-level factor, i.e. we 
consider 

T _ 9 2 r f d D k y\ p» - 

Note that from the invariance of (|188l) under rescalings p\ — ► \\p\ and p 2 — > X2P2, one can de- 
duce that it depends on p^ and p 1 ^ through the dimensionless variable cosh 7 = (p± ■ P2)l \Jv\v\ 
only. Geometrically, 7 corresponds to the cusp angle formed by the momenta p\ and p 2 in 
Minkowski space. 

Now we would like to extend the integration region from £ to the full phase space. In 
doing so, the integral in (|188p acquires additional, artificial ultraviolet (UV) divergences. In 
dimensional regularisation, the UV and IR divergences cancel, and the integral over the full 
phase space is zero, 

T Cf I W^dI TT^TT 4r^ Z V^ = AlR + AUV = - (189) 

2 J (2ir) u (-pi ■ k + iO) {-p2 ■ k + iO) 

So instead of studying the IR divergences of (|189p at k <C m, we can equivalently study its UV 
divergences at k 3> m. These UV divergences are well-understood, because they correspond to 
cusp divergences of Wilson loops, as we will see presently. By using the identity 



1 

-1 



00 

dsie i S1 (- Pl -k+io) ; ( 19Q ) 



-pi ■ k + i0 j 
and similarly for —p2 ■ k + iO, we can rewrite ()189p as 

/ = -{igfc F y o d S1 ds 2 J -^p^0^y H ~ sipi ~ S2P2+tO) 

j /"OO /"OO 

= o(^) 2 Cf / dsi / ds 2 PiP2 D^ v (pisi +P2S2) , (191) 
1 Jo Jo 

where in the second line we introduced the Fourier transform of the gluon propagator, 

D,u{x) = J £^D^(ky Hx+i0) ■ (192) 
Equation (|191|) is just the one-loop contribution to a Wilson loop El 

W(C) = (0|TrPexp (ig j dx^{x)^j |0) , (193) 
32 Strictly speaking, it is a Wilson line, but one can imagine it is closed at infinity. 
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defined by the contour C = C% U C2, where C\ = {sp\,s £ [— 00, 0]} and C2 = {tp2,t £ 
[0, 00]}. This means that the quarks were effectively replaced by Wilson lines along their classical 
trajectory. 




(a) 



(6) 



Figure 25: One-loop contribution to cusped Wilson loop. The double lines denote the integration 
contour. The UV divergences of the Wilson loop correspond to the IR divergences of the form factor. 

After the Fourier transform, the UV divergences of (|193p come from the region where pis\ + 
P2S2 "C m, i.e. from the integration region of small s\,S2- This corresponds to the situation 
where the gluon propagator is integrated near the cusp of the Wilson loop. Taking into account 
the relation between UV and IR divergences (|189f) . we see that the IR divergences of the form 
factor (|184p are equivalently described by the UV divergences of a cusped Wilson loop, as shown 
in Fig. [25] (a). The latter satisfies a renormalisation group equation [180] (see also section l5.1.2p . 
As explained in section 15.1.21 for N = 4 SYM this RG equation implies that the divergent part 
of the Wilson loop (and hence that of the form factor as well) takes the form 



For an explicit two-loop calculation, see [200J. 

We remark that since we need only the UV divergences of the Wilson loop shown in Fig. 
(a), which arise from the integration of gluons near the cusp, we can choose a modified integra- 
tion contour, as shown in Fig. [23 (b). In this way, the new Wilson loop has only the desired 
UV divergences, but no IR divergences. 

Another important remark is that the definition of the Wilson loop (I193P contains an unusual 
mixture of momenta and coordinates. Indeed, the 'coordinates' of the Wilson loop (|193p are 
defined by the gluon momenta (and formally have 'wrong' canonical dimension). The same is 
true for the Wilson loops that enter the duality with gluon amplitudes, as we will see in the next 
section. These Wilson loops 'live' in a dual coordinate space defined by the gluon momenta. 

6.2 Duality relation 

The conjectured duality states that in the planar M = 4 SYM theory the finite part of the 
logarithms of the gluon amplitude and of the Wilson loop are equal (up to an inessential additive 



upon the formal identification of the external on-shell gluon momenta in the amplitude with the 
light-like segments forming the closed polygon C n (the contour of the Wilson loop), 




(194) 



constant) 



F (MHV) =F? (WL) +const 



(195) 



a a a 

: = X i+1 ~ X i ■ 



(196) 
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Thus, the Mandelstam variables for the scattering amplitudes t[ = (pj + . . . £>j + j_i) 2 are related 
to the distances z?- between two cusp points on the integration contour of W(C n ) are follows, 

4^/4? := x i,i+j/ x k,k+l ■ (197) 

The divergent parts of the scattering amplitudes and the light-like Wilson loops are also related 
to each other but the relationship is more subtle since the two objects are defined in two dif- 
ferent schemes (infrared regularisation for the amplitudes and ultraviolet regularisation for the 
Wilson loops), both based on dimensional regularisation. From the previous section (see also the 
more detailed discussion in [§T]) we know that the leading IR divergence of the amplitude (the 
coefficient of the double pole in equation (|114p ) coincides with the leading UV divergence of 
the Wilson loop (the coefficient of the double pole e^y in equation (|176p ). since both are con- 
trolled by the universal cusp anomalous dimension r cusp (a). One can also achieve a matching 
of the coefficients of the subleading simple poles corresponding to the (non-universal) collinear 
anomalous dimensions G(a) and T(a). |£j To this end one relates the parameters of the two 
different renormalization schemes as follows, 

^Wuv := t? ] /p&erW , e UV := -e m e**W . (198) 

Here the functions 7(a) and 5(a) are chosen in a way to compensate the mismatch between G(a) 
and r(a), without creating extra /^-dependent finite terms |91j . It should be stressed that this 
procedure is not analogous to comparing two different renormalisation schemes for the computa- 
tion of the same divergent object. It is rather a change of variables (regularisation parameters) 
which allows us to compare two different objects computed in two different schemes. Another 
such change of variables is the identification (|196p of the particle momenta with the light-like 
segments on the Wilson loop contour. 

The duality relation (|195p was inspired by the prescription of Alday and Maldacena for 
computing gluon scattering amplitudes at strong coupling [79]. This prescription essentially 
recasts the amplitudes into light-like Wilson loops in the dual variables (|196|) . as discussed in 
section HI A priori, one would expect the strong coupling relation between gluon amplitudes and 
Wilson loops to receive 1/y/X corrections, which might spoil the relation at weak coupling [79] . 
Nevertheless, in |83| it was found that at one loop and for four points the Wilson loop and the 
gluon amplitude agree, which lead to the idea that the duality might also be true perturbatively. 
In [S3j the duality at one loop was shown to apply to n-point amplitudes as well. These one-loop 
results will be reviewed in section 16.31 

6.3 Duality at one loop 

We already discussed results of perturbative calculations for gluon amplitudes, in section 13.1.41 
Here, we review the Wilson loop calculation at one loop. Let us start by computing for n = 4 
points, following [83J. The perturbative expansion of the Wilson loop 

W{C A ) = ^(0|Tr P exp (i j> dx^A^x^j |0) (199) 

to the lowest order in the coupling is given by 

W(C 4 ) = l + \(ig) 2 C F I dx^l dy» D^(x-y) + 0(g i ), (200) 
1 JCi Jc 4 

33 We thank Paul Heslop for turning our attention to the incomplete discussion of this point in the first version 
of [91| . We are also grateful to Lance Dixon for a discussion of the different physical interpretations of the IR and 
UV simple poles [201] . 
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where D^ v is the free gluon propagator and C F = (N 2 — 1)/(2N C ). Its expression in the 
Feynman gauge was given in (|171|) . Recall that the contour C4 consists of four light-like segments 
x ii+l> x 1i+i = (with the cyclicity condition i + 4 = i for the indices). Hence, the Wilson 
loop can only depend on the two available Lorentz invariant quantities, x 2 % and x^- Through 
(|196|) they correspond to the Mandelstam variables s and t of the four-gluon amplitude, i.e. 



13 



(pi + P2) = s and a^ 4 = (P2 + P3) = 






(c) 



Figure 26: The Feynman diagrams representation of the integrals contributing to (|200p . The double line 
depicts the integration contour C and the wiggly line the gluon propagator. Figure adapted from [83 . 



The Feynman diagrams contribution to (I200p are shown in Fig. [26j In the Feynman gauge, 
diagram (c) vanishes due to the light-likeness of the edge the gluon is attached to. The vertex 
type diagram (a) was already calculated in section 15.4.21 In application to the four-point case 
we find from (|173p 

g 2 C F T(l - e) 



V(pi ,P2) + V(p 2 ,P3) + V(p 3 ,p 4 ) + V(p 4 ,Pl) 



4vr 2 



[(" 



~1E 



[X s 



-IE 



(201) 

Let us now turn to diagram (b), in which the gluon is attached to segments P2 and p 4 . This 
diagram is finite for e — > 0, so we evaluate it in D = 4 dimensions 



Integration yields 
I{P2,Pa) 



9 2 C F 



dx dy {p2 ■ Pi) 



g 2 c F 



4tt 2 Jo (P2X + P4,y + P3) 2 8ir 2 J tx + sy - (s + t)xy 



dx dy (s + t) 



(202) 



9 2 C F 
8tt 2 



dx 



x 



s 
s+t 



, s . 1 ~ 

in — h In 

t x 



9 2 C F 

167T 2 



[ln 2 (s/t)+vr 2 ] , (203) 



where we assumed that s, t < 0. The other diagram of the same topology, with the gluon going 
from segment p\ to p% gives the same contribution, since it is obtained from interchanging s and 
t. Finally, summing up all contributions and expanding in powers of e, we find 



1 + 



9 2 C F 
4vr 2 



{~2 K-^nT + (-/^n + \ m 2 ( J) + y } ( 204 ) 

In order to compare with the corresponding gluon amplitude, we write the finite part of In W{C 4 ) 
at one loop in the planar limit, which amounts to replacing C F — > N/2, i.e. g 2 C F /(Ai: 2 ) — > a, 



\nW(C A ) = F A 



(WL) 



2 



13 



L '21 



+ y 



+ 0(a 2 ) 



(205) 



Comparing to the corresponding expression for the four-gluon amplitude (|125p and the one-loop 
value of the cusp anomalous dimension (|177p . we see that (I205p is indeed in agreement with the 
duality relation (j!95j) . Let us stress again that the additive constant in the finite part of (|204|) 
is scheme dependent. 



To get the n-point result, one has to evaluate the finite diagram (b) in more general kine- 
matics |84| . The result is again in agreement with the duality relation (|195|) , 
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6.4 Checks of the duality at two loops and beyond 



We remark that the one-loop Wilson loop calculations presented in the previous section did not 
use special properties of N = 4 SYM. Indeed, only the free gluon propagator D fMU (x) entered 
the calculation, and hence the result in other gauge theories would have been the same. 

This motivates to extend the checks of the duality (11950 beyond one loop, where one can 
see a difference between e.g. QCD and N = 4 SYM. The remaining sections of this thesis are 
devoted to this. We report on results that confirm the duality (I195p in several non-trivial cases 
at two loops and beyond. 

In [86], Drummond, Korchemsky, Sokatchev and the present author carried out the two-loop 
computation of the Wilson loop for n = 4 points, which we will present in section Technical 
details of this calculation can be found in appendix [Bj where the evaluation of each Feynman 
diagram is shown explicitly. In section [8] we present the results of the publication [87] , where we 
derive all-order Ward identities for the Wilson loops (jl]) . In section [9j we perform a calculation 
at six points and two loops and compare to the BDS ansatz and a recently available two loop 
calculation of the six-gluon MHV amplitude. 

7 Two loop tests of the duality 

7.1 Rectangular Wilson loop 

In this section we present the result of the calculation of the rectangular Wilson loop 



to two loops. The details of the calculation can be found in the appendix. 

In order to compute the two-loop Feynman diagrams shown in Fig. [57] we employ the tech- 
nique described in detail in Refs. [1581 1159} 85, 190J. Furthermore, the rectangular light-like 
Wilson loop under consideration has been already calculated to two loops in |85} 1190] in the 
so-called forward limit x^ 2 = ~ x u^ or equivalently x\ z = — x| 4 , in which the contour C4 takes 
the form of a rhombus. 

The vertex-like diagrams shown in Fig. [27Tc). [27Td) and 127(f) only depend on the distance 
and hence we can take the results from [85 } 1190] . The Feynman diagrams shown in Fig. 127(h) 
-12711) are proportional to scalar products {x\i -2:34) and/or (223 -£14) and, therefore, they van- 
ish in the limit x^ 2 = — £34 due to x? i+1 = 0. For generic forms of C4, these diagrams have to 
be calculated anew. Similarly, we also have to reexamine the contribution of diagrams shown 
in Figs. [2"7Tel and 127(g). Finally, the calculation in Ref. [85|, 1190] was performed within the 
conventional dimensional regularisation (DREG) scheme. To preserve supersymmetry, we have 
to use instead the dimensional reduction (DRED) scheme. The change of scheme only affects 
the diagrams shown in Figs. l27Tc) and !27( k) which involve an internal gluon loop [16] , 

The results of our calculation can be summarised as follows. Thanks to the non-Abelian 
exponentiation (|182p . the two-loop expression for the (unrenormalised) light-like Wilson loop 
can be represented as 




(206) 



( 




) 



2 



hxW(C A ) = i^C F wW + 



C F Nw i2) + 0{g 6 ) . 



(207) 
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(i) (j) 00 (1) 

Figure 27: The Feynman diagrams contributing to lnPF(C4) to two loops. The double line depicts 
the integration contour C4, the wiggly line the gluon propagator and the blob the one-loop polarization 
operator. 



According to (|204p the one-loop correction u/ 1 ) is given by 

«,« = -1 [(-x? 3 + (-4, ^Y] + \ In 2 ( J) + y + 0(6) . (208) 

The two-loop correction is given by a sum over the individual diagrams shown in Fig. [27] 
plus crossing symmetric diagrams. It is convenient to expand their contributions in powers of 
1/e and separate the UV divergent and finite parts as follows 

« W =£[(-4iM 2 )* + (-^mT] {^2+iAg + ^+^}+4 >+0( e ), (209) 

where the sum goes over the two- loop Feynman diagrams shown in Fig. 1271(c)— (1) ■ Here 
(with < n < 4) are dimensionless functions of the ratio of distances x\^jx\^. Making use of 
(|209p . we can parametrise the contribution of each individual diagram to the Wilson loop by 

the set of coefficient functions . 

• UV divergent 0(l/e 4 ) terms only come from the two Feynman diagrams shown in Figs. 127(d) 
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andETtf) 



A 



(d) 



1 

16 ' 



A 



(f) 



1 

16 



(210) 



UV divergent 0(l/e 3 ) terms only come from the two Feynman diagrams shown in Figs,l27lfc) 
andETtf) 

*3=5. 4'» =-i (211) 



UV divergent 0(l/e 2 ) terms only come from the Feynman diagrams shown in Figs. [2"?Tc) 

mis) 



i(«0 



4' 
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7T 
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(e) 



7T 

24 ' 
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(f) 

-2 



4 + 96* ' 



A 



(g) 



7T 

48 

(212) 



UV divergent 0(l/e 1 ) terms come from the Feynman diagrams shown in Figs. l27T c) 
EHh)J23k) andE^l) 
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~2 ~ " 


A« = 




A?} = 


^M 2 , 


A« = 


\m 1 , 


A« = 


-\ Ml 



48 + 24^ 3 ' 



:M 2 



(213) 



Finite O(e ) terms come from all Feynman diagrams shown in Figs. l27( c)- f2"7T l) 



A 



(c) 



A, 



(e) 



A, 



(g) 



A, 



(i) 



1 00 



7T 4 1 A 

2 + - + -C3, 
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24 



M\ - -M 3 + vr 4 
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1 

24 



M 2 , 



M 1 + -M 2 , 
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A, 



(f) 



A, 
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A, 



0) 



7T , 

2880 ' 

12 2880 6 s 
^-M 2 + |m 3 + ^-Afi , 



-iM 2 
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7T 2 1 1 

-Mi H Mi M 2 M 3 

24 2 8 



(214) 



Here the notation was introduced for the integrals Mj = Mi(xf 3 / X24) 



Mi 

M 2 
M 3 



dp 


H 


'a/3 


P-a 






dp 


In | 


'a/3 


P-a 




v^a/3 


dp 


In | 


'a/3 


P-a 







1 

~~ ~2 

Hpp) , 

ln 2 (/3/3) 



vr 2 + In 2 



13 



21 



(215) 



with (3 = 1 — p., a = 1 — a and a/a = xl^/x^- We do not need the explicit expressions for 
the integrals M 2 and M3 since, as we will see shortly, the contributions proportional to M 2 and 
M3 cancel in the sum of all diagrams (for completeness, their explicit form can be found in 
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Appendix IB, 3D . Note that the integrals (I215p vanish in the forward limit x 2 3 = —xi, A . 

We would like to stress that the above results were obtained in the Feynman gauge. Despite 
the fact that the contribution of each individual Feynman diagram to the light-like Wilson loop 
(or equivalently, the A^— functions) is gauge-dependent, their sum is gauge- invariant. 

Next, we substitute the obtained expressions for the coefficient functions, Eqs. (|210p - (|214p . 
into (I209p and finally arrive at the following remarkably simple expression for the two-loop 
correction, 

W ^ = [(-4, ^ + (-4, m [e~ 2 ^ + e" 1 ^} - YA In 2 (J) - ^ + 0(e) . (216) 

We verify that in the forward limit, i.e. for x 2 3 = —x 2 ^, this relation is in agreement with the 
previous calculations of Refs. |85l 1190] . The following comments are in order. 

Arriving at (|216|) we notice that the leading UV divergent 0(l/e 4 ) and 0(l/e 3 ) terms cancel 
in the sum of all diagrams, in agreement with (1176f) . According to (I212p . the coefficients in 
front of 1/e 2 are given by a sum of a rational number and it 2 — term. The rational terms cancel 
in the sum of all diagrams. As a consequence, the residue of the double pole in e of in 
Eq. (|216p is proportional to £2- In a similar manner, the residue of at the single pole in e 
is proportional to £3 and this comes about as the result of a cancelation between various terms 
in (|213p containing rational numbers, it 2 — terms as well as the integrals Mi and M2. The most 
striking simplifications occur in the sum of finite O(e ) terms (I214D . We find that the integrals 
M2, M3, M 2 as well as the rational corrections and the terms proportional to ir 2 and £3 cancel 
in the sum of all diagrams leading to — 755 tt 4 + ^tt 2 M\. 

We would like to stress that the two-loop expression (|216p satisfies the "maximal transcen- 
dentality principle" in M = 4 SYM [136] . Let us assign transcendentality n to functions Li n P 4 ! 
and transcendentality 1 to a single pole 1/e. Then it is easy to see from (|216p that the coefficient 
in front of l/e n (including the finite 0(e ) term) has transcendentality equal to 4 — n. In this 
way, each term in the two-loop expression (|216p has transcendentality 4. For the same reason, 
the one-loop correction to the Wilson loop, Eq. (I208P has transcendentality 2. Generalising this 
remarkable property to higher loops in planar M = 4 SYM, we expect that the perturbative 
correction to the Wilson loop (|167l) to order 0(g 2n ) should have transcendentality 2n. 

Notice that in our calculation of the two-loop Wilson loop we did not rely on the multi- 
colour limit. In fact, due to the special form of the maximally non-Abelian colour factors, 
c( n ) = CpN n ~ l , relation (|2U7|) is exact for arbitrary N. As was already mentioned, these colour 
factors become more complicated starting from n = 4 loops, where we should expect terms 
subleading in N. 

7.2 Pentagonal Wilson loop 

In [87], the calculation of the pentagonal Wilson loop at two loops was carried out in order 
to further test the duality (|195p . The Feynman diagrams are similar to those shown in Fig. 
[271 Individual diagrams can depend on ratios of the five Lorentz invariants available from 
five light-like distances, x 2 3 , x| 5 , x\ x and x\ 2 - We derived (multiple) parameter integral 
representations for all diagrams. The pole part of all diagrams could be evaluated analytically, 

34 E.g. ln(a;) = Lii(l - a;) and <2 = tt 2 /6 = Li 2 (l) are assigned transcendentality 1 and 2, respectively. 



and for the finite part we found within good numerical precision the following result: 



r -. 9 1 T 1 5 2^2 „2 4 

(2) 1 7T 1 7 ^ 2 2\2e , ^ i /' X i,i+2 \ n / X i+l,i+3\ V" 

We would like to mention that the calculation of the pole part in [87] was carried out analytically 
for arbitrary n at two loops, in agreement with the general formula (|176H . As explained in 
that paper, one can define 'auxiliary' diagrams to subtract the pole terms, thereby providing 
a definition of the finite part of the two-loop n-cusp Wilson loop in terms of a set of finite 
parametric integrals. 



7.3 Check of the duality at two loops 

Let us introduce the 
in the planar limit: 



Let us introduce the following notation for the perturbative expansion of In W(C n ) = Z n +F^ L ^ 



F^=aFiY 1 h) +a 2 F%V + .... (218) 
Then we see from (I216p and (J2T7J that we have 



Comparing (pl"9l) . (f220l) to (fl"25j) . (fl"26j) it is straightforward to see that the duality relation (fT95|) 
is verified for n = 4, 5 at two loops. 



One may wonder, independently of the duality with gluon amplitudes, why the two- loop 
results for (the finite part of the logarithm of) the rectangular and pentagonal Wilson loops 
(j219[) and (|220p are so simple. For example, one could have expected nontrivial multi-variable 
functions to appear in the 5-point case, such as harmonic poly logarithms. In section [8] we will 
see that the simplicity observed here is explained by conformal symmetry. There, we will derive 
all-order broken conformal Ward identities that have stringent implications for F^ L \ 



8 Conformal symmetry of light-like Wilson loops 

In this section we review the derivation of all-order broken conformal Ward identities for the 
Wilson loops W(C n ) [86l [87] . 

Recall that C n is a polygonal contour with n light-like edges. The decisive observation is that 
such contours are stable under conformal transformations. That is, the contour C' n , which is the 
image of the contour C n upon a conformal transformation, is also made of n light-like segments 
with new cusp points x'^ obtained as the images of the old ones xf. This property is rather 
obvious for translations, rotations and dilatations, but we have to check it for special conformal 
transformations. Performing a conformal inversion P^l. x^' = x^/x 2 , of all points belonging to 
the segment 

x^r i )=r i x^ + (l-r i )xf +1 , (221) 
with r £ [0, 1], we obtain another segment of the same type, 

x'^tI)=tIx'? + (1-t>)x% 1 , (222) 

35 A special conformal transformation (boost) is equivalent to an inversion followed by a translation and then 
another inversion (cf. section [2]). 
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with (x' l+1 - x[Y = and r[ = n /[n + (1 - ^(xtf /(x' i+1 ) 2 ]. 

Were the Wilson loop W(C n ) well defined in D = 4 dimensional Minkowski space-time, then 
it would enjoy the (super)conformal invariance of the underlying M = 4 SYM theory. More 
precisely, in the absence of conformal anomalies we would conclude that 

W{C' n ) =W(C n ). (223) 

However, as we already pointed out in section 15.41 the light-like cusps of C n cause specific 
ultraviolet divergences to appear in the Wilson loop (|167|) , For this reason we use dimensional 
regularisation with D = 4 — 2e and e > 0, which breaks the conformal invariance of the action, 
as we will see presently. In the dimensionally regularised N = 4 SYM theory, the Wilson loop 
W(C n ) = (W n ) is given by a functional integral 



(W n ) = [ VAVXVcft e^--^ Tr Pexp (if dx ■ A(x) 



(224) 



where the integration goes over gauge fields, A, gaugino, A, and scalars, (ft, with the action 

S e = 2 2e J d D xC(x), C = TV [— 5-^1/] + gaugino + scalars + gauge fixing + ghosts. (225) 

Here fi is the regularisation scale and we redefined all fields in such a way that g does not 
appear inside the Lagrangian C(x). This allows us to keep the canonical dimension of all fields, 
in particular of the gauge field A^(x), and hence to preserve the conformal invariance of the 
path-ordered exponential entering the functional integral in (|224p . However, due to the change 
of dimension of the measure f dPx in (|225|) the action S e is not invariant under dilatations and 
conformal boosts, which yields an anomalous contribution to the Ward identities. 

8.1 Anomalous conformal Ward identities 

The conformal Ward identities for the light-like Wilson loop W(C n ) can be derived following 
the standard method [202], 11341 198j . by acting on both sides of (I225|) with generators of confor- 
mal transformations.. The expressions of the generators of conformal transformations acting on 
fields are given in equation (fTT 



Let us start with the dilatations and perform a change of variables in the functional integral 
(|224p . (j)'j(x) = <f)](x) + el])<fti(x). This change of variables could be compensated by a coordinate 
transformation x^' = (1 — s)x^. We recall that the path-ordered exponential is invariant under 
dilatations, whereas the Lagrangian is covariant with canonical weight Ac = 4. However, the 
measure f d D x with D = 4 — 2e does not match the weight of the Lagrangian, which results in 
a non- vanishing variation of the action S e , 

SnS £ = [ d D x C(x) . (226) 

#V e J 

This variation generates an operator insertion into the expectation value, (5oS e W n ), and yields 
an anomalous term in the action of the dilatation generator on (W n ), i.e., 

n lie f 

MW n ) = V (Xi ■ di) {W n ) = / d D x {C{x)W n ) . (227) 

In a similar manner, the anomalous special conformal (or conformal boost) Ward identity is 
derived by performing transformations generated by the operator K w , Eq. (|17p . on both sides 
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of (I224p . In this case the nonvanishing variation of the action 5kv S t again comes from the 
mismatch of the conformal weights of the Lagrangian and of the measure f dPx. This 
amounts to considering just the term in (fTTl) with A& = Ac — D = 2e, and hence 



n a ■ r 

K"(W n ) = V(2<^ • di - a$d?){W n ) = -- ^- / d D x x v (C(x)W n ) . (228) 
~1 9 V J 



The relations (|227|) and (|228|) can be rewritten for In (W n ) as 



MW.-^j^:^, (229) 



To make use of these relations we have to evaluate the ratio (C(x)W n ) / (W n ) obtained by in- 
serting the Lagrangian into the Wilson loop expectation value. Due to the presence of e on the 
right-hand side of (|229p . it is sufficient to know its divergent part only. 

8.2 Dilatation Ward identity 

As we will now show, the dilatation Ward identity can be derived by dimensional arguments 
and this provides a consistency condition for the right-hand side of (|229H . By definition (I224p . 
the dimensionally regularised light-like Wilson loop (W n ) is a dimensionless scalar function of 
the cusp points x\ and, as a consequence, it satisfies the relation 



\J2( x i ■ Q i) - I HW n ) = . (230) 

In addition, its perturbative expansion is expressed in powers of the coupling g 2 [i 2e and, there- 
fore, 

^(W n ) = 2eg 2 ^(W n ) = --^ j d D x (C(x)W n ) , (231) 

where the last relation follows from (|224|) . Recall that according to (|175|) . the Wilson loop (W n ) 
can be split into the product of divergent and finite parts 

In (W n ) = Z n + F n . (232) 

Notice that the definition of the divergent part is ambiguous as one can always add to Z n a 
term finite for e — > 0. Our definition (I176p is similar to the conventional MS scheme with the 
only difference that we choose the expansion parameter to be a^ 2e instead of a. The reason 
for this is that Z n satisfies, in our scheme, the same relation fid n Z n = 2eg 2 d g 2Z n as (W n ), see 
(|23ip . Together with (|232[) and (|230p , this implies that the finite part of the Wilson loop does 
not depend on the renormalisation scale, i.e. 

fid^Fn = 0(e) . (233) 

Using (|232p and the explicit form of Z n in (|176p . the relation (|23ip leads to the following 
dilatation Ward identity 37 ! 

n 

Y / (x i -d Xi )F n = 0. (234) 
i=i 



36 Another source of non-invariance of the action S e is the gauge-fixing term (and the associated ghost term 
of the non-Abelian theory) which is not conformally invariant even in four dimensions. However, due to gauge 
invariance of the Wilson loop, such anomalous terms do not appear on the right-hand side of 

37 In what follows, we shall systematically neglect corrections to F n vanishing as e — > 0. 
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Xi-l 



(a) (b) (c) 

Figure 28: The Feynman diagrams contributing to (C(x)W n ) to the lowest order in the coupling. The 
double line depicts the integration contour C n , the wiggly line the gluon propagator and the blob the 
insertion point. 



Adding to this the obvious requirement of Poincare invariance, we conclude that the finite part 
F n of the light-like Wilson loop can depend on the dimensionless ratios xfj/x^ only. In partic- 
ular, for n = 4 there is only one independent ratio, i.e. F± = F4 (^13/^24)- 



Furthermore, making use of (|232f) . (11760 and (I234p we find that the all-loop dilatation Ward 
identity for (W n ) takes the form 

n n / (Q \ 

Bln(W n ) = J> • di) HW n ) = - ± 5>' £ M-M+i^ 2 ) + • (235) 

i=l 1>1 i=l V / 

This relation provides a constraint on the form of the Lagrangian insertion on the right-hand 
side of (p29l) . 

8.3 One-loop calculation of the anomaly 

The derivation of the dilatation Ward identity (|234|) relied on the known structure of cusp sin- 
gularities ()176l) of the light-like Wilson loop and did not require a detailed knowledge of the 
properties of Lagrangian insertion (C(x)W n ) / (W n ) . This is not the case anymore for the special 
conformal Ward identity. 



To start with, let us perform an explicit one-loop computation of (C(x)W n ) / (W n ) . To the 
lowest order in the coupling, we substitute (W n ) = 1 + 0(g 2 ) in the denominator and retain 
inside C{x) and W n terms quadratic in gauge field only. The result is 



{C(x)W n 

(W n ) 



1 

'iiV 



Tr [{d,A v {x) - d u A a {x)f] Tr \{i§ c dy ■ A(y)f ) + 0(g 



(236) 



The Wick contractions between gauge fields coming from the Lagrangian and the path-ordered 
exponential yield a product of two gluon propagators (11701) . each connecting the point x with 
an arbitrary point y on the integration contour C n . Gauge invariance allows us to choose e.g. 
the Feynman gauge, in which the gluon propagator is given by (|17ip . To the lowest order in 
the coupling constant, the right-hand side of (|236j) receives non- vanishing contributions from 
Feynman diagrams of three different topologies shown in Figs. [281( a) - (c). 



In the Feynman gauge, only the vertex-like diagram shown in Fig. [281(a) develops poles in e. 
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Performing the calculation, we find after some algebra 38 ! 

2 2i 2e {C tl Wn) = (-^Ui + i^y\^ D Hx - Xi ) + e^){x\ Xl . lM+1 ) + 0(e c 
9 A 4 \"n/ i=1 L 



(237) 

where a = g 2 N/(87r 2 ) and the notation was introduced for 



TT^ ^ (x| , Xi , Xj-i-i) 



(is 
s 



5 {D) (x -xi- sxi-i-i) + 5 (D) (:r -Xi + sx iti+1 ) - 25 {D \x - x t ) 



(238) 

We see that the leading double-pole singularities are localised at the cusp points x = x%. The 
subleading single poles are still localised on the contour, but they are 'smeared' along the light- 
like edges adjacent to the cusp. 

Substitution of (l237jl into (f229|) yields 

n 

Bln(W n ) = ~aJ2 {-xti^T^ 1 + 0(a 2 ), (239) 

i=l 

n 

K"ln(W„) = -2a^xr(-x|_ lii+lM 2 )V 1 + 0(a 2 ) . 



i=l 



Notice that X^ILi ^^H^l^-i' x «> does not contribute to the right-hand sides of these rela- 
tions by virtue of 



d D xT^{x\x i . 1 ,Xi,Xi +1 ) =0, (240) 

J d D xx u T ( - r) (x\xi-i,Xi,x i+ i) = (xi-i + x i+1 - 2xi) u . 

As was already mentioned, the right-hand sides of the Ward identities (1239 p are different from 
zero due to the fact that the light-like Wilson loop has cusp singularities. 

We verify with the help of (|177p that to the lowest order in the coupling, the first relation 
in (|239p is in agreement with (I235p . 

8.4 Structure of the anomaly to all loops 

To extend the analysis of the special conformal Ward identity to all loops we examine the all- 
loop structure of divergences of (C(x)W n ) / (W n ) . They arise in a very similar way to those of 
the Wilson loop itself, which were discussed in section 15.4.21 

It is convenient to couple C{x) to an auxiliary 'source' J(x) and rewrite the insertion of the 
Lagrangian into (W n ) as a functional derivative 



(C(x)W n )/(W n ) = -i-j—\n{W n )j 
dJ(xj 



(241) 

j=o 



where the subscript J indicates that the expectation value is taken in the N = 4 SYM theory 
with the additional term J d D xJ{x)C{x) added to the action. This generates new interaction 



38 It is advantageous to perform the calculation by taking a Fourier transform with respect to x and later take 
the inverse transform of the final result. 
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vertices inside Feynman diagrams for {W n )j, but does not affect the non-Abelian exponentiation 
property (|182p l 39 l The only difference compared to (|182|) is that the webs w^ k ' now depend on 
J(x) through a new interaction vertex. Making use of the non-Abelian exponentiation we obtain 



(jC(x)W n )/(W n ) 



k=i v y 



. 5w^ 
''8J(x) 



(242) 



j=o 



As for ln(W n ), the webs produce at most double poles in e, and are localised at a given cusp. 
This allows us to write down a general expression for the divergent part of the Lagrangian 
insertion, 



2ie (C(x)W n 



SV £ (Wn) 



(243) 



1>1 i=l 



1 / r 



(0 

cusp 



2 \ le 



+ r« 5 (D) (x - Xi ) + TW(x; ar<, x m ) I + O(e) 



which generalises (]237p to all loops. The following comments are in order. 

In Eq. (|243|) . the term proportional to 8^ D '{x — Xi) comes from the double pole contribution 
to the web w^ k \ which is indeed located at short distances in the vicinity of the cusp xtf. The 
residue of the simple pole in the right-hand side of (|243|) is given by the cusp anomalous di- 



mension. This can be shown by substituting (j243H into the dilatation Ward identity (|227f) and 
comparing with 



The contact nature of the leading singularity in (|243p can also be understood in the following 
way. The correlator on the left-hand side can be viewed as a conformal (n + l)-point function 
with the Lagrangian at one point and the rest corresponding to the cusps. In the N = 4 SYM 
theory the Lagrangian belongs to the protected stress-tensor multiplet, therefore it has a fixed 
conformal dimension four. The Wilson loop itself, if it were not divergent, would be conformally 
invariant. This means that the n cusp points can be regarded as having vanishing conformal 
weights. Of course, the presence of divergences might make the conformal properties anomalous. 
However, the conformal behaviour of the leading singularity in (|243|) cannot be corrected by an 
anomaly! 40 ! Then we can argue that the only function of space-time points, which has conformal 
weight four at one point and zero at all other points, is the linear combination of delta functions 
appearing in (j243H 3 

As was already emphasised, the residue of the simple pole of 5w^ / 8J(x) at the cusp point 
Xi depends on its position and at most on its two nearest neighbours Xi-± and Xj + i, as well as 
on the insertion point. It gives rise to a function (x;Xi-i,Xi,Xi+i), which is the same for all 
cusp points due to the cyclic symmetry of the Wilson loop. 



Notice that in (|243p we have chosen to separate the terms with the collinear anomalous 
dimension from the rest of the finite terms. This choice has implications for the function 
T"\x;Xi-i,Xi,Xi + x). Substituting the known factor Z n (I176p and the particular form of (|243H 



39 We recall that the non-Abelian exponentiation is not sensitive to the form of the action. 

40 Such an anomalous contribution should come from a 1/e 3 pole in the correlator with two insertions of the 
Lagrangian, but repeated use of the argument above shows that the order of the poles does not increase with the 
number of insertions. 

41 The mismatch of the conformal weight four of the Lagrangian and D = 4 — 2e of the delta functions does not 
affect the leading singularity in 
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into the dilatation Ward identity (I227D . we derive 



/ d D x rW(x;x i - 1 ,x i ,x i+1 ) = . (244) 
i=i ^ 

We can argue that in fact each term in this sum vanishes. Indeed, each term in the sum is a 
Poincare invariant dimensionless function of three points Xi and x%+\. Given the light- 

like separation of the neighbouring points, the only available invariant is xf_ li+1 . Further, 
Xi—i, Xi, Xi+i) cannot depend on the regularisation scale because n always comes in the 
combination a/j? e and thus contributes to the 0(e) terms in (|243p . The dimensionless Poincare 
invariant J dPx (x; Xi, a^+i) depends on a single scale and, therefore, it must be a 
constant, after which (|244|) implies 

J d D x T®(x;x i -. 1 ,x i ,x i+1 ) = . (245) 

For I = 1 this relation is in agreement with the one-loop result (|24Up , 



8.5 Special conformal Ward identity 

We are now ready to investigate the special conformal Ward identity (|228jl . Inserting (I243p into 
its right-hand side, and integrating over x we obtain 

n 

K"lnW„ = ^2{2x%Xi-di-xfd?)lnW n (246) 

i=l 

_ / r (0 \ n n 

E «' nr + rW E M-m+i^ 2 ) < - 2 E Tl ^-i> **> + > 

1>1 V / i=l t=l 



where 



^(xj-i, g?j, gj+i) = a' / d D xx^T (/) (x;Xi_i,Xi,x i+ i) . (247) 
l>i J 



Next, we substitute lnW„ = Z n + F„ into (|246p . replace lnZ n by its explicit form (|176p and 
expand the right-hand side in powers of e to rewrite (|246p as follows: 

K' / F ra = lr cuS p(a)^ln * M+2 - 2 ^ T^- 1; g< , x m ) + O(e) . (248) 

i=i x i-i,i+i j=i 

Note that the quantities T iy (xj_i, X{, #i+i) are translation invariant. Indeed, a translation under 
the integral in (I247p only affects the factor x 1 ' (the functions T^ l \x; Xi-i, X{, Xj+i) are translation 
invariant), but the result vanishes as a consequence of ([245 1) . Furthermore, T u (xi-\, Xi, a?i+i) 
only depends on two neighbouring light-like vectors xf_ li and x^ i+1 , from which we can form 
only one non-vanishing Poincare invariant, x?_ li+1 . We have already argued that are 
independent of (i, and so must be T u . Taking into account the scaling dimension one of T u , we 
conclude that 

T u (x i - 1 ,Xi,x i+1 ) = ax\_ Xi + 0x v iji+1 , (249) 

where a, (3 only depend on the coupling. The symmetry of the Wilson loop W n under mirror 
exchange of the cusp points translates into symmetry of (I249p under exchange of the neighbours 
Xi-x and Xi + \ of the cusp point X{, which reduces (|249p to 

T u {xi-i,Xi,x i+l ) = a (x\_ x + Xi +1 - 2x") , (250) 
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with q = a + 0(a 2 ) according to (I240p , Substituting this relation into (I248D we find 

n 

5^T 1, (x i _ 1 ,Xi,x i+ i) = 0. (251) 
i=l 

This concludes the derivation of the special conformal Ward identity. In the limit e — > it takes 
the form (cf. |86j) 

n 1 n x 2 

Y,{1x v iXi ■ d t - x\d\) F n = -r cusp (a) ln • (252) 

i=l i=l ^i-l.i+l 

8.6 Solution and implications for F n 

Let us now examine the consequences of the conformal Ward identity (|252|) for the finite part 
of the Wilson loop W n . We find that the cases of n = 4 and n = 5 are special because here 
the Ward identity (|252|) has a unique solution up to an additive constant. The solutions are, 
respectively, 



Fi = -r cusp (a)ln 2 ^ + const , (253) 




F 5 = --r cusp (a) ^ln [ ) In [ -f^ ) + const , (254) 

i=i 



as can be easily verified by making use of the identity 



^4 = 2(^ + ^)4. (255) 



We find that, upon identification of the kinematical invariants 

x l,k+r : = (Pk + ■ ■ ■ + Pk+r-l) 2 , (256) 

the relations (|253|) and (|254|) are exactly the functional forms of the ansatz of [39J for the finite 
parts of the four- and five-point MHV amplitudes (or rather the ratio of the amplitude to the 
corresponding tree amplitude). 

The reason why the functional form of -F4 and F§ is fixed up to an additive constant is 
that there are no conformal invariants one can build from four or five points x% with light-like 
separations xf i+1 = 0. Such invariants take the form of cross-ratios 

. (257) 

x ik x ji 

It is obvious that with four or five points they cannot be constructed. This becomes possible 
starting from six points, where there are three such cross-ratios, e.g., 

22 22 22 

ry ^ ry* ry* ty* nr* ry**-' 

„ _ x 13 x 46 „ _ x 24 x 15 „ _ x 35 x 26 / 9 ro\ 
U\ — —J 5~, U 2 — —2 5~, ^3 — — o 2~ ' 

x 14 x 36 x 25 x 14 x 36 x 25 

Hence the general solution of the Ward identity at six cusp points and higher will contain an 
arbitrary function of the conformal cross-ratios. 

Nevertheless, one can prove that the functional form of the ansatz of [39] still provides a 
particular solution to the Ward identity (|252p |86j . The ansatz of [39] for the logarithm of the 
ratio of the amplitude to the tree amplitude reads 

ln^( MHV ) = Z n + Fi BDS ) + C n + O(e) , (259) 
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and 

4 BDS) = 5 Wa)^ (260) 

where Z n is the IR divergent part, F n is the finite part depending on the Mandelstam variables 
and C n is the constant term. At four points the proposed form of the finite part, written in the 
dual notation (|196|) . is 

•^ = ^ln 2 (#)+4C 2 , (261) 



while for n > 5 it is 



\ n - 2 - ^* 2 2 ^ 

= ^ y~] 9n,i , 9n,i = - ^2 ln( 2 ' ,t+T j lnf t + 1 ' t+r+1 j + D n ^ + L n j + -( 2 ■ (262) 

i=l r =2 X i,i+r+l X i,i+r+l J 1 

The functions D Ut i and L n j depend on whether n is odd or even. For n odd, n = 2m + 1, they 



arc 



m_1 xh^x 2 



r= 2 X i,i+T-+l X i~l,i+r 
1 - / X f,i+m \ , ( X 1+l,i+m+l 



L n i = —- In In 



x i,i+m+l Jj i+m,i+2m 



For n even, n = 2m, they are 

"1-2 2 2 2 2 

Dn,i = -2^Li 2 M -- 2 - 2 1- -2 J, (264) 



r —2 ' X i,i+r+l X i—l,i+r' ~ * Uy i,i+Tn^i-l,i+m—l 



1 



L„ i = - In 



We have already seen that at four points and five points the general solution to the Ward iden- 
tity coincides with (|259[) . We now show that the ansatz (|259|) is a solution of the Ward identity 
for arbitrary n. 

First we observe that the dilogarithmic contributions in (I263p and (|264p are functions of 
conformal cross-ratios of the form (|257p . They are therefore invariant under conformal trans- 
formations and we have immediately 

K^D n .i = . (265) 
For the logarithmic contributions we use the identity (I255p . When n is odd we then find 

m— 1 

K^9n,i = —2 [xf +r ^ +r+1 (lnx 2 +1 ^ +r+1 — lnx 2 i+r+l ) — xf i+1 (lnx 2 i+r — lnx 2 j +r 2 )] (266) 

r=2 

~ X i+m,i+m+lO Jil X i+l,i+m+l ~ m x i+m,i+2m) ~ ( x i+l,i+2m ~ X i+m,i+m+l ) ( m x i,i+m ~ ^ nx i+m+l,i) • 

Changing variables term by term in the sum over i one finds that only the lnx 2 i+2 terms remain 
and indeed (I252p is satisfied. The proof for n even goes exactly the same way except that one 
obtains 

- 1 

K^Ffi = [in x j,j+2 ( x j + x i+2 ~ % x i+l) + 2 ^ nX i,i+m( x i+m-l,i+m+l ~ X i-l,i+l)] (267) 
i=l 

and one has to use the fact that n = 2m to see that the ln xf i+m term vanishes under the sum. 

Thus we have seen that the BDS ansatz for the n-gluon MHV amplitudes satisfies the 
conformal Ward identity for the Wilson loop. 
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9 Hexagon Wilson loop and six-gluon MHV amplitude 



In section [7J we saw that the duality © is valid at two loops for n = 4, 5 points (or gluons). 
For n = 4, 5, the Ward identities derived in section [8] for the Wilson loops provide a possible 
explanation of the duality to all orders in the coupling constant, assuming the BDS ansatz is 
correct in these cases. This speculation receives support from strong coupling, where the result 
of [79] is in agreement with the BDS ansatz for n = 4. 

Recall that the authors of [80J found a disagreement with the BDS ansatz at strong coupling 
for the number of gluons n very large. Given this result it seems natural to ask whether the 
ansatz already breaks down at some perturbative level and for a finite number of gluons. If we 
assume the duality with Wilson loops, then n = 4 and n = 5 are completely fixed by conformal 
symmetry. In fact, one might even suspect that the observed validity of the duality relation 
for n = 4, 5 is true only because both objects have the same conformal symmetry. To put it 
differently, the proposed duality might be reduced to the weaker (but still highly non-trivial) 
statement that the MHV gluon amplitudes are governed by dual conformal symmetry. The 
first real test of the stronger form of the duality © is provided by the case n = 6. In order 
to distinguish between these possibilities we performed a two-loop calculation of the hexagonal 
Wilson loop. @ 



9.1 Finite part of the hexagon Wilson loop 

The finite part of the hexagon Wilson loop, Fg^\ does not depend on the renormalisation 
scale and it is a dimensionless function of the distances a;?-. Since the edges of Cq are light-like, 



= 0, the only nonzero distances are xj i+2 an(a x fi+3 (with i = 1, . . . , 6 and the periodicity 
condition = xi). We saw in section [8] that it has to satisfy the conformal Ward identity 
(|252|) . Specified to n = 6, its general solution is given by [86J 

F 6 (WL) = F 6 (BDS) + R 6 ( Ul ,u 2 ,u 3 ; a) . (268) 
Here, upon the identification pi = Xi + \ — Xi, 



p (BDS) _ l r / x v-^ 
r 6 ~ 7 1 cusp W 



i=l 



, f X i,i+2\-, / x i+l,i+3 
1- 2/ X i,i+3 \ 1 T . X i,i+2 X i+3,i+5 



(269) 



+ Li 2 1 , 

4 Vx i+l,i+4 7 Z K X i,i+3 X i+2,i+5 

while Rq(u\, U2, U3; a) is an arbitrary function of the three cross-ratios El 

22 22 22 

.. _ x 13 x 46 _ x 24 x 15 _ x 35 x 26 (17(X\ 

Ui — —5 5-, 112 — ~n 9-, ^3 — —9 9- • l 2 ' U J 

X 14 X 36 x 25 x 14 x 36 x 25 

These variables are invariant under conformal transformations of the coordinates xf, and there- 
fore they are annihilated by the conformal boost operator entering the left-hand side of (|252p . In 
addition, the Wilson loop W{C§) is invariant under cyclic (xj — ► Xj+i) and mirror (xj — » xg_j) 
permutations of the cusp points [87J. This implies that R(ui,U2,U3) is a totally symmetric 
function of three variables. 



42 This section is based on the publications [90] and [91] . 

43 The last term in (|269[l is a function of cross-ratios only, but we keep it in Fg BDS ', because it is part of the 
BDS conjecture. 
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The one-loop Wilson loop calculation of [H] has shown that the 'remainder' function Rq is 
just a constant at one loop, which is a confirmation of the Wilson loop/scattering amplitude 
duality going beyond the scope of conformal symmetry. However, one might suspect just a low 
loop-order 'accident'. The point is that the function Rq must satisfy a further, rather powerful 
constraint, coming from the collinear limit of gluon amplitudes |52} 156] . It could be that due 
to the limited choice of loop integrals at this low perturbative level, the function (|269|) made of 
them is the only one satisfying both the conformal Ward identity (|252l) and the collinear limit. 
If so, at some higher perturbative level new functions with these properties might appear which 
could spoil the BDS ansatz and/or the Wilson loop/scattering amplitude duality. 

Combining together (|195p and (I268p . we conclude that were the BDS conjecture and the du- 
ality relation (I195p correct for n = 6, we would expect that R(ui, 112,113) = const. The explicit 
two-loop calculation we report on here shows that this is not true. 



For the sake of simplicity we performed the calculation of W(Cq) in the Feynman gauge. In 
addition, we made use of the non-Abelian exponentiation property of Wilson loops |176Lll94lll95j 
to reduce the number of relevant Feynman diagrams. In application to R(ui, U2, U3) this property 
can be formulated as follows (the same property also holds for Fg WL ^) 



r = £- 2 c fR w + ^ 2 c f nr^ + o^ 



(271) 



where Cp = (N 2 — 1)/(2N) is the Casimir in the fundamental representation of SU(N). The 
functions R^ and R^ do not involve the colour factors and only depend on the distances 
between the cusp points on Cq. At one loop, R^>(ui, 112, U3) is in fact a constant 



As explained in [86] , the relation (|27ip implies that in order to determine the function .Fg WI ^ 
at two loops (and hence f^) it is sufficient to calculate the contribution to W{Cq) from two-loop 
diagrams containing the 'maximally non-Abelian' colour factor C F N only. All relevant two- loop 
graphs are shown in Fig. The finite part of W(Cq) receives contributions from all Feynman 





Figure 29: The maximally non-Abelian Feynman diagrams of different topology contributing to Fq WL \ 
The double lines depict the integration contour C§, the dashed lines depict the gluon propagator and the 
blob stands for the one-loop polarisation operator. 
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diagrams shown in Fig. 1291 We derived parameter integral representations for them, which are 
denoted by A^\ so that the two-loop contribution to Fg WL ^ = Yl p >i aP -^6^p ^ s gi ven by 

4T = 2 E A o a) - ( 272 ) 

a=a,...,u 

The integrals Aq are difficult to compute analytically and so we evaluated them numerically, 
as described in the next section. 



9.2 Numerical evaluation 

For a given set of kinematical invariants, 

7^_r2 2 2 2 2 2 2 2 2 \ f97Q^ 
n — i x 13' x 14> ^15' x 24> x 25' x 35' x 26' x 36> X 46J > l^'^J 

defined as the distances between the vertices of the hexagon C6, it is straightforward to evaluate 
the sum on the right-hand side of (|272p numerically. We should take into account, however, that 
in D = 4 dimensions there exists a relationship between the distances (|273p . such that only eight 
of them are independent. This relation reflects the fact that the five four-dimensional vectors 
Pi = x i+i ~ x i ( w hh i = 1, . . . , 5) are linearly dependent, and hence their Gram determinant 
vanishes 

G[K) = det \\(pi ■ pj)\\ = , (i,j = l,...,5). (274) 

Using pi = Xi + \ — xi, the entries of the matrix can be written as linear combinations of the 
distances (|273p . 



We would like to note that the sum on the right-hand side of (|272p defines a function of xf- 
even for configurations K which do not satisfy (I274p . This function can be viewed as a particular 
continuation of the finite part of the hexagon Wilson loop off the hypersurface defined by (|274p • 

Let us now consider the relation between the hexagon Wilson loop and the BDS ansatz, 
Eq. (|268p . To two-loop accuracy it takes the form 

Rw(ui,u 2 ,u 3 ) = F$ L) - Fg DS) , (275) 

where i?w and F^ DS ^ denote the two-loop contributions to the remainder function (|27ip and 

to the BDS ansatz (|269p . respectively. We recall that the functions and Fj^^ satisfy 

the Ward identity (I228p . and hence R\y is a function of the three conformal cross-ratios (|270l) 
only. 

The simplest way to check relation (I275p is to evaluate the difference Fq^ L ^ — F^ DS ^ for 
different kinematical configurations K and K', related to each other by a conformal transfor- 
mation of the coordinates (with i = 1,...,6). Since F^ L ^ and -Fg.^ ^ are dimensionless 
functions of the distances xfj, they are automatically invariant under translations, Lorentz rota- 
tions and dilatations of the coordinates xf. The only non-trivial transformations are the special 
conformal transformations (boosts), which are combinations of an inversion, a translation and 
another inversion. 



Let us start with a kinematical configuration K = K(xfj), Eq. (|273p . and perform an inver- 
sion of the coordinates, xf — ► xf/xf, to define the new configuration 

K< = K{x%l(x\x 1 j j). (276) 
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Since the variables u (I270P are invariant under such transformations, u a [if] = u a [K'], the dif- 
ference F^ L ^ — i ? g.2 DS ' ) should also be invariant. 

As an example, let us consider six light-like four-dimensional vectors pP, 

Pi = (1,1,0,0), p 2 = (-l,p,p,Q), p 3 = (l,-p,p,0), (277) 

Pa = (-1,-1,0,0) , p 5 = (l,-p,-p,0) , p e = (-l,p, -p,0) , 

with p = l/y/2 and X^jPf = 0- These vectors define the external momenta of the gluons in the 
six-gluon amplitude. Applying the duality relation pi = x,+i — Xj, we evaluate the corresponding 
distances (12731) . 

: x 2 4 = xf 5 = X24 = X25 = —2 , x§ 6 = —2 — 2y/2 , 

x 13 = x 35 = x 26 = x 46 = — 2 — V% , (278) 
and the conformal cross-ratios (|270|) . 

«i = «3 = \ + 3^/2, «2 = 1. (279) 

By construction, this kinematical configuration satisfies the Gram determinant constraint (|274p . 
To define the conformal transformations (|276|) . we choose an arbitrary reference four- vector 
x^ = (x^, x}, x 2 , xf) and reconstruct the remaining x— vectors according to xf +1 = xf +pf (due 
to translation invariance, the Wilson loop does not depend on the choice of x^). Then, relation 
(|275p implies that the function Ryj evaluated for the kinematical configuration (I276|) should be 
the same as for the original configuration K. 



As an example, we choose x^ = (1, 1, 1, 1) and apply the conformal transformation (|276|) to 
defined in (|2?8j) to obtain the new kinematical configuration 

K {h) . T 2 _ 2 _ 2 _ 2 _ _ 1 _ 1 Hy 2 _ _-, _ 3 /n 

i\ . X 14 — X 15 — X 2 4 — X 2 5 — 2 4 V Z ' X 36 — 1 4 VZ ' 

r 2 -_3_,/o r 2 __5_7./o r 2 

x 13 — 2 VZ > x 35 — 2 4 VZ ' ^26— 4 8 VZ ' 

»46 = -f-?V^. (280) 

The results of our numerical tests are summarised in Table [TJ They clearly show that i^g^j and 

i ? g.2 DS ' 1 vary under conformal transformations, whereas their difference i?w = F^ L ^ — ^e^^ 
stays invariant. 



We recall that in four dimensions the kinematical invariants (I273D have to verify the Gram 
determinant constraint (|274j) . This relation G[K] = is invariant under the conformal trans- 
formations (I276p . simply because the conformal boosts map six light-like vectors pf into an- 
other set of light-like vectors. There exist, however, certain kinematical configurations K' for 
which itifX] = iti[-K"'] but G[-fC'] 7^ 0. Since the difference function only depends on the 
u— variables, its value should be insensitive to the Gram determinant condition 0. For example, 
consider the following kinematical configuration 

: x 2 4 = x 2 5 = — 1 , X24 = X25 = —2 , = —2 — 2\/2 , 

x 2 3 = -l-l/v/2, x| 5 = x 2 6 = x 2 g = -2-^2. (281) 

The corresponding conformal cross-ratios (|270p are given by (379)1 . but G[K^] ^ 0. We verified 
numerically that # w [# (a) ] = Rw[K {b) ] = Rw[K (c) ] with accuracy < 10~ 5 (see Table!]). This 
observation allows us to study the function i?\y( u i, u 2, U3) without any reference to the Gram 
determinant condition. 

44 We recall that the functions entering (|275[l can be defined for configurations K' satisfying G[K'\ 7^ 0. 



71 



Kinematical point 


-^6;2 


r 6;2 






-5.014825 


-14.294864 


9.280039 




-6.414907 


-15.694947 


9.280040 




-5.714868 


-14.994906 


9.280038 



Table 1: Two-loop contributions to the hexagon Wilson loop, K. 2 , to the BDS ansatz, 

i*g^ D , and to their difference, -R\y> evaluated for three kinematical configurations (12780 . (|280p 
and ([28 ip corresponding to the same values of iti, U2 and it3, Eq. (|279j) . 



9.3 Collinear behaviour 

We recall that for the six-gluon amplitude M.§ depending on light-like momenta, Ym=X Pi = ® 
and pf = 0, the collinear limit amounts to letting, e.g. p§ and pg be nearly collinear (see e.g. 
|52j for more details), so that (p$ + Pe) 2 ~^ and 

p£-*i**, p£-(l-z)P", (282) 

with P 2 = and < z < 1 being the momentum fraction. Using the identification pf = 
x^ +1 — xf, we translate these relations into properties of the corresponding Wilson loop W(C%). 
We find that the cusp at point 6 is 'flattened' in the collinear limit and the contour Cq reduces 
to one with five cusps. In terms of the distances x 2 j, the collinear limit amounts to 

2 2 , /-i \ 2 

x 36 ~^ zx 13 "+~ \ l ~ z ) x 35 i 

x 2 6 - (1 - z)x 2 5 , (283) 

while the other distances xf 3 , X24, X25, £35 remain unchanged. For the conformal cross-ratios the 
relation (|283|) implies 

ui — > n , U2 — > , us — > 1 — u , (284) 

with ti = zx\§/ (zx\§ + (1 — ^)x| 5 ) fixed. As was already mentioned, the relation ([268|) is con- 
sistent with the collinear limit of the six-gluon amplitude provided that, in the limit ([2840 . the 
function R(ui, U2, 1*3) approaches a finite value independent of the kinematical invariants. The 
same property can be expressed as follows (we recall that the function R(ui,U2,us) is totally 
symmetric) 

R{0,u,l-u) = c, (285) 

with c being a constant. Using our two-loop results for the finite part Fq, we performed thorough 
numerical tests of the relation (|285|) for different kinematical configurations of the contour Cq. 

We found that, in agreement with (|285p . the limiting value of the function i?( 2 )(7, u, 1 — u) 
as 7 — * does not depend on u. Since the duality relation (|198p is not sensitive to the value of 
this constant, it is convenient to subtract it from R^(pf,u, 1 — u) and introduce the function 

, u, 1 — u) = c — i? (2) (7,-u, 1 -u) (286) 

which satisfies .R (2) (0, u, 1 - u) = 0. To summarise our findings, in Fig. [30] we plot the function 
w- 2 ^ (7, u,l — u) against 7 for different choices of the parameter < u < 1 and in Fig. [31] the 
same function against u for different choices of the parameter 7. The important region for the 
collinear limit is where 7 is close to zero. We also give numerical values for a range of values of 7 



x 15 ~> 1 



' 16 
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0.2- 



0.05 




Figure 30: The 7— dependence of the function w 2 ^(^/,u, 1 — u), Eq. (|286p . for different values of the 
parameter u — 0.5 (lower curve), u = 0.3 (middle curve) and u — 0.1 (upper curve). 



0.25- 



0.15 




0.05- 



Figure 31: The u— dependence of the function R^(-f,u, 1 — u), Eq 
parameter 7 = 0.001 (lower curve), 7 = 0.01 (middle curve) and 7 = 0.1 (upper curve) 



for different values of the 



such that one can see how the function R( 2 \ui, 112, W3) 

U\ = 7, U% = U, U3 = 1 — U. 



varies in the particular parametrisation 



So, in conclusion, our explicit n = 6 Wilson loop calculation [90] has shown that at two loops 
there exists a non-trivial 'remainder' function Rq satisfying the conformal Ward identities and 
the collinear limit condition. The crucial test then was to compare the results of our Wilson loop 
calculation with a parallel two-loop six-gluon amplitude calculation, in order to check whether 
the proposed duality between Wilson loops and gluon amplitudes continues to hold at this level 
(which, if true, automatically implies the breakdown of the BDS ansatz). The results of the 
six-gluon calculation have become available very recently [92], and we present the comparison 
in section 19. 
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9.4 The hexagon Wilson loop versus the six-gluon MHV amplitude 



The recently completed two-loop six-gluon amplitude calculation [92] found that indeed the BDS 
ansatz fails at that level, i.e. 

Ra = F^ UV) - if 2 DS > . (287) 

with a nontrivial 'remainder function'. We now compare numerically the 'remainder' function 
Ra with the analogous quantity defined for the hexagonal Wilson loop in (I275|) . In order to 
test the duality relation (|195j) . we have to show that for general kinematical configurations K 
in (|273l) . we have 

R W [K]-R A [K] =c w - (288) 
Since i?w is a function of u±, U2, U3 only, this relation would imply that so is Ra- 

To get rid of the constant cw, which has lower numerical precision than the evaluations 
for generic kinematics, we subtract from (|288p the same relation evaluated for some reference 
kinematical configuration K^°\ 

R A [K] - R A [K {0) ] = R W [K] - R W [K^] . (289) 

The numerical tests of this relation for different kinematical configurations are summarised in 
Table [2J 



Kinematical point 


(ui,u 2 ,u 3 ) 


Rw — Ryy 


-Ra — R^ 




(1/4,1/4,1/4) 


< 10~ 5 


-0.018 ±0.023 


K m 


(0.547253, 0.203822, 0.88127) 


-2.75533 


-2.753 ±0.015 


K m 


(28/17, 16/5,112/85) 


-4.74460 


-4.7445 ± 0.0075 


K w 


(1/9,1/9,1/9) 


4.09138 


4.12 ±0.10 




(4/81,4/81,4/81) 


9.72553 


10.00 ±0.50 



Table 2: Comparison of the deviation from the BDS ansatz of the Wilson loop, Rw, an d of 
the six-gluon amplitude, Ra, evaluated for the kinematical configurations (|290p . Here, R™ = 
# w (l/4, 1/4, 1/4) = 13.26530 and R$ = fi A (l/4, 1/4, 1/4) = 1.0937 ± 0.0057 denote the same 
quantities evaluated at the reference kinematical point . The numerical results for Ra and 
are taken from Ref. |92j . 
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K K 1 : 


2 

x i,j- 


^2 


— 1 J x i,i+3 — 


-2; 












r 2 
x 13 


= 


-0.7236200, 


r 2 
x 24 


= -0.9213500. 


r 2 
x 35 


= -0.2723200 , 


r 2 - 
x 46 — 


-0.3582300 




r 2 


= 


-0.4235500, 


r 2 
x 26 


= -0.3218573. 


r 2 
x 14 


= -2.1486192, 


r 2 - 
x 25 — 


-0.7264904 




2 
x 36 


= 


-0.4825841 ; 
















r 2 
x 13 


= 


-0.3223100, 


r 2 
x 24 


= -0.2323220. 


r 2 
x 35 


= -0.5238300 , 


r 2 - 
x 46 — 


-0.8237640 




r 2 


= 


-0.5323200 , 


2 

x 26 


= -0.9237600 . 


r 2 
x 14 


= -0.7322000 , 


r 2 - 
x 25 — 


-0.8286700 




x 36 




-0.6626116; 
















x a,i- 


h'2 


= — 1 , x 2 4 




1/2, xL = - 


•5/8, 


xL = -17/14; 






is (4) 


2 


1-2 


--1 x 2 - 

^ ! X J,J- 


f3 = 


-3; 










tf(5). 


x 2 - 
x l,«- 


f2 


— 1 ) x i,i+3 ~ 


-9/2. 











Among the six configurations in equation (|290|) only the first four verify the Gram determinant 
condition 12741 Also, the configurations and are related to each other by a conformal 
transformation. This explains why the first entry in the third column of Table [2] is almost zero, 
and also reflects the high precision of the numerical evaluation for the Wilson loop. Comparing 
the numerical values for the six-gluon amplitude and the hexagon Wilson loop, we observe that 
their finite parts coincide within the error bars. Therefore, we conclude that the duality relation 
(|195jl is satisfied to two loops at least. We consider this very strong evidence that the duality 
should hold to all orders in the coupling, although the 'remainder' function R n is likely to receive 
corrections at each loop order. 



10 Conclusions and outlook 

In this report, we studied a new duality between Wilson loops and gluon amplitudes. The latter 
was motivated mainly by work of Alday and Maldacena, which showed that there is a relation 
between gluon amplitudes and Wilson loops at strong coupling, using the AdS/CFT correspon- 
dence. Unexpectedly, the duality seems to hold also at weak coupling. Previous studies of the 
duality were restricted to one loop. In this thesis we computed the expectation value of the 
rectangular and pentagonal Wilson loops at two loops. A comparison with the corresponding 
MHV gluon amplitudes shows that the duality holds. 

As was emphasised in this report, broken (dual) conformal symmetry plays an important role 
on both sides of the duality. We showed that the integrals contributing to the loop corrections 
of the four-gluon amplitude up to four loops are all 'pseudo-conformaP in the sense explained in 
section 13.21 The origin of the dual conformal symmetry for gluon amplitudes remains puzzling. 
For the Wilson loops, however, it is completely understood and controlled by an all-orders Ward 
identity, which we derived in this thesis. 

The case n = 6 provides a test of the duality relation beyond (dual) conformal symmetry. 
What makes the hexagonal Wilson loops and the six-gluon MHV amplitude agree? It seems 
natural to conjecture that (dual) conformal symmetry is part of a larger symmetry group, which 
is yet to be discovered. Such a speculation is further supported by the fact that both Wilson 
loops and gluon amplitudes are intimately related to the cusp anomalous relation, which in turn 
is governed by integrability. 
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The results obtained give strong support for the duality between Wilson loops and gluon 
scattering amplitudes in N = 4 SYM. In our opinion, given the evidence presented here, it is 
likely that the duality is true to all orders in the coupling constant, and for an arbitrary number 
of points/gluons. This duality gives rise to several intriguing possibilities. 

The calculation of the pole part of the n-cusp Wilson loop in [87] was carried out analytically 
for arbitrary n at two loops, in agreement with the general formula (I176p . As explained in [87] . 
one can define 'auxiliary' diagrams to subtract the pole terms, thereby providing a definition 
of the finite part of the two-loop n-cusp Wilson loop in terms of a set of finite parametric inte- 
grals. We studied the latter numerically in the n = 5, 6 discussed in this report, and a 
generalisation to n > 6 should not be difficult. 

It would be very interesting if the arguments of Alday and Maldacena that related gluon 
scattering amplitudes and Wilson loops at strong coupling could be extended to arbitrary val- 
ues of the coupling constant. This would constitute a proof of the duality via the AdS/CFT 
correspondence, or conversely, such a derivation would be an impressive test of the AdS/CFT 
correspondence. 

So far, the duality is limited to maximally helicity violating (MHV) amplitudes. It is natural 
to wonder whether or how the duality extends to amplitudes with different helicity configura- 
tions, as for example NMHV. An obvious question is then how to match the helicity structure of 
the NMHV amplitudes to the Wilson loops, which are helicity independent. One may speculate 
that there could exist suitable modifications of the Wilson loops such that they match those 
amplitudes. For example, one could consider Wilson loops with local operator insertions. 

A related question is what the dual conformal properties of NMHV amplitudes are. This 
will be the subject of a forthcoming paper together with J. Drummond, G. Korchemsky and E. 
Soktachev. 

Another interesting possibility, in our opinion, would be to consider a different point of the 
moduli space of N = 4 SYM, such that some particles become massive. One can ask the question 
whether scattering amplitudes of massive on-shell gluons are still dual to the corresponding 
Wilson loops, whose sides are no longer light-like (because of the mass). On the other hand, for 
finite masses the conformal symmetry that was important in our analysis would be lost for such 
a configuration. 

11 Acknowledgements 

First of all, it is a pleasure to thank my advisor, Emery Sokatchev, for initiating me to this 
interesting area of research and for sharing his knowledge. Special thanks are also due to Andrei 
Belitsky, James Drummond and Gregory Korchemsky for numerous discussions and encour- 
agement. I have also profitted from discussions with Costas Bachas, Benjamin Basso, Niklas 
Beisert, Massimo Bianchi, Fawzi Boudjema, Lance Dixon, Vladimir Dobrev, Burkhard Eden, 
Johanna Erdmenger, Giovanni Feverati, Laurent Gallot, Francois Gieres, Jean-Phillipe Guillet, 
Paul Heslop, Candide Jarczak, Dieter Muller, Eric Pilon, Jan Plefka, Eric Ragoucy, Paul Sorba, 
Yassen Stanev, Kelly Stelle, Raymond Stora, Frank Thuiller, Ivan Todorov and Arkady Tseytlin. 
I am grateful to Nans Baro, Livia Ferro, Charlotte Grosse-Wiesmann and Konstantin Wiegandt 
for pointing out typos in the draft to me. 

I would like to thank the Laboratoire d'Annecy-le-Vieux de Physique Theorique (LAPTH) 
and its director, Patrick Aurenche, for providing a stimulating working environment during the 



76 



last three years. The secretaries, Nicole Berger, Veronique Jonnery, Dominique Turc and Vir- 
ginie Malaval have helped considerably to simplify the administrative part of this thesis. It 
has been great to share the experience of doing a PhD with the fellow students Nans Baro, 
Samuel Belliard, Christophe Bernicot, Candide Jarczak, Ninh Le Due, Gregory Sanguinetti, 
Wessel Valkenburg. They and the other members of LAPTH made the three years much more 
enjoyable and worthwile. 

This research was supported in part by the French Agence Nationale de la Recherche under 
grant ANR-06-BLAN-0142. Partial financial support from the Universite de Savoie through an 
Explora'Doc grant and travel money through an APS grant (2005) is gratefully acknowledged. 
It is a pleasure to acknowledge the warm hospitality extended to me by the Theory Group of 
the Dipartimento di Fisica, Universita di Roma "Tor Vergata", where part of this work was 
done (from April 2007 to September 2007). Special thanks are due to Yassen Stanev for his 
hospitality and numerous enlightening discussions. 

I would like to thank the theoretical physics sections of the Universita di Roma "Tor Ver- 
gata", the Universita Torino, the Univerisita Milano Bicocca, the ETH Zurich, and Imperial 
College London, where part of this work was presented, for their seminar invitations. In par- 
ticular I would like to thank the participants of these seminars for their comments and for 
discussions, which certainly improved the way the material is presented in this thesis. 

Finally, this thesis would not have been written had it not been for the university education 
I received. At the Universitat Augsburg I would like to thank in particular Reinhard Schertz 
and Gert-Ludwig Ingold for their excellent lectures that aroused my interest in mathematics 
and theoretical physics, and similarly Daniel Decamp at the Univerite de Savoie for initiating 
me to the world of elementary particle physics. Finally, I would like to thank Luc Frappat, 
Frank Thuiller and Jean-Claude Le Guillou for their help in organising my Erasmus year at 
the Universite de Savoie, and Fraingois Delduc for the organisation of the theoretical physics 
master's programme at the Ecole Normale Superieure de Lyon. 

A Alternative proof of <5>W = ^r(3) 

using the Mellin— Barnes rep- 
resentation 

In this section we show how the above identity between the off shell triple box and tennis court 
integrals derived in section 12.61 can also be obtained by means of the method of Mellin-Barnes 
(MB) representation. 

A.l Introduction to the Mellin— Barnes technique 

This method is one of the most powerful methods of evaluating individual Feynman integrals. 
It is especially successful for evaluating four-point Feynman integrals. For massless off-shell 
four-point integrals, first results were obtained by means of MB representation in [143] 1144] . 
In the context of dimensional regularisation, with the space-time dimension d = 4 — 2e as a 
regularisation parameter, two alternative strategies for resolving the structure of singularities 
in e were suggested in [203] 1204] where first results on evaluating four-point on-shell massless 
Feynman integrals were obtained. Then these strategies were successfully applied to evaluate 
massless on-shell double [2031 EH EH EH EH EH [209] and triple [2101 EEH EH] boxes, with 
results written in terms of harmonic polylogarithms [212] . double boxes with one leg off shell 
[23112131 and massive on-shell double boxes [251 EES EEZ1 EE] (see also chapter 5 of [219]). It 
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A Im(z) 



-A /3 



X- 



-X- 



Re(z) 



Figure 32: Position of the integration contour (dashed line) and poles of the integrand (crosses) 
in formula (|29ip . for real A. 



is based on the MB representation 

1 1 



(X + Y) x r(A)2vri 



/3— joo 



X A + 2 



T(X + z)T(-z)dz, 



(291) 



applied to replace a sum of terms raised to some power by their product to some powers. The 
integration in (|29ip goes along a straight line parallel to the imaginary axis, with — Re(A) < 
[3 < 0, so that the poles of T(A + z) lie on its left and those of T(—z) on its right, as shown in 
Fig. [32j Formula (|291l) is easy to understand. If Y < X, we can close the integration contour 
on the right. By the residue theorem, we will pick up contributions from the poles of T(—z), 
leading to 



1 1 rP+ioo y z 



T(A) 2m 



P—ioo 



T(A + z)T(-z) dz = X- x J^(-l) 

n=0 



r(A + n) fry 
' r(A) \x) ■ 



(292) 



The right-hand side of (|292p is indeed the series expansion of (X + Y) X for Y < X. However 
equation (|291|) is also valid for X < Y, which can be seen by closing the integration contour on 
the left. 



The first step of the method is the derivation of an appropriate MB representation. It 
is very desirable to do this for general powers of the propagators (indices) and irreducible 
numerators. On the one hand, this provides crucial checks of a given MB representation using 
simple partial cases. (For example, one can shrink either horizontal or vertical lines to points, 
i.e. set the corresponding indices to zero, and obtain simple diagrams quite often expressed in 
terms of gamma functions.) On the other hand, such a general derivation provides unambiguous 
prescriptions for choosing integration contours (see details in |219j). 

A. 2 Example: One-loop integral 

As an example, let us come back to the one-loop integral $W. Using the relation to the three- 
point integral (|73p it is straightforward to obtain the following parameter integral representation: 

v) - (2m)-* C d MMM^ -01-02- &») (m) 
or, after a change of variables, 

,i ,i A 



$W(u,v) = (2iriy z I dx I dX-. . 

o jo X{xu + xv) + Xxx 



(294) 



where x = 1 — x and A = 1 — A. Then, one uses the Mellin Barnes formula (|29ip twice in order 
to factorise the denominator in (|293l) . i.e. 



a + b- 



/zoo ficyo 
d Zl \ dz 2 T{-z l )T{-z 2 )T{l + z l + z 2 )a^b^c i -- 1 '^- z ^ 
-ioo J —ioo 



(295) 
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As usual, the integration contours are straight vertical lines, and they are chosen between the 
poles of the T functions, i.e. 



Re(zi)<0, Re(z 2 )<0, Re(*i) + Re(z 2 ) > -1 , (296) 

so one could choose Re(zi) = —0.2 and He(z 2 ) = —0.4 for example. The integration over the 
parametric integrals is now straightforward, and we obtain the MB representation [143] (see also 



dzt / dz 2 u zl v Z2 T 2 (-z 1 )T 2 (-z 2 )T 2 (l + z 1 + z 2 ). 

-ioo J— ioo 



(297) 



For < u, v < 1 one can close the integration contours in (|297p on the right and use the residue 
theorem to compute the integrals. This yields the result as a series expansion in u and v. 

The representation (|297|) is also useful for studying the limiting behaviour of As an 

example, consider the OPE limit x± 2 — ► 0, i.e. u — ► 0, cf. ([771) . Since Re(zi) < the limit 
cannot be taken naively in the integrand of (|297| . Therefore, we deform the contour such that 
< Re(zi) < 1. In doing so, we cross a pole at u = and hence pick up a contribution from 
the residue of u zl F 2 (— zi)T 2 (l + z\ + z 2 ). Neglecting terms that vanish as u — > 0, we obtain 

/ioo 
v Z2 T 2 (-z 2 )T 2 {\ + z 2 )dz 2 
-ioo 

/ioo 
v Z2 T 2 (-z 2 )T 2 (l + z 2 ) [tf (1 + z 2 ) + 1E ] dz 2 + O(u) . (298) 
-ioo 

In the second line je is the Euler constant and for integer v the expression ^(1 + v) + 7^ is 
equal to an harmonic sum SV, = Y2n=l ^/ n - Computing the integrals in (|298j) by closing the 
contour on the right we pick up residues from v = 0, 1, . . . 

R eSz2=n [v Z2 r 2 (-z 2 )T 2 (l + z 2 )] = v n lnv (299) 
R eSz2=n [v Z2 T 2 (-z 2 )T 2 (l + z 2 )(^(l + z 2 )+ lE )] = 2v n + n) + 5 n lnv] . (300) 

The resulting series can be summed. The first is just a geometrical series, and the second one 
can be done as well (e.g. using Mathematica) . We obtain 



7T 2 



1-V 

or equivalently (still for < v < 1), 



In (u) ln(u) - 2 ln(i>) ln(l - u) H 2Li 2 (u) 

3 



+ O(u) , (301) 



$ w (u,v) = ^— - [In («) ln(v) + 2Li 2 (1 - «)] + O(u) , (302) 

in agreement with (|77p . Although this expression was obtained for < v < 1, it can be analyt- 
ically continued to all v > 0. 

Many other examples of the application of the Mellin-Barnes technique to the evaluation of 
Feynman integrals can be found in the book [219] . 

A. 3 Alternative proof of the magic identity at three loops 

Let us now use the MB technique to find an alternative proof of the 'turning identity' $( 3 ) = \E'( 3 ) . 
We consider the off shell triple box and tennis court labelled as shown in Figs. [331 and I34[ with 
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Figure 33: Labelled triple box. 
Pi 1 T ^3 




Figure 34: Labelled tennis court. 



general powers of the propagators and one irreducible numerator in the tennis court chosen as 
[{h + h) 2 ]~ ai1 ■, where are the momenta flowing through lines 1 and 3 in the same direction. 

Experience shows that a minimal number of MB integrations for planar diagrams is achieved 
if one introduces MB integrations loop by loop, i.e. one derives a MB representation for a one- 
loop subintegral, inserts it into a higher two-loop integral, etc. This straightforward strategy 
provides the following 15-fold MB representations for the dimensionally regularised off-shell 
triple box and tennis court with general indices: 



Ti(ai,.. . ,a 10 ;s,t,pl,pl,pl,pl;e^ - 



{ilt d l 2 f (-l) a (- S ) 6 - a - 



3c 



nj=2,4,5,6,7,9 r ( a i) r ( 4 - 2e - a 4,5,6,7) 



j r+ioo 15 



dz 



(-pj) Z12 (-pi) 213 (-p|) 24 - 9 ' 14 (-p!)^ 5 - 10 - 15 (-ty 

( g^Z4,5,9,10,ll, 12, 13, 14, IB 



^ r(a 9 + ^ii,i2,i3)r(a 7 + zi,2,3)r(2 - e - a 5 , 6 ,7 - zi,2,4)r(2 - e - 04,5,7 ~ 21,3,5) 
T(ai - z 2 )r(a 3 - z 3 )T(4 - 2e - 01,2,3 + 21,2,3) 
r(a 5 + zi,4,5)r(a4,5,6,7 + e - 2 + ^1,2,3,4,5)^(2:11,14,15 ~ gg) 
T(a 8 - z 7 )r(ai - 2 8 )r(4 - 2e - a 8j9i i + 2 6 , 7 , 8 ) 
xT(2 - e - a 8j9 + z 6 , 7 - Zn,i2,i4)r(2 - a 2 ,3 - e + 21,3 - 2 6 ,8,io) 
xr(a 8) 9,io + e — 2 + ^n, 12,13,14,15 — ^6,7,8)F(2 — e — ag,io + ^6,8 — 211,13,15) 
xr(a 2 + z 6)7i8 )r(2 - e - ai, 2 + z 1>2 - 26,7,9) 

xr(z 6 , 9 ,io - 2i)r(ai, 2 ,3 + e - 2 - Zi, 2 ,3 + 2 6 ,7,8,9,10) \{ T(-Zj) ; 

7=2,3,4,5,7,.. .,15 



(303) 
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T 2 (a 1 , . . ., a n ;s,t,p 1 ,p 2 , p 3 ,p A ;e 



(i7r d / 2 ) 3 (-l) a (- S ) 6 - a - 3t 
n j= 2,4,5,6,7,9 r ( a i) r ( 4 - 2e - 04,5,6,7) 



— r'°°TTd (-pj) 212 (-pi) zl3 (-pi) 25 - 10 - 14 (-pp' 15+Z8 (-*) 211 rr rf x 

15 / . 11 i (_ tAZ5, 8, 10, 11, 12, 13, 14, 15 11 ^ 3' 



^ r(a 9 + -Zn,i2,i3)r(a7 + ^i,2,3)r(2 - a 5 , 6 ,7 - e - 2i, 2 , 4 )r(2 - a 4 , 5 ,7 - e - 21,3,5) 
r(oi - z 2 )T(a 3 - z 3 )r(4 - 2e - 01,2,3 + Zl,2,3)r(ai - 27) 

^ r(a 5 + ^l,4,5)r(Q4,5,6,7 + £ ~ 2 + Zl,2,3,4,5)r(2 - a 2 , 3 ~ 6 + Zl,3 ~ 2 6 ,8,lo) 

T(8 - 4e - a - 2 5 , 6 , 8 ,i )r(a 8 - 24,9)1X01,2,3,4,5,6,7,11 + 2e - 4 + 24,5,6,7,8,9,10) 
xT(6 — a + aio — 3e — 25,6,7,8,10,11, 12,i4)r(a + 3e — 6 + 25,6,8,10,11,12,13,14,15) 
xT(a 2 + 2 6 , 7 , 8 )r(2 - e - ai, 2 + 21,2 - 2 6 ,7 i9 )r(6 - 3e - a + a 8 - 24,5,6,3,9,10,11,13,15) 
xr(z 6 , 9i i - 2i)r(ai, 2 ,3 + e - 2 - 21,2,3 + ^6,7,8,9,10) 

xr(ai, 2,3,4,5,6, 7,11 + 2e — 4 + 24,5,6,7,8,9,10,11,14,15 ) • (304) 

Here 04,5,6,7 = 04 + 05 + a 6 + 07, a = ^ aj, 211,12,13 = 211 + 212 + 213, etc. Moreover, the letters 
s and t denote the usual Mandelstam variables s = (pi + P2) 2 and t = (p\ + P3) 2 . 

These representations are written for the Feynman integrals in Minkowski space. (This is 
rather convenient, in particular this allows one to put some of the legs on-shell.) The corre- 
sponding Euclidean versions are obtained by the replacements — s — > s, — t — ► t, —p\ — > p\ , . . . 
and by omitting the prefactors (— l) a and i 3 . 

To calculate the triple box we need, i.e. if = T\ (1, . . . , 1) at d = 4, we simply set all the 
indices a, to one. We cannot immediately set e = because there is T(— 2e) in the denominator. 
The value of the integral is, of course, non-zero, so that some poles in e arise due to the 
integration. To resolve the structure of poles one can apply Czakon's code [220J, which provides 
the following value of the integral in the limit e — > after relabelling the variables by 210 — > 

22, 214 -> 2 3 , 215 2 4 , 2n 2 5 , 212 ~> 2 6 : 



,(o)_ K) 3 /' +too A^.(-p?) 26 (-pi)~ 1 ^ 6 (-P3)" 1 - Z5 ' 6 (-p!) 26 (-i) 



r (oj _ i"' ; / tt d 



zr, 



T(l + 2 3 , 4 )r(l + 2i - 2 3 ,4,5)r(2 2 ,3,4,5 - 2l)r(2 4 - 2 6 ) y-r 

1 + 2 4 - 2 6 )r(i + 22,4 - 2 6 )r(2 + 2 lj6 - 2 2 , 4 )r(2 + 23,5,6) 1/ Zj 
xr(2 2 , 4 - 2 6 ) 2 r(i + 21,6 - 2 2 , 4 ) 2 r(i + 2 5)6 )r(i + z 3i5>6 ) . (305) 



x r 



To calculate the tennis court we need, i.e. if = T 2 (l, . . . , 1, —1) at d = 4, we proceed like in the 
previous case. Czakon's code provides the following integral (after relabelling 210 — > 2 2 ,2i 4 — > 

23, 215 -> 2 4 , 2n -> 2 5 , 212 2 6 ): 



r (0) = f +l °° A , (-p 2 )^(-p 2 )- 1 -^-^(- p 2 ) -l- 2 5-.6 ( _ p 2 r6( _ tr 

2 (2vri)6 7_ loo 1 = 1 Zj (-s)l-** 



T(l + 2 3 , 4 )r(l + 2i - 2 3 ,4, 5 )r(2 2 ,3,4,5 ~ 2l)r(2 4 - 2 6 ) -pr 

x r(i + 2 4 - 2 6 )r(i + 21 - 2 2 , 3 ,5, 6 )r(2 + 2 3 ,5, 6 )r(2 + 22,3,5,5) ^ Zi 

xr(2i - 2 2 , 3 ,5, 6 ) 2 r(i + 2 5>6 )r(i + 2 3 ,5, 6 )r(i + 2 2 , 3 ,5, 6 ) 2 • (306) 

Now the simple change of variables 2 2 — > — 2 2 + 21 — 23 — 24 — 25 in (I306p leads to an expression 
identical to (|3U5|) up to a factor of u and we obtain the identity 2~f = slf , which corresponds 
to the identity $^ 3 ^ = of section [2.6.21 (Observe that the factor s here appears because the 
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general integrals (j303H and (|304p are denned without the appropriate prefactors present in the 
definitions of $( 3 ) and 

Let us stress that one can also apply the technique of MB representation in a similar way 
in various situations where a given four-point off-shell Feynman integral cannot be reduced to 
ladder integrals. 

B Two-loop calculation of the rectangular light-like Wilson loop 

In this appendix we present the detailed calculation j95[ EE] of the two-loop four-point Wilson 
loop. 

B.l Computation of individual diagrams 

Due to non-Abelian exponentiation, we only need the diagrams containing then non-Abelian 
colour factor ~ CfN c . The two-loop correction w^ 2 ' is given by a sum over the individual 
diagrams shown in Fig. [27] plus crossing symmetric diagrams. It is convenient to expand their 
contributions in powers of 1/e and separate the UV divergent and finite parts as follows 

^ (2) =E[(-^3A 2 ) 2e + (-^4A 2 ) 2e ] {^+^+^+^}+4°>+0( e ), (307) 

where the sum goes over the two-loop Feynman diagrams shown in Fig. [271(c)— (1). Here 
(with < n < 4) are dimensionless functions of the ratio of distances x\^jx\^. Making use of 
(|307p . we can parameterise the contribution of each individual diagram to the Wilson loop by 
the set of coefficient functions A_^. In (|307p fx is defined by jl 2 = /i 2 e 7E 7r, where 7e is the Euler 
constant. The gluon propagator in the Feynman gauge is 

D^(x) = rf v D{x) , D(x) = - * 3 (e' 7E A 2 ) £ (-x 2 + i0)~ 1+e . (308) 

For simplicity we will consider the Wilson loop in the region where x 2 3 > 0, x\^ > 0, which 
allows us to drop the +i0 prescription in the gluon propagator (|308|) . Let us also introduce some 
notations that will be used throughout this section: 

a = x 2 3 /(x 2 3 + X24) , a = 1 — a , 7 = a/a. (309) 

In addition, we will sometimes use the notation 

Pi =<+!-<■ (310) 

Finally, we will drop a factor of (— fi 2 ) 2 " in intermediate steps of the calculation, since it can be 
easily reinstated at the end. 

B.l.l Contributions independent of the kinematics 

In [85], the two-loop computation was carried out in the special kinematics x 2 3 = —x\^ = 1. 
Since diagrams /( c ),/( d ) and 7^ only depend in the kinematics through trivial prefactors like 
(— x 2 3 fl 2 ) 2e , we can use the results of |85j . and find the following contributions^! to the coefficient 

45 For diagram 1^ we take into account the difference between the DREG scheme used in [85] and DRED 
scheme used here, see [16j . 
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functions denned in (|307|) . 



Ac) 

A -4 


= 0, 


Ac) _ 1 Ac) _ 1 

-3 - o > A -2 ~ a i 
o Q 




1 

= 2 + 


7f2 A (c) 
4X ' 


vr 2 

= 2 + — 

1 9 


+ k 3 . 




(311) 


Ad) 

A -4 


1 

lb 




vr 2 

yo 




= --Cs 

24 


_ 

/i — 


7 ^ 

7T . 

nnon 


(312) 


AS) 

A -4 


1 

"16' 


,(/) _ 1 ,(/) _ 

A -3 _ g > A -2 — 


1 

"4 + 


5 ^ 

96 ' 


A<3 = - 


1 ^ 2 

2 48 + 




(313) 


Af) 
A o 


= -2 - 


12 2880 6 s 












(314) 



Let us now consider the diagrams containing two crossed gluon propagators, i.e. 1^ , I^ h ' , 1^' , 1^' . 
Their group theory factor is Tr (t a tbt a tb) = Cf (Cf — l-^c)> of which we keep the part — ^CfN c 
only because of non-Abelian exponentiation. 

B.1.2 Diagram I^ e) 

We consider the crossed diagram [pi,P4] [^4,^3] (the square brackets indicate the two segments 
to which the gluon is attached) shown in Fig. l27T e). The other diagrams of the same topology 
will be taken account of by including a symmetry factor 4 for the contribution to . The 
contribution of the diagram takes the form 

/( e ) = -^C F N c (ig) 4 [ dh [ dt 2 [ dst [ ds 2 

Jo Jti Jo Jo 

x (Pi • P4)(P3 ■PA)D{pis 2 +P4^)-D(P4(1 - h) +P3S1) ■ 

Taking into account the light-likeness of the edges of the Wilson loop are light-like, i.e. pf = 
x\ i+l = 0, we find 

= "^^F r2(1 " ^lA^y f dt ^ J dt * J ds ^ J ds * ( fi2 ^(l " h) Sl 

" 11 " ° (315) 
Doing the integrals, one finds that the contribution of the diagram is 



-1+e 



,<«> _ -S^r^i - c) (AAe-^r\ fi - ££±2 



2(87r 2 )2" r ~ e )( x i3 x 24 e y g4 



(316) 



It develops a double pole in e. Including the symmetry factor, its contribution to is 



A {e) 2 = --, A^ = ks, ^--Mi-iV. (317) 
-2 24' _1 2 12 720 v 7 



B.1.3 Diagram /CO 

The contribution of this diagram takes the form 



= -IjV c Cjr(^) 4 / dtl f dt2 f dsi f ds2 



/0 Jti Jo JO 
x (Pi ■Pi)(P2 ■ Pa)D{pis 2 + Pit 2 )D(pi +p 2 Si +PAt\) 



Integration over si,S2 and £2 yields 

/W = i A r cC 7 F ^^r 2 (l- e )(e- 2 ^x| 4 ) £ ( a ; 2 3 + x 2 4 ) 
1 Z" 1 dti(l-tf) 



6 JO x 13 fa 13 + x 24)^l 



[fa 2 3 (l-ti)) e -fa| 4 ti) e ] (318) 
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This expression in not symmetric under a — > 1 — a. Its symmetrisation gives (upon substitution 
x — > 1 — x) 



j(h)(sym) 



^ CFNc ii^{h 



[(4.) 2e + 0&) 2e ] [ — In (?) In (xx) 
J x — a Vox/ 



(319) 



+- 



1 



dx /(\.r- 
In 



f In 2 (xx) - \ In 2 (* 
2 2 \ a 



iq x — a \ax 
and in the notation of the basic integrals of section IB, 31 

2 



1, , /ax 
- In — 

2 Vax 



+ 0(e) 



j(/i)(sym) 



3,, 1 , o /a\ , , 1 



+ 



M 3 -hn 2 (-)m 1 + =-M 9 
2 2 Va/ 2 



+ 0(e] 



Including the combinatorial factor corresponding to this diagram, c^ 1 ) = 8, we have 



= \m 2 



A 



CO 



7T- 



; Mi + -M 3 H Mi . 



(320) 



(321) 



B.1.4 Diagram JW 

The contribution of this diagram takes the form 



/'" -iiV c C F (ig) 4 / dt x [ dt 2 [ ds 2 [ dst 

JO Jti JO J s 2 

x (P2 • Vi) 2 D(pi + p 2 si + Pih)D(p 1 + p 2 .S2 + P4*2) 

It is finite for e — > and, therefore, the calculation can be performed in D = 4 dimensions 



rW — _I /v (7t7 - / 



with 



J = / dti / dt 2 / ds 2 / dsi [(asi + ati - siti)(as 2 + at 2 - s 2 t 2 )\ 1 

J Jtl J ./S2 

Then, after integrating over s 2 and ii and rescaling s\ — ► asi and i 2 — ► od,i one gets 

/'/'■ /■'/» d Sl dt 2 ( i-t 2 \ ( ,.i-»i 

It seems natural to simplify the integrand by changing integration variables 



x = S\- 



1 ~*2 
*2 



y = *2 : 



1 - s 



which leads to 
J = 



dxdyf{x,y), f( x ,y) = —^—^—— 
E x(l+x)y(l + y) 



ln(l + x) ln(l + y) 



where the integration domain S is defined by 



1 — xy 1 

< < 1 + - 

1 + y 7 



1 — xy 

0< -— - < 1 + 7. 

1 + x 



a s 
7 = ^ = t>0. 

a r 



(322) 



(323) 



(324) 



(325) 



(326) 
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This can be written more conveniently as 

1 



V < - 

x 



(x > 0) 

7^ 



V < -(1 + 7 + -) 

x 

1 

y > 



(-Kx< 0) 



1 + 7 + 72 

Thus the integration over £ can be parametrised as follows: 



(327) 
(328) 
(329) 



J 



dx 



l/x 



dyf(x,y) + 



dx 



dyf(x,y) + 



1+7+7^ 



dx 



-(1+7+7/2:) 



dyf(x,y) . 



1+7+7^ 



Ji 



J2 



J34 



(330) 

The integration regions lie in the first, second and third and fourth quadrant (clockwise). The 
first integration can be easily done using 

1 



y hi + x) 

dzf(x,z) = 

X(l + X) 



(1 + x) ln^(l + y) - Li 2 (-y) 



(331) 



The final integration is harder, but all integrals can be expressed in terms of poly logarithms 
[22 lj. After some algebra we find 

(332) 



Ji 
■h 



7T 

180 

1Q7T 4 

+ — In* (7) + — ln 4 ( 7 ) + — ln^l + 7) + - ln( 7 ) ln 3 (l + 7) 
360 12 v ,; 24 KU 12 v ; 6 y n v ' 



7T 9/ s 1 , a, . 7r , 9/ , 1 , . , , 1 , a 



12 



ln 4 (l + 7) 



*~ 1. 2 , , 1 



_--ln^( 7 ) + -Li 2 (- 7 ) 



In 



1 + 7 



[Li 3 (-7) + Cs] -Li 4 



Li 2 (-7) 
1 



1 + 7 



Li 4 



1 + 7 



(333) 



J: 



31 



^ - ^ ln a (l + 7) - 1 111(7) ln 3 (l + 7) + 7^ ln 4 (l + 7) 



+ 



+ ln 



ir 1. 2 , , 1 



— --lr/( 7 ) + -Li 2 (- 7 ) 



7 



1 + 7 



[Li3(-7) + Cs]+Li 4 



Li 2 (-7) 
1 



1 + 7 



+ Li 4 



7 



1 + 7 



(334) 



Here we used frequently polylogarithm identities like those given in appendix IB.4I in order to 
simplify the individual expressions. Summing them up, we see that most terms cancel between 
J 2 and J34, and we obtain the simple result 



J=^[vr 2 + ln 2 ( 7 )] : 



Including the symmetry factor 2, we find 



a« = __L Ml 2 
24 1 



(335) 
(336) 



B.1.5 Factorised diagram JW') 

The contribution of this diagram takes the form 



jCj) 



±C F iY c (i<7) z 



X (P2 -P4)(Pl -P3)D(p 2 +P1S1 +p 3 ti)D(pi +P2S2+PAt 2 ) 
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and the integral factorises into the product of finite one-loop integrals, 



X 24 X + X 13y ( X 13 + x 2a) X V- 



(337) 



Evaluating the integral we find 

1 



B.1.6 Self-energy diagram 1^ 

The contribution of this diagram reads 



dxdy (xj 3 + x| 4 ) 

I 

= (338) 



/(*) = -g 2 C F (p 2 ■ p 4 ) f dsx f dt 1 DW( Pl + P2 s 1 +p 4 t 1 ), 
Jo Jo 



(339) 



where the gluon propagator in the Feynman gauge, with the one-loop correction included, is 

(340) 

with D = 4 — 2e. Here we added the contribution both of e— scalars in the DRED scheme and 
of n s scalars. Then, 

' ik) = 2~MirD C ' F l (3D - 2 - - 2( ° - 2) "' - "J ( D -4KD- 3 ")(P-1) JW ' < 341) 

with 

j( fe ) = (x^ + x 2j /" efei /" [_ x 2 3Sl _ ^ + (x 2 3 + x l A ) Sltl f-D (342 ) 

JO j0 

The integration is straightforward. In N = 4 SYM, for nf = 4N C and n s = GN C , 



1 / 2 \ 2 

/(fe) = 8\h) C FN c {e- l [{xl^ + {xl A )^]M l +AM l + 2M 2 + 0{e)} 
The combinatorial factor corresponding to this diagram is = 2 , so we have 



(343) 



A™ = -M x , Af =M l+ X -M 2 . (344) 
B.1.7 'Mercedes-Benz' diagram 1^ 

We assume that two gluon legs are attached to the segment p 2 and the third leg is attached to 
the segment p^ 

= \g A C F N c f 1 dh I' dt 2 f 1 dt 3 (345) 
^ Jo Jo Jt 2 

xp^p^ 2 p^ 3 T^ W2fM3 (-id Zl ,-id Z2 ,-id Z3 ) / d D z A \\D(zi - z A ) , 

J i=i 

where 

z i = ~Pi ~ Pih , 4 = P2 f 2 , 4 = P&s , (346) 
and the three-gluon vertex is given by 

P2 2 P2 3r W 2 w(Hd Zl , -id Z2 ,-id Z3 ) = i(p 2 ■ Pi){p 2 d Z2 - p 2 d Z3 ) = i(p 2 ■ Pi){d t2 - d t3 ) . (347) 



Then, 



I® = ^g 4 C F N( P2 ■ PA ) f dh f dt 2 f dt 3 (d t2 - d t3 )J{z 1: z 2 , z 3 ) , (348) 
1 Jo Jo Jt 2 

where the three-point integral J(z\, z?, z 3 ) is defined in appendix IB. 21 By virtue of z| 3 = 



^1,-22,23) 



32tt Z) 4 - D 



dr (rf ) D/2 " 2 [(-Z! + rz 2 + fz 3 ) 



(349) 



Using the explicit expression for the z— coordinates in (I346p and equation (j349H we find 



j(0 



. 2 _ 2D5 4 C F iVr(D-3) /- 1 cft x (-* 2 3 - x| 4 ) f 1 



1287T 13 (4 - D) 2 J x? 3 ti - a&ti 



'24 ''I JO 



dr(rr 



n/2-2 



1 



/ 2 r 
(s 13 ti) 



! 21 



4-D 



+ 



1 1 



r r 



, 4--D 



Notice that in the second term, one can integrate by parts using 

{Tf) D/2-2(l _i\ {D/2 _ 2 yld D /2 _ 2 

\t t J dr 

After some algebra one arrives at 

g*C F N r 2 (l- £ ) o _ 2e7E 



(350) 



j(0 



64^ 4 (2e) 3 (l-2e) 



dt 
a — t 



2c 



2 \2e,2c 



r(2 -2e) , , 



I*(l-e) 



where we replaced D by 4 — 2e. Expanding in powers of e we find 

/ (0 = -^?{ £ [(40* + O^)*] (2M a + M 2 ) + (l - 1 . 2 ) Mx + ^M 2 + I M 3 + 0(e) } , 

(351) 

where we used the basic integrals of section IB. 31 The combinatorial factor corresponding to this 
diagram is = 4 , so we have 



A 



(0 



1 



-Mi - \m 2 , 



A, 



(0 



7T 1 1 

-Mi H Mi M 2 M 3 

24 2 8 



(352) 



B.1.8 T-shaped diagram 



We assume that two gluon legs are attached to the segment p 2 and the third leg is attached to 
the segment p 4 



/(<?) 



^g 4 c F N c I ,7/, I dl. 



efts 



(353) 



3 

Pi 2 P 2 3r Mi^2w (-idzi ,~id Z2 ,-id Z3 ) / d D z A ]J - z 4 

i=l 



where 

Zl = -p 4 t 1 , Z2 = Pl t 2 , 2 3 = + p 2 t 3 
and three-gluon vertex is given by 



(354) 



1 3 — Pi P 2 1 MIA 

= • P4)(P2 • {d Z2 ~ dzi)) + i(Pi ■ P2){pa • (d Z3 - d Z2 )) + i(p 2 ■ Pa){pi • {d Zl - d Z3 )) . 
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Translation invariance of the ^—integral can be used to simplify the vertex to 



F3 = -i{prPA)dt 3 -i{prP2)d tl -2i(pi-p 4 )(p2 ■ d Zl )+2i(p 1 -p 2 )(p 4 ■ d Z3 )+i(p 2 -p 4 )(pi ■ (d Zl - d Zz )) 

(355) 

Let us rewrite the integral as 1^ = l[ 9 ^ + , where 
M ; " " " 



r(fl) 



-g c F N c / dh / dt 2 / dt 3 {2X 2A d t3 + 2 -x l3 d tl ]V(z 1 ,z 2 ,z 3 ) 
1 Jo Jo Jo 



g 4 c F N c [ dh [ dt 2 f dt 3 
Jo Jo Jo 



x [-^24((P2<9 Z1 ) +d t3 ) + x 2 13 ((p 4 ■ d Z3 ) - d tl ) - \{x 2 13 + x 2 2i )(pi ■ (d Zl - d Z3 ))] V(z 1 ,z 2 ,z 3 ) 
The first contribution l[ 9 ^ can be split into a divergent and a finite part, l[ 9 ^ = l[ 9 J + l[ 9 j) 



T (g) 
l id 



j(9) 



= t;9 4 cfN c dt 2 - \x\ 3 I dt 3 V(0,z 2 ,z 3 ) - \x\ A I dt 1 V(z 1 , z 2 ,Pi) 

* Jo l Jo Jo 

= ^-g 4 c F N c dt 2 +^xf 3 dt 3 V(-p 4 ,z 2 ,z 3 ) + \x\± \ dt 1 V(z 1 ,z 2 ,p 1 + p 2 ) 

2 Jo L Jo Jo 



Divergences come from the two integrals in the first line only since the three arguments of the 
V— functions in second line cannot be light-like. The divergent integrals can be easily evaluated 
as 



(xj 3 ) 2e + {x 2 2A f e 



1 



r(l + e)r(l -2e) 

r(i-e) 



(356) 



Next, we examine the finite part of I\. it can be written as 



x 



21 



'fo(p2t2 + Plh)x 2 13 + PlhP 2 t 2 X 2 4 

Here [d/3] 3 = d(3\d(3 2 d(3 3 8(l — [3\ — (3 2 — [3 3 ). After changing variables, we have 

T { 9 ) _ g 4 C F N c {g) 
J i/ " 4-64tt 4 l f ' 



+ [x\ 4 Uxl 3 ]. (357) 



where 



J 



1 r-1 r-1 r-1 

I dx j d\ I dt 2 I dt\ 



(g) 



o Jo X(xt 2 + xii)7 + \xxt\t 2 



+ (7 1/7) 



(358) 
(359) 



After some algebra one finds 

r 2 



J 



(<?) 
1/ 



^vr 4 + ^ln 2 ( 7 ) + lln 4 ( 7 ) - ^vr 2 In ( 7 ) In (1 + 7 ) - ^ln 3 ( 7 ) In (1 + 7 ) 
+ Vln 2 (1 + 7) + \\n 2 ( 7 )ln 2 (1 + 7) - ^n 4 (1 + 7) + y Li 2 (-7) 



+2 In (7) Li 3 
+ (7 1/7) • 

r(ff) 



7 



1 + 7 



2 In (7) Cs + 2Li 4 (-7) - 2Li 4 



1 



1 + 7 



2Li 4 



1 + 7 



Let us now study I 2 y> . For D = 4 — 2e it is given by 



r(9) 



9 A C F N C (g) 
"4-64tt 4 2 ' 



J 2 (9) = -(1 - 2e)r(l - 2e) rfti C dt 2 C dt 3 f ~ [dpMfafofoY 

Jo Jo Jo Jo 



{x\ 3 + x 2 4 ) [PM(3 2 - 2(3 3 )x 2 24 + (3 3 t 3 ((3 2 - 2f3 l )x 



[Mf3 2 i 2 + /3it"i)x 2 3 + PMfch + (3 3 i 3 )x 2 2A 



2 l2-2e 
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Then, we change the integration variables according to 

(3i = Xx, = X , (3 3 = Xx , 
and integrate over t 2 , which leads to 



J, 



GO 



I' 1 I' 1 

(x 2 3 Hl + 7 ) / dhdh / dXX 



Jo xtu - xt 3 



X < X 



-l+2e 



[xt 3 Ati + £i(A + Axi 3 )7]' 



-l+2e 



-1+2<e 



[t 3 (A + Axti) +xtiXt 3 -f\ 



-l+2e 



We notice that integration over small x and x produces a single pole in e. We split up j!f^ into 
a simplified divergent integral and several convergent integrals, in which we neglect 0(e) terms. 



J. 



GO 

2d 



(x 2 3 ) 2e (l + 7 ) I dhdh I dXX~ e I 
Jo Jo Jo 

,-l+er+_ \x , + ,\ , \z\„/|-l+2e 



x - 1+£ [t 3 Ati + ti(A + At 3 )7]~ i+ " £ - x~ 1+e [t 3 (A + Xh) + h\i 3 j] 



l-l+2e 



J, 



GO 

2/ 



(1+7) / rfi lC Z£ 3 f dX f 
Jo Jo Jo 



dx 



- 2 

+ 2 

1 

x 
1 



A t 3 + ti7 
A xii7 — xt 3 
1 



x 



xii7 ~~ 



xt 3 Aii + ti(A + Axt 3 )7 t 3 (A + Xxt\) + xtiAi 3 7 

*17 *3 



xt 3 Ati + £i(A + Axi 3 )7 t 3 (A + Axii) + xtiAi 3 7 
1 1 



xi 3 Ati + ti(A + Axt 3 ) 7 t 3 A*i + *i(AAt 3 )7. 
1 1 



z L £ 3(A + Xxti) + xtiXt 3 j t 3 (X + Xti) + tiXt 3 j 



To summarize, 



4 9) = 4i + 4f + O(e) , 



where only the first term produces poles in e. Calculating J$ we get 
rGO 1 



J, 



2</ 



— X 



13 



2 \ 2e / 1 dx , xa , 
+ x 24 J / m — m(xx) 



dx 



x — a 



q x — a xa 

ln(xx) (ln 2 (xa) — ln 2 (xa)) + In 3— I h lnxlnx 

xa \ 6 



(360) 



(361) 



or equivalently, writing j!fj = j!fj d + J^l , 



J. 



GO 



J, 



2dd 
(9) 



2c 



i(xi 3 r + (xi,r] m 2 



(362) 



2df 



ln(aa)M 2 



dx 
x — a 



ln(xx) (in 2 



ln 2 (xa)) + In I h lnxlnx )(3j63) 

xa \ 6 / 



xa 



where the integral M2 is defined in (|377p . From (|356p and (|362p we can read off the contribution 
of 7^0 to the pole coefficients (taking into account the combinatorial factor 4 for this diagram), 



A ( 9 ) _ ZL 

- 2 - 48 ' - 1 



1 1 
--M 2 + -C 3 . 



(364) 
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Let us now turn to the finite part. All integrals appearing in it can be evaluated in terms of 
poly logarithms. After a lot of algebra, one finds 

4$ = {^vr 4 + ^ln 2 ( 7 ) + lln 4 ( 7 )-^ 2 ln( 7 )ln(l+ 7 )-iln 3 ( 7 )ln(l + 7 ) 

+2vr 2 ln 2 (1 + 7) + In 2 ( 7 )ln 2 (1 + 7) - ^ln 4 (1 + 7) + 3 In (7) Li 3 (-7) 

6 In ( 7 ) Li 3 (^J\ - 3 In ( 7 ) Cs - 6Li 4 - 6Li 4 (^-) } (365) 



4f = |^ 2 + in 2 ( 7 ) + ^ln 4 ( 7 )-^ 2 ln( 7 )ln(l + 7 )-iln 3 ( 7 )ln(l+ 7 ) 



and 



2 f ~ \120 ^2 w; ^8 wy 2" U J 2 

+2vr 2 ln 2 (1 + 7) + in 2 ( 7 )ln 2 (1 + 7) - ^ln 4 (1 + 7) + 3 In (7) Li 3 (-7) 

+6 In ( 7 ) Li 3 - 3 In ( 7 ) (s ~ 6Li 4 (j^) " 6Li 4 (y^j } (366) 

Summing up (I360p . (l365p and (j366H we find 

" 4f + 4 9 J f + 4f = + Hvr 4 - M 3 . (367) 

Thus, taking into account the symmetry factor 4 and a constant contribution coming from 
equation (|362p . we find the total contribution to the finite part of this diagram to be 

j&) = L M l - -M 3 + — vr 4 . (368) 
24 1 4 360 v ' 

This completes the calculation. Summing up all contributions, we find formula (|307p . 
B.2 Useful formulae for diagrams with three-gluon vertex 

In configuration space, these diagrams involve three propagators joined at the same point z 
which is integrated out 

V(z 1 ,z 2 ,z 3 ) = ^~ 2e [ d D zD(z- Zl )D(z-z 2 )D(z-z 3 ), (369) 



with D = 4 — 2e and D(z) being the gluon propagator in the Feynman gauge 

r 1+e 

Shifting variable z — > z + z\ and performing the integration one gets 

,-2e f-oo 



D(z) = - l — r (Tvn 2 Y I dss- e e- isz \ (370) 







V(zi,z 2 , z 3 ) = ^n;(vr// 2 ) 2e / dsids 2 ds 3 - 2 3 — r^— exp (iA(si,s 2 , s 3 )) , (371) 

647T 4 J (si + s 2 + s 3 ) z 6 



with 



a, n ( S 2^21 +S 3 Z 31 ) 2 2 2 SiS 2 z\ 2 + S 2 S 3 Z% 3 + 838^ 

A(si,s 2 , s 3 ) = ■ ■ s 2 z 2l - s 3 z 3l = ■ f 2_ . (372 

Sl+S 2 + S 3 Si + s 2 + s 3 



Changing integration variables according to s, = A/3j with ^ /3j = 1 and < < 1, one obtains 
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Let us consider the integral in the kinematics (z 2 — Z3) 2 = 0, where it can be further simplified. 
In that case, 



^(si,s 2 ,s 3 ) 



si(s 2 + s 3 ) 



si + s 2 + s 3 
Then, we introduce standard parametrization 



S 2 Z 2 + S3Z3 
S2 + S3 



pfii, /?i = 1 - A , [5 2 = Ax , /3 2 = A(l - x) 



and find 



■l+4e 



64vr 4 



2 . 2e r(l-2e) /- 1 



dx (xx) 



[(X«21 +XZ 31 ) 2 ] 1 -2 £ 



(374) 

(375) 
(376) 



B.3 Basic integrals 

Here we give explicit expressions for integrals encountered in the two-loop calculation. 



Mi 

M 2 

M3 



dx . rax 



In 



x — a V ax 
1 dx 



1 r 
2 



vr z + ln z ( - 

a 



In (™) ln( 



x — a V ax / 



xx 



7r /a 
-— ln(aa) + 2Li 3 (l - Li 3 -- 
2 V a 



Li 3 



O: 



a 



ln^-J [Li 2 (a)-Li 2 (a)] 



dx /ax 
r m — 

x — a V ax 



ln 2 (xx) 



(377) 
(378) 

(379) 



.il^ _ .^2 ^2 (a) + 61n(a) ln((5) + ln 2 (s) ] 



-— [ln 4 (a) +ln 4 (a) +41n(a)ln 3 (a) - 181n 2 (a) ln 2 (a) + 41n(a) ln 3 (a)] 



-4 In (~J [Li 3 (a) - Li 3 (a)] + 8 [Li 4 (a) + Li 4 (a) 



Ma 



1 dx , o /ax 
In 13 l — 



x — a Vax/ 

where a = 1 — a and Li n (z) (with n = 2,3,4) are poly logarithms [221J. 



(380) 



B.4 Identities for poly logarithms of related arguments 

Most (or all) of the following polylogarithm identities can be found in Lewin's book [221] . They 
can be easily proved by differentiating (the integration constant can be fixed by evaluating the 
identity for a particular value of z). 



Li 2 i—z) + Li 2 



1 



7T 



7T 



Li 2 (z) + Li 2 (l-z) = -ln(z)ln(l - z) + — 

o 

It follows from (|382[) . by putting z —> 1/(1 + z), 



Lio 



1 



+ Li 2 



.1 + 2 

Similar useful identities are 



1 + z 



In 



l + z 



In 



1 + z 



Li 2 (1 - z) + Li 2 1 



dn 2 (z) 



z>0. 



< z < 1. 



7T 



(381) 
(382) 



+ — , 2 > . (383) 




z > 0. 



(384) 
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and 



Li 2 (-z) + Li 2 



l + z 



--In 2 (l + z), z>0. 



(385) 



Using (|381j) . (|383p and (I385p the functions Li 2 (—5) , L12 (iTij 1^2 \TTz) can au ' 3e ^expressed 
in terms of Li 2 {—z) and logarithms. The inversion identity for Li 3 reads 



There is also 



Li 3 (-z) - Li 3 



Li 3 (-z) + Li 3 



~hi(z)-hn 3 (z), z>0. 
6 6 



(386) 



l + z 



+ Li 3 



l + z 



C3 



[ ! 7T 2 

-ln 3 (l + z) - -ln(z)ln 2 (l + z) - — ln(l + z) 

o Z D 



For Li4 there is just the inversion identity: 



1 



7T" 



7vr 4 



Li 4 (-z) + Li 4 — = ln 2 (z) ln 4 (z) , 

1 z ) 12 v 7 24 w 360 



z > 0. 



z > 0. 



(387) 



(388) 
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